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ACTIVATION ENERGY FOR CREEP TIN! 


ABSTRACT 


The literature values ranging from 8000 26,000 cal per g-atom for 
the activation energy for the creep tin. The present investigation analyzed 
the results constant stress creep tests during steady state creep, and 
means temperature cycling. value approximately 9000 cal per g-atom for 
the creep high purity tin temperature range 300 350° was obtained. 


INTRODUCTION 


The activation energy concept applied creep has received considerable 
attention the past few years. The activation energies reported for creep 
tin vary widely the temperature range 300 500° Thus, below 360° 
values from 8000 per g-atom are reported. However, above 
approximately 360° all investigators agree activation energy the 
vicinity 21,000 cal per g-atom. Table shows the values which have been 
reported, and also the values obtained this investigation. 

Dorn and co-workers (1957, 1955) claim that the activation energy below 
approximately 369° the same the activation energy above this tem- 
perature. Breen and Weertman (1955, 1956, 1957), however, have shown 
clearly their investigations that over the temperature range 298° 
the activation energy 11,000 cal g-atom. This value agrees quite 
well with the analysis the data Tyte (1938) which was conducted 
Breen and Weertman, and also the authors. Tyte’s data value 
about 8990 cal g-atom. 

The present investigation, therefore, was conducted the temperature 
range 300° 350° 


EXPERIMENTAL PROCEDURE 

Spectrographically pure tin (99.998%) and zone-refined tin 
estimated) were extruded into diameter wire and tested using Andrade- 
Chalmers (1932) constant stress, beam-type creep machine. Thermostated 
water baths maintained temperature within +0.5°C. Extensions were 
measured directly with cathetometers accurate 0.01 and 0.001 cm. 


1Manuscript received June 10, 1958. 
Contribution from the Department Metallurgical Engineering, University Toronto, 
Toronto, Ontario. 
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RESULTS 
After the samples had entered the secondary stage creep, the creep rate 
was determined. The logarithm the creep rate was then plotted against the 
reciprocal the absolute temperature. The slope the line thus yielded 
the activation energy divided the gas constant (Fig. 1). 


-1.6r 


Spectrographically Pure Tin 


350°K 


Fic. energy plot for creep tin. 


The method temperature cycling developed Frenkel, Sherby, and 
Dorn (1955) was also employed. typical creep curve resulting from this 
procedure shown Fig. The activation energies could then calculated 
according the formula 


activation energy, 
creep rate before change temperature, 
creep rate after change temperature, 
temperature before change, 
temperature after change, 
gas constant. 


Only one run was completed successfully the zone-refined tin owing 
difficulty with excessive grain growth which occurred room temperature 
after extrusion. 

The experimental results are recorded Table 
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Run 
Pre-load psi 
Stress 605 psi 
Spectrographically Pure Tin (99.998 Sn) 
Temperature cycled shown 
4 
nse ns 0 6 6&6 CO CO WO NO 1S 1S 190 
Time, hours 
Fic. Typical creep curve determined the temperature cycling method. 
TABLE 
ACTIVATION ENERGY FOR THE CREEP TIN 
Activation 
energy, Temperature, 
cal/g-atom Reference Remarks 
23,300 320-336 tests, 99.4% Sn. Polycrystals 
22,700 351-394 tests, 99.4% Sn. Single crystals 
21,000 298-478 
19,000 Grain boundary shear. Analysis 
and 
11,000 99.9% Sn. Break curve 
26,000 99.9% Sn. Break curve 
8,000 292-343 Analysis authors and Break 
for 
17,400 
19,000 298-373 Quoted Grain boundary relaxa- 
tion 
11,000 298-403 Single crystals. Break curve. 99.99% 
22,000 403-490 Single crystals. Break curve. 99.99% 
12,000 298-393 Single crystals compression. Break 
curve. 99.99% 
24,500 Single crystals compression. Break 
curve. 99.99% 
9,600 (This 300-350 Sn. Log vs. 1/T 
investigation) 


11,000+3,000 


9,000 
1,500 


Sn. Cycling. Average 
determinations 

Zone-refined Sn. Log vs. 1/T 

Zone-refined Sn. Cycling. Average 
determinations 
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DISCUSSION 


The four activation energies determined this investigation are satis- 
factory agreement. The slightly higher mean value for the cyclic temperature 
tests pure tin can partially accounted for the data Tyte 
(1938) and Breen and Weertman (1955), shown Figs. and are 
considered. break the vs. 1/T curves obvious each figure. This 
break occurs over temperature range with the temperature the 
break appearing stress dependent. Figure this investigation also 


Increasing 


Stress 
400° K 350° 300° 


1 8 


Fic. energy for creep tin after Tyte (1938). 


shows slight upward curvature one end. This upward curvature could 
explain the high mean value 11,000 cal g-atom obtained for the spectro- 
graphically pure tin when the cyclic temperature method used obtain 
the activation energy. The investigations Breen and Weertman (1956) and 
Weertman (1957) also show break the vs. 1/T curves. 

The activation energy for creep not measurably affected minor 
impurities, which accord with the work Sherby and Dorn (1952). 


CONCLUSIONS 


The activation energy for the creep high purity tin the temperature 
range 300° 350° approximately 9000 cal g-atom. 
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Increasing 
Stress 


Fic. energy for creep tin after Breen and Weertman (1955). 
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THE ELECTRON NEUTRINO ANGULAR CORRELATION 
THE BETA DECAY THE FREE 


ABSTRACT 


The electron neutrino angular correlation coefficient has been measured for 
the beta decay the free neutron measuring the momentum spectrum the 
electrons emitted restricted angular range from the direction the recoil 
protons. Experiments were made for two different angular ranges and yielded 
best value for the angular correlation coefficient +0.07 with error 
+0.12. This error includes the statistical error the form its standard 
deviation and allowance for all known systematic errors. 


INTRODUCTION 


Before the discovery that parity not conserved the beta-decay process, 
measurements the electron neutrino angular correlation were essentially the 
only experimental approach our understanding the form the inter- 
action occurring the beta-decay process. The discovery parity non- 
conservation has provided new experimental approaches this problem but 
measurements the electron neutrino angular correlation still provide the 
most widely applicable method determining the interaction. This paper 
describes measurement the angular correlation coefficient for the beta 
decay the free neutron. Preliminary results were reported the 1955 
meeting the Royal Society Canada Toronto, Ontario, and abbre- 
viated form were reported elsewhere (Robson 1955). Since those reports, 
additional data have been accumulated and more careful investigation has 
been made the systematic errors inherent the measurements. The final 
result including all data very similar the value previously quoted though 
the fuller appreciation the systematic errors has resulted the final error 
being slightly larger. 

Since the neutrino cannot readily detected necessary infer the 
electron neutrino angular correlation from measurements based the detection 
the electron the recoil nucleus. There are several possible methods 
including direct measurement the momentum distribution the recoil 
nuclei 1957; Good and Lauer 1957; Maxson al. 1955), 
the angular distribution between the electrons and recoil nuclei (Rustad and 
Ruby 1955) and the momentum distribution the electrons within limited 
angular range from the direction the recoil nuclei (Alford and Hamilton 
1957). All these measurements are dependent the electron neutrino angular 
correlation coefficient and therefore give estimates its value. The present 
experiment was based the last the above-mentioned techniques. This 


received July 1958. 
Contribution from Physics Division, Atomic Energy Canada Limited, Chalk River, 
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Issued A.E.C.L. No. 690. 
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method was chosen since afforded the best method reducing the effects 
unavoidable background radiation. 

The full theoretical expressions for the momentum distribution the 
electrons from the beta decay the neutron have been given Jackson, 
Treiman, and Wyld (1957) for parity non-conservation and possible time 
reversal violation. Since this experiment unpolarized neutrons were used, 
great simplifications occur the formula and further simplify using 
the notation represent the full expressions for the coupling 
constants corresponding parity non-violation and parity violation the 
ith interaction, the formula becomes 


where the ratio the electron velocity the velocity light, 
the angle between the electron and neutrino, 
the element solid angle between the electron and neutrino, 
the electron neutrino angular correlation coefficient. 


a= 


and would equal for pure scalar interaction, 
for pure vector interaction, 
for pure tensor interaction, 
for pure axial vector interaction. 


this expression for the momentum distribution, the so-called Fierz 
interference term has been neglected. Experiments the ratio capture 
positron emission the decay (Sherr and Miller 1954), and the 
shape allowed spectra indicate that this interference term small. The 
effect its possible presence will discussed the discussion the end 
this paper. 

The corresponding expression for the momentum distribution emitted 
fixed angle from the proton recoils can readily obtained applying the 
conservation linear momentum between the electron, the neutrino, and 
the proton. becomes: 


sing 


where fo(p) and are now functions the electron momentum and the 
angle Full expressions for and are given the Appendix the 
end this paper. 

the apparatus detects events which have angle between the directions 
the proton and the electron with efficiency the observed momentum 
spectrum the electrons coincidence with the protons will be: 
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where and are functions the electron momentum and are given 
the Appendix. F(p) the Coulomb correction term for 

Thus, provided g(¢) not independent measurement the momentum 
distribution the electrons coincidence with the protons can compared 
with the above expression for this distribution and value derived 
the least squares method fitting the experimental points the theoretical 
curve. Two experiments were carried out with different experimental arrange- 
ments giving different functions for the first, referred hereafter 
series was such that coincidences were recorded for between 
approximately 150° and 170°. the second, referred hereafter series II, 
was such that the coincidences were recorded for between approximately 
165° and 180°. The actual variations with which resulted from 
these two experimental arrangements will discussed later. 

Section this paper will describe the experimental equipment used for 
these measurements with general description followed detailed 
descriptions particular aspects the apparatus Sections 2(f). 
Section describes the experimental procedure, Section the theoretical 
spectrum shapes, and Section the experimental results. These are then 
discussed Section 


DESCRIPTION THE APPARATUS 
2(a). GENERAL DESCRIPTION 
The source neutrons was collimated beam from the NRX reactor 
Section The collimated beam entered large vacuum tank, passed 
through center electrode assembly surrounding the part the beam from 
which decay events were recorded, and then proceeded finally stopped 


1 FOOT 


BEAM CATCHER 


NRX 
REACTING CORE 


LEAD Z STEEL 
RY WAX MIX BISMUTH 
CONCRETE GRAPHITE 


Plan view the collimation arrangements used obtain the beam neutrons 
from the NRX reactor, the experimental apparatus, and the beam catcher. 
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beam catcher. The details the vacuum tank and its associated apparatus 
can seen from Fig. which shows plan view the equipment used 
for the series runs. The apparatus used for series differed only the 
arrangement the baffles the beta spectrometer and has been described 
elsewhere (Robson 1955). 

The apparatus the left side Fig. lens type beta spectrometer 
which selected the beta particles the required momentum. also defined 
the region the beam from which decay electrons could recorded and 
also the angles which these electrons had subtend its axis. defined 
the size the effective source neutrons and played dominant role 
determining the form the function 

The center electrode assembly surrounding the effective source was held 
potential relative the grounded vacuum tank. was equipped 
with guard rings maintain uniform potential inside the assembly that 
the protons resulting from neutron decay events the effective source recoiled 
essentially field free region. Those which recoiled with angle lying 
within about +10° the projection the spectrometer axis were able 
escape from the center electrode assembly through gridded aperture in. 
diameter and were then accelerated and focused onto the first surface 
electron multiplier. Thus, electron was recorded the counter 
the beta spectrometer and the corresponding proton the electron multiplier, 
the geometry the apparatus ensured that the decay event had occurred 
that part the neutron beam located near the center the center electrode 
assembly and that the angle between the electron and proton lay the 
restricted range approximately 150° 170° for series and 165° 180° 
for series more detailed discussion these angles will given Section 

The main advantage this experimental approach that the two counters 
are located considerable distance away from the neutron beam and are 
therefore less affected the high gamma radiation intensity which accom- 
panies collimated neutron beams. 


NEUTRON BEAM 

The neutron beam was obtained with the aid series collimators 
located horizontal beam hole the NRX reactor shown Fig. The 
hole through the reactor shield about diameter and equipped 
with large iron stopcock known the inner gate, which, when closed, 
placed total thickness in. iron the path the beam. This reduced 
the intensity the beam level which was possible work the 
apparatus without radiation hazard. Behind this gate bismuth plug was 
located the extreme inner end the beam hole adjacent the reactor 
core. The plug was in. thick with in. thick cylindrical skirt. attenuated 
the direct reactor gamma rays large factor and resulted reduction 
the number gamma rays the final collimated beam factor 
about three. Its presence only reduced the thermal neutron intensity the 
beam about 15%, presumably replacing those neutrons from the core 
which are scattered out the beam with neutrons scattered from the reflector 
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into the beam. Its effectiveness reducing the fast neutrons the beam 
not known though presumably considerably reduces the very high energy 
neutrons inelastic scattering. 

Outside the inner gate series collimators tapered diameter 
the entrance the vacuum tank. These collimators looked 
diameter the bismuth plug. Various materials were used the 
collimators changing from ordinary concrete the innermost region 
boron lead oxide wax mixture the intermediate rezion and finally 
unlined boron carbide casting the outer region near the experimental 
apparatus. the point where the collimators emerge from the main reactor 
shield slot was located into which boron carbide shutter could inserted 
shut off the thermal neutrons the beam without appreciably attenuating 
the fast neutrons gamma rays. This slot shown Figs. and 

The measured composition the beam the effective source region 
listed Table 


TABLE 
COMPOSITION THE COLLIMATED BEAM 


= — 4 
Beam intensity 
Component Inner gate open Inner gate closed 
Neutrons above 200 kev 
Gamma radiation r/hr 


THE BETA SPECTROMETER 

The spectrometer was the lens type using two thin electromagnets with 
their fields the same direction along the axis. Its function was threefold: 
(i) define the region the neutron beam from which decay events could 
detected, (ii) accept electrons whose direction lay only within defined 
limits, and (iii) select electrons with known momentum. The first two 
requirements were met using two annular ring entrance apertures which 
geometrically defined the region from which electrons could enter the spectro- 
meter and the directions which they must travelling. The last requirement 
was met the normal focusing properties spectrometer. The first the 
entrance apertures was located the central electrode assembly 
second about inches away from the neutron beam. 

The spectrometer was operated with poor resolution was possible 
subject the requirements that its behavior had accordance with 
the theoretical prediction and that could accurately calibrated. The 
poor resolution was required obtain high counting rate possible 
and the accurate calibration was required avoid the introduction serious 
errors the least squares fitting the experimentally observed spectrum 
the theoretical shape. These conditions were best met using ring focusing 
series and point focusing series II. minimize the possibility 
detecting scattered electrons incorrect momentum, the vacuum tank was 
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considerably larger radius than the maximum radius electron trajectory 
and the detector was arranged respond only electrons approximately 
the correct momentum focused the spectrometer. 

The electron detector was scintillation counter using anthracene mounted 
tapered lucite light pipe, oil-sealed selected DuMont type 6292 
photomultiplier. The photomultiplier was selected for low dependence the 
gain with counting rate and for resolution this particular arrangement. 
was shielded from the spectrometer magnetic field double iron-pipe 
magnetic-shield assembly. The pulses from the photomultiplier were passed 
through single channel pulse height selector before being used for coinci- 
dences with the proton pulses from the electron multiplier. The window 
the pulse height selector was arranged cover the range pulse amplitudes 
corresponding electrons the momentum selected the spectrometer 
and was adjusted for each setting the spectrometer. This pulse height 
selection served the double purpose greatly reducing the unwanted back- 
ground counting rate due gamma rays and further reducing the possibility 
detecting scattered electrons from neutron decay events that did not have 
the correct momentum corresponding the spectrometer setting. The setting 
the window was made two methods: (i) using electrons transmitted 
through the spectrometer from source mounted the sensitive volume 
the center the central electrode assembly and (ii) using electrons emitted 
from the surface the central electrode assembly stray radiation from the 
reactor beam. Figure illustrates the differential spectrum the pulses 
obtained using radioactive source and Fig. 3(6) the differential spectrum 
using the reactor beam. the latter case the continuous spectrum not associ- 
ated with the peak corresponding the electrons correct momentum was 
subtracted extrapolation under the peak. both Figs. and the 
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(a) 


300 400 300 400 


Fic. Differential pulse height curves for the electron counter. shows the spectrum 
from radioactive source and the spectrum obtained using electrons scattered from the 
surface the central electrode assembly with the inner gate open. 
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lower and upper settings the window are shown and might 
thought that settings obtained the first method would more reliable 
than those obtained the second; however, the second method has the 
advantage that carried out under the identical conditions operation 
used the actual coincidence runs and, after became apparent that the 
settings achieved the two methods were very closely the same, more and 
more reliance was placed the second method. 

The obvious advantages the use pulse height selection the output 
the scintillation counter are obtained the expense certain amount 
uncertainty how constant was the fraction the pulses within the 
window all the pulses due correct momentum electrons. Using curves 
the type shown Fig. the lower setting was chosen the basis 
constant fraction 92% the total area under the curve excluding the 
noise pulses. However, more reliance was placed the second method 
setting the window there was small uncertainty the actual fraction 
achieved. This uncertainty was considerably less than the statistical error 
the final counting rates and has been allowed for the estimate the 
final error. 

Before the experiment was set the reactor face, consideration was 
given system degaussing eliminate the effect the earth’s magnetic 
field. was discovered, however, that the earth’s field was less than 0.1 gauss 
over the spectrometer when located the reactor face. This low value was 
due the steel girder construction the building, which reduced the field 
generally throughout the building subnormal values, and the presence 
large steel platform located about feet above the apparatus and 
connected the thick steel face the reactor. Over the range energies 
involved this experiment this low residual field caused negligible distortion 
observed spectrum, and further degaussing was not employed. 

Sources RdTh, and were used establish momentum 
calibration for the spectrometer and measure its resolution and effective 
solid angle. these measurements the sources were mounted the center 
the sensitive volume and area approximated the cross-sectional area 
the effective neutron source. They were held the same potential the 
central electrode assembly that the 7-kv decelerating potential difference 
was present the same manner during the actual runs neutron decay 
events. Details the performance the spectrometer set for series 
and series are listed Table II. Figure shows the electron distribution 
from the dominant parity change transition plotted the 
ratio the observed counts per unit momentum interval divided 
where the electron momentum, the neutrino momentum, the Coulomb 
correction factor, and the unique shape correction factor (Rose 
al. 1953). The small contribution from the ground state transition and the 
background have been subtracted and the mean the points normalized 
unity. Perfect operation the spectrometer would give for all points; 
can seen that departures from this are small and considerably less than 
the errors the rates obtained from neutron decay events which lay the 
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TABLE 
CHARACTERISTICS THE BETA SPECTROMETER 


Series Series 
Geometric angles 
Measured solid angle 2.5% 0.8% 
Calculated solid angle* 2.6% 1.0% 
Measured resolution 7.7% 8.5% 

(full width half maximum) 

Magnification 1.0 0.65 
Diameter counter (inches) 3.0 1.5 


*The calculated solid angle is less than would be expected from the geometric 
angles of the entrance baffle system owing to the presence of structural supports 
not shown in Fig. 2 but required to locate the central parts of the various baffles. 


Cc} 

0.9} 


200 250 300 350 400 450 500 
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Fic. Electron distribution from the dominant parity change transition 
plotted the ratio the counts per unit momentum interval 


region 10%. Below 200 kev the points rise because finite source 
thickness and the neighborhood the end point they rise because errors 
the correction for the poor resolution the instrument and the imperfections 
the calibration. These checks indicated that the spectrometer was adequately 
linear the time the checks were made; allow for possible variations between 
calibrations, allowance was made the final error cover non-linearities 
and calibration errors which might have occurred during the long periods 
which the instrument was use. 

The use two entrance apertures the spectrometer defined the region 
the neutron beam from which decay events could detected. This was 
investigated experimentally mounting small source device which 
would allow its location relative the intercept the spectrometer axis 
and the neutron beam axis varied. this way the response the 
spectrometer could investigated function the location the source. 
was found that the effective neutron source was limited length approxi- 
mately three centimeters, and other directions the size the neutron 
beam. The form the limitation length shown Figs. and 14; 
further discussion these source dimensions will found Sections and 
this paper. 
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The current for the spectrometer coils was electronically stabilized 
comparing against standard cells the potential developed the current 
passing through cooled constantan resistor. The comparison was made 
400-c.p.s. chopper whose output after amplification controlled the generator 
supplying the current the coils. The short term and long term stability 
achieved this technique was better than 0.1%. The pressure the vacuum 
tank was maintained less than mercury throughout all 
runs means oil diffusion pump equipped with water-cooled baffle 
followed Freon-cooled baffle held —40°C. liquid-nitrogen trap 
was also employed and maintained full liquid nitrogen for the duration 
the experiment. was found necessary defrost and clean the Freon baffle 
and the liquid-nitrogen trap about once month. Apart from this trouble 
was experienced maintaining the required vacuum over the extended 
duration the experiment. 


THE PROTON COLLECTION 


The proton collection system was designed accept all the protons from 
the source volume that recoiled directions enabling them escape from 
the field-free region established the center electrode assembly. The escape 
direction was defined gridded circular aperture the axis the spectro- 
meter but the side the assembly furthest away from the spectrometer 
that the escaping protons had angles near 180° from the beta-spectrometer 
axis. The system was designed accept all protons from zero recoil energy 
the maximum possible recoil energy (755 electron volts) and was designed 
keep the proton counter far from the beam possible. With the geometric 
arrangements used the protons this experiment ranged recoil energy 
from 100 700 with most them the range 400 650 ev. The 
system can seen Fig. and will now described. 

Those protons which are produced the source volume with direction 
recoil which enables them emerge through the gridded are 
accelerated through across the gap between the central electrode and 
the gridded end the first circular tube held kv. They are then focused 
the symmetrical two-tube lens, formed the gap between the two circular 
tubes, through the entrance aperture the proton counter. 

The trajectories the protons through the first gap can accurately 
calculated since essence parallel plate condenser. With voltage 
difference across it, all protons the maximum possible recoil 
energy from any point the source volume can reach the second gap without 
striking the walls the first tube. higher voltage across this gap would 
course have had advantages far the further focusing the protons 
concerned but would have led undesirable perturbation the field 
inside the central electrode, and might have caused the acceptance protons 
the wrong recoil direction. 

The trajectories through the second gap and the proton counter were 
calculated using the data for symmetrical two-tube lenses given Klemperer 
(1953). The calculations were made for various proton recoil energies and 
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directions and allowance was made for spherical aberration. Spherical aberra- 
tion considerable help this application since reduces the focal length 
the lens for trajectories which depart considerably from the central axis. 
This tends compensate for the increase focal length which would otherwise 
occur owing the higher proton recoil energy causing the departure the 
trajectories from the central axis. The calculations indicated that with 
across the lens, all protons which reach the second gap and whose initial 
recoil energy lies between zero and the maximum recoil energy 755 volts, 
should pass through the entrance aperture the proton counter. Thus with 
across the first gap and across the second gap, the electrostatic 
system does not limit the source volume and can accept all protons for which 
the beta spectrometer can accept the corresponding electron. The calculations 
also indicated that because the large entrance aperture the proton 
counter the system should not all critical the voltage across the 
second gap. 

order test this system, coincidence runs were made with the voltage 
across the first gap fixed and with various voltages across the second 
gap. The runs were made with the beta spectrometer set 400 kev. this 
electron energy the proton recoils should have the widest spread energy 
and should thus provide the best test the system. The results these 
runs are shown Fig. from which can seen that the coincidence rate 
was essentially independent the voltage across the second gap between 


T 

W 4 


ACROSS TWO TUBE LENS 


Fic, Coincidence rate betwcen the proton and electron counter plotted against the voltage 
across the electrostatic focusing lens the proton collection system. 


limits 4kv and Below the rate dropped, presumably because 
the higher energy proton recoils were beginning miss the entrance aperture 
the counter because insufficient bending the second gap. Above 
trouble with breakdown occurred owing the high radiation field. The 
constancy between and gives evidence that the system 
forming predicted and was accepting all the protons for which the beta 
spectrometer could accept the corresponding electron. 

The proton detector was 14-stage electron 
plier using electrodes fabricated from commercial annealed beryllium copper. 


ROBSON: BETA DECAY 1461 


The arrangement the electrodes shown Fig. and was arrived 
trial and error using the steel ball and rubber membrane method obtaining 
the trajectories the electrons between the stages. The arrangement 


Fic. The electrcn multiplier used counter. 


similar that described Snell and Miller (1948) and the over-all behavior 
the multiplier was similar the behavior described these workers. 
The multiplier was activated induction heating the entire assembly 
about 510° pressure 200 microns normal untreated tank oxygen. 
The resulting stage gain the multiplier was about 3.3 per stage when 
operated 400 volts per stage and was very stable with time except during 
the first few hours after the voltage was applied. could made relatively 
gamma insensitive keeping the potential the first electrode ground 
potential that the secondary electrons from the iron shield were not attracted 
onto the first surface. The gamma sensitivity could, however, greatly 
increased raising the potential the first electrode about 200 volts 
above the grounded iron shield that the secondaries from the shield were 
attracted onto this first electrode. Magnetic shielding the electron multiplier 
was provided the form mild steel housing with additional soft iron 
shield near the first electrode. This shielding was sufficient make the gain 
the multiplier independent the current through the beta spectrometer 
coils. 
2(e). ELECTRONICS 

The electronics were conventional vacuum tube design and are shown 
block form Figs. and For all the runs series and for the first 
set runs series II, hereafter referred series Ila, the arrangement 
shown Fig. was used. With this arrangement the coincidences between 
the protons and electrons were obtained with the aid time sorter 
that their distribution time could examined. For this arrangement the 
proton side consisted wide-band linear amplifier, fast discriminator, and 
pulse shaper producing pulse 0.1-microseconds duration which was fed 
all the coincidence units. The electron side consisted the single channel 
pulse height analyzer and pulse shaper producing pulse 0.20 micro- 
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seconds which was fed through delay cables the coincidence units. 
minimum five and maximum eleven identical coincidence units were 
used with delay the electron pulse input 0.18 microseconds between 
each unit. For most the runs six units were used which formed time 
sorter covering the range 0.50 microseconds about 1.70 microseconds. These 
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Fic. Block diagram the electronics used record the coincidences series and Ia. 
Block diagram the electronics used record the coincidences series IIb. 


delays are such that the pulse from the proton counter must later than 
the pulse from the electron counter amount the above-mentioned 
range order record coincidence. Thus the first unit the time sorter 
recorded events which the proton pulse was between 0.50 and 0.80 micro- 
seconds later than the corresponding electron pulse. parallel with this 
time sorter single master coincidence unit recorded all events where the 
proton pulse was between 0.50 microseconds and 1.70 microseconds later 
than the electron pulse. 

further set runs was taken with different arrangement the electronics 
which the time sorter was not used but five separate discriminators and 
coincidence units recorded the coincidences different proton discriminator 
biases. For this set runs, referred hereafter series the geometrical 
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arrangement the apparatus was identical with that used series Ila. The 
electronic arrangement shown Fig. which, with the discussion 
mind, self-explanatory. 

Considerable effort was devoted making the circuits stable and 
reliable possible; particular the high voltages supplied for the counters 
were obtained from chopper type stabilizing units using standard cells 
reference potentials. 


2(f). NEUTRON BEAM MONITORING 


Three beam monitors were used, all recording independently. These were: 
(a) plutonium-fission ionization chamber mounted the beam catcher, 
(6) halogen-filled Geiger counter wrapped cadmium with high stability 
feedback d-c. amplifier monitoring the current from its collector wire, and 
(c) halogen-filled Geiger counter with standard pulse counting circuits. The 
halogen-filled counters were located close the beam the shielding. Monitors 
(a) and were used the actual analysis the runs with (c) being used 
remote indicating device and reserve monitor. The measurement the 
d-c. current from monitor (6) proved remarkably reliable method 
monitoring the beam intensity. 


EXPERIMENTAL PROCEDURE 
3(a). RANGE ELECTRON MOMENTA INVESTIGATED 

Before discussing the experimental procedure discussion will given 
the range electron momenta over which data can usefully taken. study 
Fig. will show that the electron—proton angle fixed and the electron 
momentum less than the neutrino momentum there only one proton 
recoil momentum which will satisfy the conservation linear momentum, 
whereas exceeds there are two possible proton recoil momenta which 
will conserve momentum. These two will referred the high momentum 
recoils and the low momentum recoils. For the neutron decay for 
1.07, equivalent electron kinetic energy 237 kev. For given 


(a) 


(b) 


Fic. The conservation linear momentum between the electron the neutrino and 
the recoil with fixed angle 160 degrees between and The top figure indicates the 
situation where less than resulting only one possible value the lower figure shows 
the two possible values which can cccur larger than 
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electron momentum possible calculate the corresponding ranges 
proton recoil momenta for the geometric arrangements used. These can then 
combined with the limits the distance which the protons travel the 
field-free region obtain the spread the transit times the protons the 
field-free region. The proton transit times the rest the proton system 
can then added obtain the total transit time limits for the proton from 
the decay events their arrival the electron multiplier. Representative 
values these total transit times are listed Table 


TABLE III 


TRANSIT TIMES IN MICROSECONDS OF THE PROTONS CORRESPONDING TO VARIOUS ELECTRON 
MOMENTA, Pp. HIGH AND LOW REFER TO THE HIGH MOMENTUM RECOIL GROUP AND THE LOW 
MOMENTUM RECOIL GROUP 


Series Series 
High Low High Low 


can seen from Table III that the spread proton transit times 
less than one microsecond for below 1.07 and above about 1.42. the 
region between 1.07 and 1.4, however, the spread for the low momenta recoil 
group quite large. the actual experiment was impractical record 
data where the spread transit times exceeded one microsecond since the 
required long resolving time the coincidence unit would have allowed too 
many random coincidences recorded. However, the range 1.07 
1.23 was calculated that the fraction the low momenta group was 
less than three per cent the high momenta group; consequently some data 
were taken this range and corrected for the loss the low momenta group. 
This correction was considerably less than the statistical error the data. 

The situation can summarized follows: 

(i) Below 1.0 the number events too few warrant the time 
necessary accumulate statistically useful data. data were taken 
this range electron momenta. 

(ii) From 1.0 1.07; only high momenta recoils occur and the result- 
ing spread transit times less than 0.5 microseconds. Data were 
taken this range. 

(iii) From 1.07 1.23; the number low momenta recoils less than 
three per cent the number high momenta recoils. Data were taken 
this range and corrected for the loss low momenta recoils. 


wil 
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(iv) From 1.23 1.45; the spread transit times the low momenta 
recoils large and their numbers are relatively large. Data taken 
this range were not used the analysis the results. 

(v) Above 1.45; the spread transit times the low momenta group 
less than 0.5 microseconds and overlaps that the high momenta 
group. Data were taken this range. 


TAKING THE DATA 


Data were taken counting the coincidences between the scintillation 
counter and the electron multiplier various settings the current through 
the beta spectrometer coils. each current setting data were taken five 
different bias settings the discriminator the electron multiplier channel. 
series and this involved taking five separate runs whereas series 
the data for the five bias conditions were obtained simultaneously. For each 
run series and both the time sorter data and the over-all coincidence 
data were recorded whereas series only the over-all coincidence data 
were available. Random coincidences were calculated from the individual 
counter rates and the resolving times the coincidence circuits, and were 
subtracted from the total coincidences give the genuine coincidences 
resulting from the neutron decay events. The resolving times were measured 
before and after each run measuring random coincidences with the aid 
radioactive sources. The electron multiplier was made more sensitive 
gamma rays for these random coincidence measurements the manner 
described Section 2(d). 

typical time sorter curve shown Fig. 1.46 run series 
Random coincidences have been subtracted and amounted about 100 
coincidences for each time sorter channel. this figure the proton transit 
time shown the center the individual channels that each point 
represents coincidences between about —0.14 microseconds and +0.14 
microseconds from the plotted abscissa. 
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Fic. 10. Time sorter results for series run taken 1.46. 
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prove that the coincidences were due neutron decay events, runs 
were made with (a) boron shutter located Fig. which eliminated 
the thermal neutrons from the beam, and with voltage the center 
electrode assembly that protons could not reach the electron multiplier. 
These checks, conjunction with the fact that the coincidences occurred 
the correct delay for neutron decay events, indicated that the coincidences 
were fact due neutron decay events. 

The total number coincidences could obtained for individual run 
series series Ila two methods: (a) dividing the area under the 
time sorter curve genuine coincidences the average resolving time 
the time sorter channels and using the number genuine coincidences 
recorded the over-all coincidence unit. These two methods agreed well 
the total number coincidences. the case series IIb only the latter 
method was available. From the individual runs set taken fixed 
spectrometer current integral bias curve was plotted against the electron 
multiplier discriminator bias. the basis large number such bias 
curves appeared that over the range bias settings used, the integral bias 
curves were well fitted exponential. best estimate for the number 
coincidences particular spectrometer setting was therefore obtained 
least square fitting the data the individual runs exponential bias 
curve and using the intercept zero bias. The purpose using this method 
taking electron multiplier bias curves each setting was overcome 
any possible long term drifts the gain the electron multiplier; during 
the experiment, however, the slopes the exponentials remained constant 
within statistical errors indicating the high stability the electron multiplier 
when operated continuously good vacuum. Figure shows one set 
runs plotted integral bias curve together with the corresponding least 
squares exponential. The normalizing the individual runs was made 
the basis the neutron monitor readings with arbitrary unit megawatts 
convenient measure the beam strength. The megawatt unit was 
established certain date correspond reactor thermal megawatt 


COINCIDENCES PER MW-HR 


ELECTRON MULTIPLIER BIASIN MILLIVOLTS 


11. Coincidence results for set series runs shown integral electron multiplier 
bias curve. The abscissa refer the pulse size the last electrode the electron multiplier. 
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but owing load changes the reactor our megawatt unit differed 
occasion from the reactor thermal megawatt per cent. general 
each run lasted for about 800 megawatt hours. 


THEORETICAL SPECTRUM SHAPES 


order obtain the theoretical spectrum shape was necessary 
estimate the efficiency g(¢) with which the apparatus could detect events 
which have angle between the direction the proton and the electron. 
The significance was discussed Section was estimated 
performing Monte Carlo type calculations with the aid the FERUT com- 
puter the University Toronto. The author deeply indebted Dr. 
Kennedy the Theoretical Physics Branch Atomic Energy 
Canada Limited for arranging and co-ordinating these Monte Carlo calcula- 
tions and the subsequent numerical integrations which led the theoretical 
spectrum shapes. 

essence the Monte Carlo calculations were performed selecting 
random point the effective source volume and then selecting proton 
recoil direction random. the resulting trajectory the proton was such 
that would detected the proton system this event was investigated 
further and corresponding electron direction was chosen fixed angle 
from the proton direction but with the plane the angle again randomly 
chosen. The trajectory this electron was then investigated see whether 
would enter the second entrance aperture the beta spectrometer and 
detected. could detected this event was recorded event, with 
certain value which could detected the system. repeating this 
sequence large number times and for different chosen values the 
electron proton angle was possible obtain the shape the variation 
with 

The input data required for the Monte Carlo calculations was the effective 
source volume. discussed Section 2(a) the effective source volume was 
defined the neutron intensity across the beam diameter plane per- 
pendicular the neutron beam and the acceptance the beta spectro- 
meter the direction the beam. The neutron density across the beam 
was measured the induced activity copper wires and shown for 
horizontal diameter Fig. 12. The vertical diameter distribution was almost 
identical. Also shown this figure are two alternative weighting functions 
used weight the random selection decay event location different Monte 
Carlo runs. The measurement the length distribution has been discussed 
Section 2(c); the results for series are shown Fig. together with four 
alternative weighting functions used the Monte Carlo runs. The corre- 
sponding results for series are shown Fig. 14. Several sets Monte 
Carlo runs were made corresponding different choices the weighting 
functions for the source dimensions. These are listed Table IV. 

The results the Monte Carlo calculation for set are shown for 
both series and series Fig. 15. These curves therefore indicate the 
efficiency the apparatus for detecting events where the angle between the 


CANADIAN JOURNAL PHYSICS, VOL. 36, 1958 


TABLE 


NOTATIONS INDICATE THE DIFFERENT MONTE CARLO 
RUNS CORRESPONDING DIFFERENT ASSUMPTIONS 
THE SOURCE DIMENSIONS 


Monte Carlo Radial Length 
set distribution distribution 
III 


1.0 
o\- 
0.5 


CENTIMETRES 


Fic. 12. The points represent the measured neutron density across the neutron beam 
the center the sensitive volume. The lines represent alternative approximations used 
the Monte Carlo runs. For fuller explanation see text. 


13 SERIES I 14 SERIES I 
05 Bos 05) 
fe) 2 3 -2 fe) ' 2 3 
2 IN CENTIMETRES 2 IN CENTIMETRES 


Fic. 13. The points represent the measured longitudinal distribution along the sensitive 
volume for series The lines represent alternative approximations used the Monte Carlo 
runs. For fuller explanation see text. 

Fic. 14. The points represent the measured longitudinal distribution along the sensitive 
volume for series II. The lines represent alternative approximations used the Monte Carlo 
runs. For fuller explanation see text. 


electron and proton with the source dimension assumptions listed 
Table IV. The results for the other Monte Carlo sets were quite similar for 
less than 170° but differed slightly the region from 170° 180°. 
From the Monte Carlo results for the corresponding functions and 
were computed desk calculator with the aid IBM 602A; these 
computations were made the computing staff the Theoretical Physics 
Branch Atomic Energy Canada Limited under the direction Dr. 
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Kennedy. Figure shows the functions and plotted against electron 
momentum for series Monte Carlo set together with plot the Fermi 
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ELECTRON MOMENTUM IN UNITS OF M,C 
Fic. 15. The efficiency the apparatus for detecting event occurring anywhere 
the source volume with angle between the direction the electron and the recoil 


proton. 
Fic. 16. The functions and obtained from the Monte Carlo runs for series under 


the assumptions the source volume listed set II. Also shown the Fermi distribution 
for the neutron decay spectrum. 


distribution for the neutron beta spectrum. The final theoretical 
spectrum shapes were obtained forming the products 


and correcting the resultant functions and for the resolution the beta 
spectrometer obtain the number electrons per unit momentum interval 


n(p) 
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where and are the resolution corrected functions and 
and the electron neutrino angular correlation coefficient. The resolution 
corrections were general less than two per cent and since the correction 
was being made theoretical curve for perfect resolution rather than the 
usual case where they are made experimental curve imperfect resolution, 
they could made with considerable confidence. The slowly varying Coulomb 
correction factor for beta-minus decay into nucleus charge was 
not applied the theoretical spectrum shapes but divisor 
the observed number events per unit momentum interval. The effect 
the Coulomb correction factor was almost negligible. 


EXPERIMENTAL RESULTS 


The experimental data accumulated the methods described Section 
were used give values for the observed number coincidences per 
unit momentum interval various electron momenta. These values were 
then compared the theoretical spectrum shape yield the 
least squares best value for and its statistical error. Least squares fits were 
made all the sets Monte Carlo runs listed Table with the exception 
that Monte Carlo set was not used for series data. Set was the first 
set tried the Ferut Computer and the particular values for 
which was calculated were not suitable for the geometric conditions 
series though they were well suited the conditions series The results 
the least squares fitting are listed Table can seen that the 
variation the resulting value for the different Monte Carlo sets 
considerably less than the error owing the statistics the counting 
rates. Since the Monte Carlo sets cover wide range assumptions about 
the dimensions the effective source can concluded that the error 
introduced not exactly reproducing the dimensions the Monte Carlo 
calculations not serious. 


TABLE 


VALUES THE ELECTRON NEUTRINO ANGULAR CORRELATION COEFFICIENT 

OBTAINED LEAST SQUARES FITTING THE EXPERIMENTAL DATA 

SERIES IIa, AND IIb THE THEORETICAL SPECTRUM SHAPES PREDICTED 


Monte Carlo 


set Series Series Ila Series 
Statistical error 


for all sets +.11 +.17 +.14 


The experimental data for series and series are plotted Figs. 
and together with the least squares best fits for the pure interactions. 
Monte Carlo set was used these figures. Figure shows both series 
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and series IIa data the same plot show the relative shapes the spectra 
for the two geometric conditions used the two series. The shift the peak 
the spectrum series corresponds the smaller mean value the 
electron proton angle. 
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Fic. 17. The experimental data series shown with the least squares best fits for the 
pure interactions. 

Fic. 18. The experimental data series shown with the least squares best fits for the 
pure interactions. 


The adopted values for from these three series runs were +.08 
for series +.01 for series Ila, and +.07 for series The errors 
these values will now discussed. 
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Fic. 19. The spectrum shapes series and series Ila shown the same scale show 
the differences resulting from the geometric conditions the two series. 


(a) Statistical errors. These were computed during the least squares fitting 
the data the theoretical curves using the statistical errors the coinci- 
dence counts. Expressed standard deviations they amount +.11 for 
series +.17 for series Ila, and +.14 for series The standard deviation 
the weighted mean series and +.11. 

(6) Error due the imperfect representation the actual source dimen- 
sions the Monte Carlo calculations. These were estimated +.03 for 
series and +.04 for the series runs from the variation for the different 
assumptions used the different Monte Carlo sets and listed Table 

(c) Error due possible beta spectrometer calibration errors. Errors 
the beta-spectrometer calibration directly affect the resulting computed value 
since they imply that the experimental points were fitted the theoretical 
spectrum shapes consistently incorrect momenta. estimate this error 
was made repeating least squares fits for various displacements the 
calibration. was found that for calibration error the resultant error 
was 0.11 for series and 0.13 for series Ila and IIb. The spectrometer 
resolution was 7.7% for series and 8.5% for series and IIb expressed 
full width half maximum. Allowance was made for possible calibration 
error 1/10th the resolution width leading possible errors +.09 
series and +.11 series and IIb. This the most serious the known 
systematic errors. 

(d) Error due possible non-linearity the spectrometer between cali- 
bration points. This was estimated making least squares fits for possible 
non-linearity 1/10th the resolution width over the range interest. 
less serious error than consistent calibration error because does not affect 
all the experimental points the same degree. amounted +.03 for 
series and +.04 for series and 

(e) Error due the window settings the pulse height selection the 
beta-counter output. This error was minimized making several independent 
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runs each spectrometer setting. was estimated allowing that 
individual run error +5% could have occurred the counting rate 
owing errors the window setting. This would lead error the 
values about +.02 all these errors were independent. This error was 
somewhat arbitrarily doubled allow for possible dependence the error 
the spectrometer setting. 

(f) Error due the Monte Carlo statistics. About 1000 successful Monte 
Carlo events were used derive the theoretical spectrum shapes for each 
Monte Carlo set. With four independent Monte Carlo sets this introduces 
standard deviation the resulting value for each series about +.05. 

combining these errors was assumed that series and series 
were taken under identical conditions. They were therefore combined form 
over-all result for series +.05 with standard deviation +.11. This 
was then combined with the other errors listed above provide final result 
for series this combination was assumed that the errors listed 
above were independent and should therefore combined quadra- 
ture. Series with its different geometrical arrangement was considered 
independently from series The final results for series and were then 
combined yield best value for the electron neutrino angular correlation 
coefficient the beta decay the free neutron +.07+.12. Table 
tabulates the errors and their combination into the final result. 


TABLE 
SUMMARY KNOWN ERRORS 


Series 
and 


Type error Series combined) 
Counting statistics +.11 +.11 
Source vclume errors 
Spectrometer errors +.11 
Window errors +.04 
Monte Carlo statistics 

+.16 +.18 


Over-all error 


DISCUSSION 


The result quoted for the angular correlation coefficient has assumed that 
the Fierz interference terms between the scalar and vector interactions the 
Fermi part the decay and between the tensor and axial vector interactions 
the Gamow-Teller part the decay are zero. From the evidence the 
case the part the decay (Sherr and Miller 1954) appears 
that the Fierz term, less than .06; the case the Fermi part the 
limit has not been well established but from the deduction Michel based 
revision analysis originally made Gerhart and Sherr the 
values transitions (Michel 1957) appears that less than .30. 
From study the spectrum Daniel and Schmidt-Rohr (1958) have 
concluded that less than .11. large value these interference terms 
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would have effect the interpretation the experimental data since 
would change the shape the theoretical spectrum shapes deduced from 
the Monte Carlo results. With limit 0.3 for the Fermi, and .06 for the 
part the Fierz term however, the effect would amount 
change less than .04. has therefore been neglected the analysis 
this experiment. 

The ratio the coupling constant the Fermi coupling 
Kofoed-Hansen and Winther 1956). Using this value the expected values 
the electron neutrino angular correlation coefficients are listed Table VII 
for four assumptions the form the beta-decay interaction. can 
seen that the experimental result consistent with beta-decay interaction 
the form but not consistent with interaction the form 
VT. 


TABLE VII 


EXPECTED VALUE FOR VARIOUS COMBINA- 
TIONS INTERACTIONS 


Form Angular correlation 
interaction coefficient 
.02 
+.47+.01 
This experiment 
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APPENDIX 
The formulae for the terms used the expressions describing the momentum 


distribution the electrons coincidence with the protons are given below. 


p/q ratio electron momentum neutrino momentum. 


The two signs arise from the possibility that the electron momentum 
exceeds the neutrino momentum there are two different proton momenta 
which will satisfy momentum conservation with electron fixed angle 


ROBSON: BETA DECAY 1475 


from the proton (see Fig. 9). These two cases are referred the text 
the high momentum recoils and the low momentum recoils. 
For both recoils considered together 
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ASYMMETRY THE CONCENTRATION FRONT 
DURING MISCIBLE DISPLACEMENT POROUS MEDIA' 


ABSTRACT 


During many feasible experiments concerning miscible displacement porous 
media, has been noted that the concentration front slightly asymmetric. 
possible that this due asymmetry the boundary conditions which 
present most practicable displacement experiments. The present paper en- 
deavors investigate the influence asymmetric boundary conditions upon the 
shape the concentration front: The diffusivity equation basic the theory 
miscible displacement has been solved for the case injection fluid constant 
concentration one end long, linear porous medium. The solution has been 
effected electrical analogue computer. Curves showing the asymmetry are 
given. 


INTRODUCTION 


The statistical approach the study flow through porous media has 
shown itself especially satisfactory the discussion miscible displacement 
processes which are important studies groundwater flow and new 
methods oil recovery. appears that the experimental evidence reported 
the literature general agreement with the predictions which the con- 
sideration statistical models leads. 

The usual solutions the differential flow equations (cf. Scheidegger 1957, 
1958) have one very important characteristic feature: They are all symmetrical 
about the mean concentration the fluid under consideration. the case 
linear displacement experiment, this means that the shape the invading 
concentration front exactly symmetrical about the half-concentration point. 
Most published data seem reasonable agreement with this pattern, 
but there appear cases where the displacement front skewed. 

One possible cause for the development skewed displacement front may 
that the boundary conditions during injection test are not really those 
that are usually postulated order arrive the symmetric shape for the 
front: The latter would occur only the case linear porous medium that 
extends infinity the positive well the negative direction. 
practice, however, porous medium cannot made infinitely long. One will, 
general, inject fluid certain kind one end and observe its 
intrusion into linear porous medium. expected that the shape 
the front will asymmetrical such experiment, simply because the 
boundary conditions are asymmetrical. 

the purpose the present paper investigate the effect the 
asymmetry the boundary condition upon the shape the displacement 
front. Only the ordinary statistical theory without autocorrelation will 

considered. The basic diffusivity equation has been solved for the prevailing 
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boundary condition electrical analogue machine. will shown that, 
for most common porous media, the asymmetry the concentration front 
due the asymmetry the boundary condition quite negligible under 
ordinary conditions. The conditions under which the asymmetry does become 
important are established and curves which must referred under these 
circumstances are provided. 


FORMULATION THE PROBLEM 
The differential equation for miscible displacement (without autocorrelation) 
linear porous medium 


9 —; 
(2.1) 


where the concentration the invading fluid, factor dispersion, 
and isa length measured from the mean position the injection front; i.e. 


(2.2) 


where the pore velocity injection and the true (length) co-ordinate 
the point under consideration. 

proposed solve eq. (2.1) for the case where fluid constant 
concentration 


injected one end into the linear porous medium. the primed 
co-ordinate system, the boundary condition (2.3) becomes 


(2.4) —vt) const. 


This problem, fact, can solved analytically introducing higher 
transcendental functions (Aronofsky 1957). However, the availability 
electric analogue computer the Imperial Oil Research Laboratory makes 
profitable obtain the solution that means, especially since the analytical 
solution still requires fair amount computation obtain actually dis- 
placement curve. 


THE ELECTRICAL ANALOGUE 

The differential equation for miscible displacement porous medium 
shown (2.1) diffusivity equation. well known that the solution 
the diffusivity equation can simulated electrical network con- 
taining resistors and condensers. order establish the relationships between 
the required solution the diffusivity equation (2.1) and the quantities 
measured electrical network, necessary analyze the latter. 

The basic electrical network which supply the solution the problem 
shown Fig. Here denotes the voltage the mth junction point, 
the mth capacitor, and the nth resistor. For our applications, all the 
capacitors and resistors are assumed identical. 
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The continuity condition for the current the 1)th junction 
point yields 
(3.1) 
The voltage across the resistor given 
Similarly, the current into the mth capacitor given 
(3.3) ic, 


Eliminating the currents from the above system and setting C,, 
(since all capacitors and resistors are equal) yields 


(3.4) 
This can written, since (the steps are equal 
and, the steps many 
(3.7) 
— on’ . 


This the diffusivity equation. 

order compare the diffusivity equation the electrical analogue with 
that describing miscible flow through porous media, one has bring both into 
dimensionless form. 


(a) Electrical Analogue 
Introducing the dimensionless variables 


(3.8) 

(3.9) 

where some constant reference voltage, yields 


which dimensionless. 


(b) Porous Media 
Introducing the dimensionless variables 


(3.11) 

(3.12) 

(3.13) 


yields 


a 
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(3.14) 
which again dimensionless. 


(c) Boundary Conditions 

Finally, one has formulate the equivalent boundary conditions for the 
two cases. 

the original, one has 


which is, dimensionless variables 
—vut 
3.16 
This signifies for the analogue 


where the number steps per unit time (i.e. the time which the 
boundary condition recedes one step). This leads the following equivalent 


vto 
3.18 
(3.18) 
(d) Comparison 


One can now compare the dimensionless forms the two cases (porous 
media and electrical analogue) and set table equivalents (see Table 
where the time the original has been denoted 


TABLE 
Analogue Original 
t/RC T/to 
V/Vo 
aRC 


Suppose now that run has been made the analogue with the parameters 
(and Vo, course) given. will shown that this can made 
correspond any porous flow problem with the values given. 

One simply has (from the very last equivalence) 


(3.19) 

hence for the time ratio analogue, original) 


and for the ratio 


n 
(3.21) 
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RESULTS 


The network indicated Fig. was set Imperial’s analogue machine. 
total 150 (RC) units was available, which half was 
used represent the displaced fluid units), the other half the displacing 
fluid (—units). Accordingly, time all the capacitors were dis- 
charged, all the capacitors were kept constant voltage. After each time 
step 1/a, one the capacitors was disconnected from the battery. 


Vv, R, V2 Vn R, Ras 


Fic. Network the analogue machine. 


Two series runs were made. series was equal seconds, 1/a 
was equal seconds/step. series II, was equal 0.8 seconds, 1/a 
seconds/step. Since data could read simultaneously only four points 
the network, several runs had made each series obtain complete 
set results. The results obtained are shown graphically Figs. 

The various figures show the skewedness the displacement front very 
well. However, the electrical analogue compared with average porous 
media, turns out that the skewedness can significant only the very 
beginning the displacement experiment. Thus, let typical linear dis- 
placement experiment performed follows: The pore velocity injection 
assume 


4 


PERCENT OF ORIGINAL VOLTAGE 


Fic. Data from series runs for short times. 
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PERCENT OF ORIGINAL VOLTAGE 


Fic. from series runs for long times. 


RC -0.8 sec 5 sec/step 
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PERCENT ORIGINAL VOLTAGE 
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Fic. from series runs. 


typical factor dispersion (Aronofsky 1957) 


series was equal seconds, equal 5.0 seconds; hence, from 
(3.20) and (3.21) 


(4.3) t/T 6.25, 
(4.4) n/x 125 


This means that the porous medium represented 125 units 
the electrical analogue, and that 6.25 seconds the analogue correspond 
second the porous medium. From Fig. appears that after 
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steps, the skewedness barely visible. This corresponds 11.2 seconds 
the miscible displacement process during which the front will have pene- 
trated 1.1 into the porous medium. 

For common values the skewedness the front should therefore 
quite insignificant. However, becomes large phenomena the 
immediate neighborhood the inflow face are under scrutiny, the curves 
shown Figs. 2-4 must referred to. 
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DISCOVERY 


ABSTRACT 


The neutron-deficient isotope has been discovered proton bombard- 
ment natural cerium, chemical separation, and mass determination. Its half-life 
hours and the end point its positron spectrum Mev. 
gamma rays were observed. 


INTRODUCTION 


the course general study neutron-deficient isotopes praseody- 
mium formed under proton bombardment cerium oxide, the extracted 
portion target exposed Mev clearly showed mass spectro- 
graph. Thus separated, this isotope repeatedly showed easily measurable 
half-life 1.5 hours. Radiations appeared limited the continuous 
positron spectrum with end point 1.7 Mev. The reason that the existence 
isotope with this half-life was questioned Handley and Olson (1954) 
that does not decay into with the 255-kev y-ray previously thought 
the main means identifying (Danby, Foster, and Thompson 
1958). 


EXPERIMENTAL PROCEDURE 


Irradiations 

The cerium oxide target was prepared heating available ammonium 
hexanitrate cerate* less than 300° avoid the formation troublesome 
insoluble oxide. Targets were prepared packing about the cerium 
oxide produced into 1/8-in. diameter tube formed wrapping 0.001-in. 
dural foil round drill rod. The free edges this foil were clamped dural 
block which turn was attached the water-cooled internal probe. This 
cooling was necessary prevent the formation the insoluble oxide during 
bombardment. Irradiations from 1/2 hours’ duration were carried out 
with 20- 40-Mev protons. 


Chemical Separation 

Cerium differs from its immediate neighbors that may exhibit valence 
well +3, the normal rare earth valence. The chemical separation 
praseodymium from cerium therefore depends upon the oxidation cerium 
its valence state and precipitation its phosphate, leaving the praseo- 
dymium solution. The detailed procedure was follows: 

The cerium oxide was dissolved thorough wetting with concentrated 
nitric acid and heating after the addition hydrogen peroxide. After disso- 
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lution, the hydrogen peroxide was removed repeated evaporation with 
nitric acid. The cerium milliliters ice-cold concentrated nitric acid was 
oxidized with potassium bromate. After dilution with cold water the cerium 
was precipitated phosphate the addition phosphoric acid. The pre- 
cipitate was coagulated heating about 50° and removed filtration. 
After the addition ammonium hydroxide the praseodymium was carried 
down with the hydroxides residual cerium and small amounts lanthanum, 
praseodymium, and neodymium impurities the target material (about 
micrograms all). This was dissolved nitric acid and reprecipitated 
hydroxide. After washing with acetone and air-drying, the precipitate was 
dissolved one drop nitric acid and transferred the ion source filament 
under heat lamp. The complete series operations normally took about one 
and one-half hours. 


Mass Separations 

The mass spectrograph Nier sector type instrument which the ion 
beam deflected along circular path 6-in. radius and through angle 
90°. The dispersion along the collector about 1.5 per mass unit mass 
150. The magnetic field 2700 oersteds supplied pairs magnetron 
magnets placed convenient locations round the magnet. 

After the loaded filament was mounted the mass spectrograph, ion 
collector the form copper foil in. X0.002 in. was carefully rubbed with 
fine emery paper, remove the oxide, and inserted into the plate holder. The 
closed system was then pumped down some minutes 
and the ion accelerating potential set 750 volts. Under these conditions the 
isotopes were focused near the high mass end the collector. The filament 
current was then slowly raised until the ion current began register the 
output meter the d-c. amplifier. After adjustment the focusing potentials 
for maximum ion current the filament temperature was again slowly raised 
until the ion current was about amp. few minutes the current 
became steady and then the accelerating potential was reduced 700 volts 
and the filament temperature set give ion current about 
amp. Under these conditions, the masses were focused near the center the 
collector and the sample was depleted half hour. 

After the copper foil collector had been removed from the spectrograph, 
no-screen X-ray film was gently pressed against its front surface and three 
pinholes were made through both foil and film they might placed 
register later. 

After suitable exposure (upwards minutes) the film was developed 
reveal lines corresponding the positions active isotopes. Estimates 
line intensities Mev and Mev are represented Fig. When dry, 
the film was again placed register with the foil and the positions active 
isotopes marked pin points. The foil was then cut into narrow strips 
containing the individual isotopes which were examined separately. 

The decay studies were made with end-on tube with 
mg/cm? mica window set lead housing and connected scaling circuit. 


wo ~ © 
MEV MEV 


Fic. line intensities (relative) for praseodymium isotopes produced the bom- 
bardment cerium protons energy and Mev. 


RESULTS 

Mass assignments were made with reference 4.5-hour the mass 
which was established Stover (1951) through its well-known daughter 
(Hollander, Seaborg, and Perlman 1953). strong 4.5-hour activity 
alone was observed under 20-Mev bombardment. 

Following 40-Mev bombardment, the same activity was observed and the 
mass the immediate neighbour with 2.0-hour half-life was confirmed. 

addition, previously unknown activity with half-life 1.5+0.1 
hours appeared mass 137. 

The decay curves for all three mass separated isotopes are shown Fig. 
The positron end-point the decay was shown 1.7 Mev from 
aluminum absorption curve analyzed the method Bleuler and Zunti 
(1946). gamma rays were observed from this isotope. 
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Fic. Decay curves for mass-separated praseodymium mass 137, 138, and 139. 
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SUCCESSIVE DECAYS FROM AND DECAY 
ADJACENT PRASEODYMIUM 


ABSTRACT 


The decay entirely 9.0-hour ground state was indepen- 
dently discovered. allowed electron capture, decays predominantly 
low-lying state with only branching 440-kev excited state. The 
440-kev transition least predominantly character. 

34.4-hour isomeric state produced the reaction decays 
254.5-kev transition the 9-hour ground state. 0.6% branching was 
found from the isomeric state directly small group apparently high spin 
states 

Studies mass-separated isotopes confirmed the existence 1.5-hour and 
The K-capture positron branching ratios observed the decay 
and show that these isotopes decay allowed transi- 
tions their cerium daughters. Characteristic gamma rays were observed. 


INTRODUCTION 


Since the discovery (Dahlstrom, Foster, and Thompson 1958) was 
not confirmed others (Handley and Olson 1954), this work was initially 
undertaken years ago discover the reasons. The clarification rests not 
only the fact that decays entirely newly observed 9-hour ground 
state Ce’, but also the recognition that the ‘36-hour’ y-ray formerly 
used identify reality arises from 254.5-kev transition from 
34.4-hour the ground state. Meantime these isomeric and ground 
states cerium were independently found Brosi and Ketelle (1955, 1956). 

addition, however, the mass separations and its neighbors 
and have permitted study positron and gamma-ray spectra 
from individual isotopes and this study has been supplemented more 
precise measurements conversion electron spectra from chemically separated 
sources praseodymium and cerium. Moreover, the positron 
branching ratios were measured for all four mass-separated praseodymium 
isotopes. each case, experiment compared with theory indicates allowed 
transition the ground state the cerium daughter. 


STUDIES MASS-SEPARATED PRASEODYMIUM ISOTOPES 
The detailed target arrangements well the methods used chemical 
separations for both praseodymium and cerium studies are given Appendix 
This also includes the methods preparing sources for the and y-ray 
spectrometers and the surface deposits for the hot filament source the mass 
spectrograph. 


received August 12, 1958. 
Contribution from the Radiation Laboratory, Department Physics, McGill University, 


Montreal Quebec, Canada. 
Accelerator Development Department, Brookhaven National Laboratory, Upton. 


Can. J. Phys. Vol. 36 (1958) 
1487 


| 


CANADIAN JOURNAL PHYSICS. VOL. 36, 1958 


(A) Mass 

Mass separation was essential the measurement X-rays and positrons 
from individual decay systems, since these radiations are common all 
isotopes investigated. Fortunately, for this purpose, ion emission from 
praseodymium about 100 times efficient from cerium. Consequently, 
each separated sample was initially almost entirely praseodymium and this 
made the analysis the decay the parent its cerium daughter much 
simpler. Owing, however, the time required prepare sources and the low 
over-all efficiency the spectrograph, these studies are fully effective only for 
isotopes greater than 1-hour half-life. 

The purity the mass-separated samples was ensured after redesign the 
collector system detailed Appendix II, where some improvements the 
efficiency the mass spectrograph are also described. 

(B) DETECTION 

The y-rays from mass-separated isotopes were detected and recorded 
28-channel (Bell and Skarsgard 1956). This instrument was also 
used the study annihilation radiation and X-rays. The Nal crystal 
detector moderate size contact with 6292 DuMont photomultiplier 
mounted lead castle lined with iron. earlier calibrations with standard 
the efficiency this detector for y-rays known function 
energy for source distances down from the crystal. This calibration was 
used later y-ray studies chemically separated sources. 

With the relatively weak samples from the mass spectrograph, however, the 
simultaneous handling four sources required the maximum counting rates 
realized placing the samples contact with the crystal cover. Each 
separated active foil was mounted between two layers scotch tape and then 
stretched across lucite ring which fitted round the crystal cover. When 
annihilation radiation was counted, additional 1/4 in. thick aluminum disk 
was placed over the foil, thus sandwiching the source between disk and crystal 
cover. 

extension the earlier calibration the detector was therefore made 
estimate the absolute efficiencies counting X-rays and positrons when the 
sources are close the crystal. Stronger sources X-rays lan- 
thanum kev) and positron annihilation radiation were counted when 
placed (1) and (2) against the crystal cover. From (1) the emission 
rates were determined and from (2) the counting efficiencies were calculated 
for the contact position. detailed discussion efficiences given Appendix 


(C) OBSERVATIONS MASS-SEPARATED ISOTOPES 


preliminary survey with the mass spectrograph was made collecting 
the active isotopes plane foil which was then placed contact with 
X-ray film. For irradiations Mev natural cerium, short contact 
exposures made immediately after mass separation revealed intense lines 
masses 139, 138, 137, and 136. This indicates that lighter isotopes, present, 
are too short-lived for this method study. 
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Following 55-Mev bombardment, separated activities were next allowed 
decay for hours and then contact exposure was made for several days. 
This showed lines masses 139, 137, 135, and 134. The absence activity 
138 and 136 was due the decay and stable cerium isotopes. 


Pr'87 

The decay annihilation radiation from this isotope consistently gave the 
1.5+0.1 hour half-life. y-rays this period were observed. upper limit 
about 10% the annihilation intensity can set for undetected y-rays. 
Such upper limits quoted this paper are for the uncorrected full energy peak 

However, the decay the X-rays showed half-life 1.5 hours growing 
into 9.0+1 hours (Fig. Evidently the 1.5-hour transition the 
previously unknown ground state which decays electron capture with 
9-hour half-life. 

After the 1.5-hour annihilation radiation had decayed insignificance, 
very weak 440-kev y-ray became clear 

The branching ratios K-capture positron emission for this and the 
following prasecdymium isotopes are presented Table 


TABLE 


Observed 


NOTE: 

(a) Handley and Olson (1954). 

(6) Dahlstrom (unreported). 

Column 2: K-capture to positron emission branching ratios. 

Column 3: positron maximum energies. 

Column 4: nuclear energy differences (Mev) calculated from the positron end points. 

Column 5: nuclear energy differences calculated from the branching ratios on the assumption of “allowed 
shape" transitions (Feenberg and Trigg 1950). 

Column 6: logarithm of the positron partial half-life, 4. = (1+/)t. 

Column 7: logarithm of the positron comparative half-life. ‘f,’ is given as a function of Wo by Feenberg and 
Trigg (1950). 


Evidently 34.4-hour does not appear the decay and is, 
course, the reason Handley and Olson did not detect the latter. 


Pr 
There are strong y-rays from this isotope. Indeed, the only indication 
rays with the known half-life 4.5 hours was possible 1150-kev peak con- 
sistently appearing with about the annihilation intensity (Fig. 
After the 4.5-hour parent had largely decayed, the 166-kev y-rays from 
140-day appeared with low intensity. 


t 
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Pr88 

displayed large number y-rays addition the annihilation 
radiation and X-rays. With half-life hours, decayed stable 
daughter, Intense y-rays were observed 300+10, 800+10, 1040+10 
kev, and 1590+20 kev (Fig. and f). y-rays still higher energy exist, 
their photoelectric peaks are less than 10% that for 1590 kev. 


Pris 

This isotope relatively short life gave low yield the mass spectrograph. 
The weak sources nevertheless gave half-life hours several 
occasions (Fig. g). This confirms previous assignment Handley and 
Olson (1954). However, y-rays were observed, and upper limit 10% 
the annihilation intensity can put any that may present (Fig. 


Additional Data from Chemically Separated Praseodymium 

Cerium was bombarded with 50-Mev protons and combined praseodymium 
activities were separated chemically. The conversion electron spectrum was 
studied spectrometer 2.1% resolution and 1.4% trans- 
mission. 

strong positron continuum and composite Auger electron peak decayed 
complex manner. 

and L-conversion lines for single transition kev were very 
prominent (Fig. The half-life was 2.1+0.05 hours and the K/(L+ M+) 
ratio was 2.2+0.1. 

Weak K-conversion lines arising from 800+7 and 160+2 kev transition also 


= 18900 
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Qa 
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Fic.2. (a) Conversion spectrum combined praseodymium isotopes chemically separated 
from cerium target following bombardment 50-Mev protons. (b) Decay 440-kev gamma 
ray daughter the praseodymium chemically separated from cerium target subjected 
40-Mev proton bombardment. 
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were present with 2.1-hour half-life. The relatively accurate value for 160 kev 
made possible convenient comparison with the known 166-kev transi- 
tion K-conversions for both 800 and 160 kev are about 3-5% the 
corresponding conversions 300 kev. 

few weak peaks observed with long half-lives will come later dis- 
cussions. They are K-conversions due (1) kev (2) 280+5 
kev and (3) somewhat stronger conversion 440+5 kev 


STUDY FORMED PRASEODYMIUM DECAY 


(a) Gamma-Ray Spectrum 

cerium target was bombarded Mev and the praseodymium 
produced was immediately separated chemically. After waiting hours for 
most the decay, the daughter cerium was extracted. 
Fortunately, only and 140-day are thus obtained with appreciable 
activities. these, the latter emits only one y-ray (166 kev). Thus the spectra 
may studied with minimum contaminating radiations. 

The cerium daughter (after repeated chemical purification cycles remove 
all traces praseodymium) revealed y-ray 440 kev and very strong 
X-ray peak. The well-known 254-kev y-ray, however, was not observed the 
cerium daughter. present, was less than the 440-kev peak intensity. 
Moreover, additional y-rays annihilation radiation were detected from the 
cerium. upper limit 0.25% the 440-kev peak intensity estimated 
for possible higher energy gamma rays 

The half-life the 440-kev transition, measured two occasions, was 
9.0+0.3 hours (Fig. 6). The X-radiation also decayed with 9.0-hour 
half-life until reached about its original activity. The mass separations 
had shown that the 9.0-hour activity was mass 137 and the chemical 
separations showed that this activity was associated with the decay cerium. 
The strong X-rays therefore indicate that decays electron capture— 
mainly the ground state with only weak branching 440-kev 
level. 

This branching ratio will considered later. Here important note 
that the 9.0-hour X-radiation fades into very weak contribution predomi- 
nantly arising from 140-day This very small background correction 
permits accurate measurement the X-ray intensity from 

Incidentally, the original cerium target does emit the 34-hour, 254-kev from 
which formed directly the target through relatively weak (p, 
reaction. 


(b) Conversion Electron Spectrum 

The praseodymium separated from cerium target which had been bom- 
barded Mev was allowed decay into cerium which was then chemically 
separated provide source (mainly for spectrometer. 

Auger electrons provided the only obvious peak. However, careful manual 
scan was made over the momentum region corresponding 440-kev total 
energy. the appropriate position, K-conversion peak was unmistakably 
observed with very low intensity. 
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(c) Branching Ratio 440-kev Level 

This branching ratio was estimated about the following manner: 

While the disintegrations the 440-kev level are obtained from the peak 
intensity and efficiency counting, only the K-electron captures are directly 
given the observed X-rays. These actually receive negligible contri- 
bution from the conversion 440 kev, but they fail indicate other types 
electrons captured and this supplement must now estimated. 

Theoretical probabilities for electron capture the and shells have 
been calculated Brysk and Rose (1955). Fortunately the fraction 
the total capture varies slowly with the disintegration energy over range 
extending considerably above the A-capture threshold. Here the 
contribution (which responsible for most the contribution) 
nearly constant fraction the From the absence positrons, 
the nuclear energy difference not large and since not critical take 

Under the above conditions, the theoretical ratios are 0.13 and 
0.01. The contribution negligible. 

For the remaining shells, estimate the contribution was 
made, assuming inverse cube variation with the principal quantum numbers. 

this way the ratio total captures captures was estimated 
1.21, Allowing for counting efficiencies and increasing the X-rays 21% 
get the total captures, was found that the ratio the disintegrations 
the 440-kev level the total disintegrations was 0.03. 


(a) Gamma-Ray Spectrum 

lanthanum target bombarded Mev (i.e., below the threshold for 
production) yielded mainly and with only 
relatively small intensity from disintegrations. This fact simplified the 
study gamma rays from 

The gamma-ray spectrometer disclosed very intense X-rays, strong y-rays 
255 and 440 kev, and much weaker group round 810 kev. 

Measured over five half-lives, the 255-kev y-ray had half-life 
hours and the 800-kev y-ray had the same period. the other hand, the 
X-rays and the 440-kev y-ray showed composite decay. After initial 
half-life which was fitted hours, they also decayed with 34.4-hour period. 
activity which could attributed the decay was detected. 

annihilation radiation was observed. The upper limit unobserved 
high energy y-rays can set 0.5% the weak 800-kev transition. 
unsuccessful search was also made for low energy rays down 10-kev limit 
detection (see, however, Brosi and Ketelle (1956)). 

After equilibrium was reached with the 34.4-hour half-life, the intensity 
ratios, averaged over several runs, for the 

X-rays 


The accuracy for the relative intensities the 255- and 440-kev y-rays 
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should equal that the counting efficiency values, about 10%. For 
the weaker 810-kev y-ray, better than 20%. The accuracy for X-ray 
intensities discussed Appendix 
(b) Conversion Electron Spectra 

Sources produced 30-Mev bombardments lanthanum were 
also scanned the lens spectrometer. Extremely strong conversion electron 
peaks due the 255-kev transition were observed with 34.4+0.3 hour 
half-life (Fig. 3d). typical conversion spectrum for the 255-kev transition 
reproduced Fig. 3a, and shows that K/L 2.9, allowing 25% for the 
M-+ contribution. 

The comparatively weak electron peaks for the 440+5 kev transition are 
also shown Fig. the K/L+ ratio 5+0.5. 

Using very strong source and counting manually over several days, the 
conversion spectrum for the groups transitions round 810 kev was obtained 
(Fig. 3c). Overlapping lines makes the assignment ratios difficult. 
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Activities the formed lanthanum target under 30-Mev proton 
bombardment. (a) and Conversion spectra from 255- and 440-kev transitions respectively. 
(c) Conversion spectra from close group transitions near 800 kev. (d) Decay 255-kev 


K-conversion line. 
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For the strongest lines, due 812+7 kev, the K/L+ ratio appeared 
The ratios the K-conversion intensities for the transitions 


255 440 812 kev 


Comparing the intensity the A-conversion peak for 255 kev with the 
absolute y-ray intensity, the K-conversion coefficient was derived, taking into 
account the known transmission the lens spectrometer. This gave 
conversion coefficient 4.8 within 30%. For and 0.5, this result 
indicates that transition (Rose and Goertzel 1955). 

Using source which gave c.p.m. the 255-kev K-conversion peak, 
manual scan disclosed low energy conversion peaks (above kev) except 
those due Auger electrons. unsuccessful search was also made for 
possible transition from some the 812-kev groups the 440-kev level. The 
upper limit for such transition less than the 440-kev 
intensity. 


Results with 180° Spectrograph 

direct method finding the element which the 255-kev transition takes 
place make accurate measurements the conversion electron energies 
and compare their differences with the corresponding differences electron 
binding energies for adjacent elements. 

The spectrograph has resolution 0.2% and gave results accurate 
about +0.2 kev (see Fig. 4). 


TABLE 
RESULTS WITH 180° SPECTROGRAPH 


Line 
Energy (kev) 214.02 248.13 254.41 
Observed 34.11 34.44 34.90 
Difference 33.87 34.26 34.70 39.2 40.2 
Difference 32.66 33.43 37.8 38.7 


From this comparison, clear that the 254.5-kev conversions take place 

Supplementary evidence for the classification the 255-kev transition 
was obtained through visual estimates the intensity ratios for the lines 
shown the microphotometer trace (Fig. 4). The resulting estimates 
L-conversion coefficient ratios now supplement the measured the lens 
spectrometer and the corresponding ratios for M3, M4, and transitions 
respectively are added for comparison (Rose and Goertzel 1955) (see Table 
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TOTALO.6% 
Mev) 


137 

STABLE 


Fic. Tentative scheme for successive decays from and 
y 


TABLE III 
2.9 4.4 3.0 1.8 
Lyu/Liu 0.25 to 0.5 0.7 0.41 0.3 


Lower order multiplicities either character are rejected the long 
half-life. The higher electric multipoles are eliminated the small contri- 
bution. comparison the observed with the calculated ratios establishes 
rather firmly the multipolarity the 255-kev transition. 


(d) Coincidences 

The 440-kev y-rays were coincident with X-rays; expected, the 255 
were not. 

Analysis Results Decay 

The decay the K-conversion electron peaks for the group transitions 
centered about 810 kev showed the simple 34.4-hour half-life. Thus, while 
decays predominantly the 254.5-kev transition weak 
branching occurs closely spaced group excited levels For this 
group, the y-ray spectrum was unresolved, and the ratios were too con- 
fused permit any assignments multipolarity. 

However, the branching the group whole was calculated from the 
intensity its unresolved y-rays, compared the total intensity the 
254.5-kev transition. The latter was obtained from the 254.5 y-ray intensity 
and the total conversion coefficient. This gave branching 0.6% the 
group whole. 
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The ratio the K-conversion coefficients the 440-kev and 254.5-kev 
transitions was obtained from the ratios their respective y-ray and 
conversion intensities. This gives 0.012+0.002 for the 440-kev transition, 
assuming the theoretical value for 254.5-kev transition 5.68. 
The interpolated theoretical values for 440 kev are 0.010 for E2, and 
0.017 for M1. 

The K/L ratio, assuming 25% the L-conversion 
peak, gave K/L 6.3+0.6 for the 440-kev transition. The theoretical values 
are K/L 6.0 for E2, and K/L 8.1 for M1, indicating the transition 
least predominantly nature. 

partial list new calculations the conversion coefficients, taking into 
account the finite size the nucleus (Sliv, privately obtained) have the for 
the 254.5 kev slightly reduced from the values Rose and Goertzel (1955). 
value about 10% less indicated, giving for the 440-kev transition 
0.013+0.002. The theoretical values Sliv for the 440-kev transition are, 


0.012+0.002, 
0.022+0.007. 


The K/L ratios will not significantly changed. Thus the assignments made 
are unaltered the Sliv calculations. 

The electron capture branching intensity the 440-kev level was also 
calculated from the 254.5- and 440-kev y-ray intensities. The total 
disintegration rate was obtained, using a7. From this, the disintegration 
rate was calculated for the two states secular equilibrium, and compared 
the 440-kev y-ray intensity, giving branching. This calculation 
independent counting X-ray intensities. 

This branching was also calculated terms the ratio the 440-kev 
y-ray the total X-ray intensity, for and secular equilibrium. 
The disintegration was expressed terms cerium X-rays, using 
ax, ar, and the fluorescent yield. The disintegration was expressed 
terms lanthanum X-rays, using the fluorescent yield and the total 
capture ratio. 

Ignoring the small percentage difference the counting efficiencies for 
cerium and lanthanum X-rays negligible effect the answer, which 
itself small), all the X-rays can added indiscriminately. From the ratio 
the 440-kev intensity the disintegration rate expressed terms 
the X-ray intensity, branching was obtained. 

earlier section, the value was also calculated for the decay the 
9-hour ground state alone. Since both calculations from X-ray intensities give 
the same result, this indicates the self-consistency the interpretation 
the decay. The value obtained independent X-ray measurements was 2%. 
However, the discrepancy not large considering the many factors involved. 
The total counting efficiencies for the X-rays, error, are most likely 
too large. This would make the branching slightly less than 3%. the 
other hand, reducing for the 254.5-kev transition, indicated the Sliv 
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results, would raise the figure somewhat. The three calculations give 
essentially the same result, about 3%. 

analyzing the decay estimates the nuclear energy 
differences involved were arrived at. positrons were observed the decay 
and upper limit the intensity can conserva- 
tively placed them. This indicates the disintegration energy the 
ground state cannot more than about 1.5 Mev. Furthermore, lower 
limit about 440 kev can set, since there branching the 440-kev 
level. 

For K-capture, the ‘‘f,’’ factor the comparative half-life can found, 
knowing the nuclear energy difference Wo, and (Feenberg and Trigg 
1950). Calculating the partial half-life the ground state, and between 
the limits energy given above, the value 4.25 log 5.5. 

For the transition from the 440-kev level, applying the appropriate 
energy limits and partial half-life, the value 


Both these results indicate allowed transitions disintegration energy about 
the high end the estimated range. 

energy systematics (Way and Wood 1954) support this result, 
predicting 1.7 Mev for ground state ground state 
transition. Taking 1.5 Mev gives about Mev for the energy difference 
between and the 812-kev level. Using the partial half-life the (0.6%) 
group whole, this gives value 


log 7.8. 


Consider the levels 743, 812, and 901 kev. Since the energy the electron 
capture transition high all cases (about Mev for the 812 kev), the 
contribution considerably smaller than the part, and 
furthermore, doesn’t vary much over the range energy involved. The 
intensity the three transitions are the same order. Thus their partial 
half-lives vary certainly less than order magnitude, and the 
logarithms will vary less than unity. Thus log applies for all three, 
indicating that first forbidden, first forbidden unique, transitions are 
probable. 


DISCUSSION 


(1) The Mass 137 Chain 

The decay 1.5-hour through 9-hour and long lived 
stable produces intense characteristic y-rays. The intense 254.5-kev 
transition which had been erroneously considered from the ground 
state completely absent this decay. 

The information obtained agrees with that Brosi and Ketelle 
They also obtained 8.7 hour half-life, decaying electron capture 
with branching 445-kev level. They assigned the transition 
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from this level either M1, which essential agreement with 
the present results. 

The 180° conversion electron spectrograph definitely established that the 
254.5-kev transition cerium. Furthermore, since the 9.0-hour 
grows into secular equilibrium with the 34.4-hour period, and since the mass 
assignment the former was directly established, the mass assignment 
the isomer also definite. 

Brosi and Ketelle (1955) also produced a-particle bombardment 
They obtained 34.5-hour half-life for the isomeric state, and 
the basis its K-conversion coefficient and critical absorption experiment 

The low intensity group 34.4-hour transitions with the main component 
812 kev indicates electron capture directly from the isomeric state 
excited states These levels not cascade through the 440-kev level. 

The spin has been measured 5/2. The single particle Shell Model 
(Mayer 1950) predicts ground state for (La), and for 

The ground state predicted probably with the not 
lying much above it. plot the energy differences between the and 
the levels for odd odd isotopes function indicates small 
energy difference for with probably the ground state (Goldhaber 
and Hill 1952). 

For odd nuclei with the number the odd nucleon between and 82, all 
states have even parity, except when the odd nucleon the excited 
state. does not decay the level, the transition involves parity 
change and the log value makes the allowed assignment unambiguous, 
since second forbidden would inconsistent with the results. Thus the pre- 
dicted transition supported the experimental observations. 

ground state for would inconsistent with observations the 
decay unless lay just few kiloelectron volts below the state. 
However, the decay supports the ground state, will shown 
later. This similar which also has ground state (Ketelle, 
Thomas, and Brosi 1956). 

Fig. 4(b) decay scheme for proposed. 

The isomeric transition from consistent with predictions based 
neighboring odd odd isotopes (Goldhaber and Hill 1952). plot 
the energy difference between the and levels, for 79, indicates 
about 310 kev the expected value for 

evidence for any decay through was found, which 
consistent with the above assignments predicting third forbidden transition 
This shows, further, that could not have 
ground state, since the transition the isomeric state would then favored 
over that the ground state. 

The independent particle model suggests ground state for with 
the level lying just above, and with additional and excited states. 
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With allowed transition from the ground state taking place, 
ground state was indicated, unless the level lay only few kiloelectron 
volts below the This supports the ground state since 
would inconsistent. Stable with measured spin 3/2, has 
ground state. 

The interpretation having ground state would consistent 
with the observed decay including allowed electron capture the 
440-kev level followed transition. However, the decay the ground 
state would then allowed. get very long half-life would require 
extremely low energy for the transition. The systematics (Way 
and Wood 1954) indicate value the order 200 kev for this transition. 
alternative interpretation, ground state for would produce 
second forbidden transition, and thus account for the long half-life (see below). 

Brosi and Ketelle (1956) produced bombarding natural cerium with 
slow neutrons. They obtained barium X-rays the decay and estimated 
the half-life years. addition they found 10-kev transition 
prompt coincidence with X-rays and the 440-kev transition. 
The intensities indicated that practically all capture events through the 
transition. Their results completely supported the level the ground 
state, the 10-kev transition being from the level the ground state. 

The low-lying levels would appear single particle levels. The 
440-kev transition least predominantly character, which could 
indicate strong transition from the L-forbidden state. support this 
interpretation, the level would decay transition the 
ground state. The half-life this transition would long compared the 
cascade through the level, which consistent with observations. the 
other hand, for initial level, predominantly transition the 
level would expected, and appears less favored. any event, the evidence 
for predominantly transition. 

The values log indicate probable first forbidden, first forbidden 
unique transitions from the three high energy levels Thus 
their parity should opposite that the state, the same the 
7/2 ground state Furthermore, high spins >9/2 are favored since 
cascade transitions through the 440-kev level occur. These levels not 
appear easily interpretable. Purely from the viewpoint the spins and 
parities favored, the results would consistent with the presence rotational 
levels the states. High rotational energies are expected near 
closed shells. However, the evidence very meager for fruitful speculation. 
produced decay, these transitions are very weak for detailed 
study, and not certain which low level the transitions occur. 


(2) The Adjacent Isotopes 

The decay 4.5-hour produces y-rays, and its daughter activity 
exhibits the 166-kev y-ray and 140-day half-life indicating decay the 
ground state. The value predicts allowed transition. The ground state 
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140-day has been established state (Ketelle, Thomas, and 
Brosi 1956). 740-kev isomeric transition from level the 
ground state was observed these workers, with half-life seconds. 
the decay transition this level would energetically allowed, 
but 740-kev radiation observed, indicating that such transition 
highly forbidden. ground state for expected, the case 
and consistent with allowed decay transition the 
state would third forbidden. the case, the absence parity 
changes makes the allowed assignment unambiguous. The level 
ground state ruled out the absence decay the isomeric state. 
Thus decays allowed transition the ground state 

decays producing very intense y-rays 300+5, 800+10, 
and 1040+10 kev. Less intense y-rays 1295+15 and 1590+20 also appear, 
and weak transition 160 kev. The value indicates allowed decay. 
For allowed transition and tentative assumption that the decay 
takes place predominantly single level the observed 
positron branching ratio gives correct value for the nuclear energy 
difference, obtained from the positron end-point energy. This agreement 
supports the assumption, and indicates the decay mainly single excited 
level, closely spaced levels. 

the ground state has its odd proton the level and its odd 
neutron the level, the neighboring odd nuclei, this would give 
and even parity for the ground state, Nordheim’s Strong Rule. 
But this favors intense allowed transition the (0, ground state 
stable The presence intense high energy y-rays plus the 
positron branching ratio opposes this assumption. seems probable that 
the ground state has high spin value and decays allowed 
transition high spin excited state, states, The beta-decay 
systematics Way and Wood (1954) indicate value about 3.7 Mev 
for the ground state This would favor the transitions pre- 
dominantly the vicinity the 1.3-Mev level which would 
consistent with observations. 

The 300-kev transition seems definitely character, which would 
indicate isomeric state with half-life the order seconds. For the energy 
and involved, the ratios vary widely with multipolarity, making the 
assignment seem unambiguous. The new calculations the conversion 
coefficients Sliv increase the coefficient for the 300-kev 
transition only about 10% above the Rose value, and the K/L ratio 
not expected change significantly. Thus the assignment still stands. 

Investigations the decay (Turchinetz and Pringle 1956; 
Glover and Watt 1957) show 810-kev (2, level The ground 
state (5, has too low energy decay the 1040 higher 
energy levels that the results these workers throw light the 
high energy levels. However, this indicates that the intense 300-kev transition, 
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and probably the weak 160-kev transition, are not fed from the 800-kev (2, 
level. Furthermore, shows that these transitions were from low-lying 
excited states rather than from cascade transitions between high energy levels, 
the decay these levels was highly forbidden spin basis, since the 
observed decay the 800-kev level has log (Glover and Watt 1957). 

purely energy basis, the 300-kev transition could cascade from the 
1590- the 1295-kev level. The 1590- and 1295-levels could have cascade 
transitions through the 800-kev (2, level well directly the ground 
state; since the resulting 790- and 495-kev transitions could concealed 
this investigation the 800-kev level decay itself, and the 511-kev 
annihilation radiation. However, the intense 1040-kev level appears have 
spin and parity such that does not cascade the 800-kev level, since 
radiation 240 kev was found. further possibility would 300-kev 
(3, level fed cascade from higher energy levels. For example, the 
1295-kev transition could cascade from the 1590 level through the 
300-kev level. The assumption the 800-kev (2, level the first excited 
state would consistent with the systematics energies (2, first 
excited states even nuclei, which indicate that about this energy expected 
for (Scharff-Goldhaber 1953). 

summary, the evidence the excited states too sketchy for 
proper organization. However, the existence for this even nucleus 
isomeric state seems fairly definite. 

Pri*6 

The value for the decay 1.0+0.15-hour indicates allowed 
transition. the previous cases, only the transition allowed does the 
branching ratio give the correct value for the nuclear energy difference. 

has its odd proton and neutron the same levels the neigh- 
boring odd nuclei, Nordheim’s Strong Rule predicts ground state 
even parity. This would give allowed transitions the (0, ground 

The beta-decay energy predictions for even nuclei Way and Wood 
would support decay level round 1100 The systematics the 
energies Ist excited (2, levels even even nuclei suggest about 900 kev 
the case (Scharff-Goldhaber 1953). 

Unlike the spectra revealed y-rays conversion 
However, the mass separated activity was weak, and the upper limits 
the intensity any unobserved y-rays cannot rule out the possibility 
high energy y-rays low intensity. Similarly high energy transitions low 
intensity could produce conversion electron peaks too small observed 
these experiments. However, the presence low energy transitions can 
ruled out, except for very low intensities. 

Thus, apparent contrast the present evidence supports low 
spin ground state for which decays predominantly the ground state 
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APPENDIX 


TARGETS, CHEMICAL SEPARATION, AND SOURCE PREPARATION 

(a) Targets 

For the production praseodymium isotopes, targets natural cerium 
were used both oxide* and forms. About the oxide was 
packed into aluminum tube diameter long made from 
0.015 foil. This target was cooled through intermediate dural clamp. 
Nevertheless, the cyclotron beam had reduced 0.5 avoid the 
relatively insoluble oxide which forms high temperatures. Alternatively, 
bare piece flattened metal from storage (to reduce oxidation) was 
mounted directly one the dural clamps. Following broadside bombard- 
ment about the front 1/16 in. was cut off and chemically separated. 

For the production targets natural lanthanum oxide were used. 
This material was handled the same manner the cerium oxide, but the 
bombardments were carried out ya. 


The Chemical Separations 

The chemical separations were based the fact that cerium can oxidized 
valence state while lanthanum and praseodymium remain the 
normal valence. 


Cerium from Lanthanum Targets 

The lanthanum metal oxide was dissolved concentrated nitric acid and 
after cooling ice bath, potassium bromate was added oxidize the 
cerium its valence. After the addition zirconium carrier, saturated 
solution potassium iodate was used precipitate the cerium with the 
zirconium. After several purification cycles dissolution and reprecipitation 
remove the lanthanum, the cerium was reduced the valence state with 
hydrogen peroxide and the zirconium was removed iodate precipitate. 
The cerium was carried from the solution with ferric hydroxide. After washing 
and drying, this precipitate was dissolved few drops and the 
iron was extracted with ethyl ether leaving the radioactive cerium small 
volume acid solution. 


Praseodymium from Cerium Targets 

The early praseodymium separations were carried out using method 
which the valence+4 cerium was precipitated phosphate (Dahlstrom, 
Foster, and Thompson 1958). However, after the development the iodate 
separation for lanthanum targets, this was also adopted for praseodymium since 
the conditions seemed less critical for the iodate than for the phosphate pre- 
cipitation. 

The cerium oxide was dissolved being heated concentrated nitric acid 
and hydrogen peroxide. After dissolution the hydrogen peroxide was removed 


*Cerium oxide was formed heating ammonium hexanitrate cerate (G. Smith Chem. 
Co., Reference Purity) about 250° The lanthanum oxide was the Spectroscopic Grade 
Johnson and Matthey. Cerium and lanthanum samples metal, each containing less than 
parts any other rare earth, were obtained from Ames Laboratory through the kindness the 
Director, Dr. Spedding. 
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repeated evaporation with nitric acid. Since the cerium metal was insoluble 
nitric acid was dissolved minimum volume hydrochloric acid. 
Repeated evaporation with concentrated nitric acid removed the chlorine. 

The cerium ice-cold concentrated nitric acid was oxidized with potassium 
bromate and precipitated the iodate the addition saturated solution 
potassium iodate. Zirconium carrier was added the supernatant solution 
remove any cerium which was not precipitated the first time. The praseody- 
mium was removed from the iodate solution hydroxide, using iron 
carrier. The iron was removed from small hydrochloric acid solution con- 
taining the praseodymium activity, which was then ready for the preparation 
sources. 

However, for mass spectrograph use, the praseodymium must applied 
nitrate the filament, and speed was essential, even the cost some loss 
efficiency. When the iodate solution was rendered basic with ammonium 
hydroxide, the praseodymium was carried down the ammonium iodate 
crystals themselves, from which was separated judicious washing. Com- 
plete removal salts was ensured dissolving drop concentrated 
nitric acid and reprecipitating with ammonium hydroxide. 


SouRCE PREPARATION 
Mass Spectrograph 
The filament electrode was placed horizontally beneath heat lamp. When 
the filament was rather warm, drop acid solution containing the activity 
was carefully applied the central region micropipette. This dried 
quickly and the source was then installed the spectrograph. 


Gamma Ray Spectrometers 

similar drop active acid solution was applied the surface thin 
teflon sheet scotch tape. The drop was dried and then protecting layer 
scotch tape was placed over the source. 


Spectrometer 

These sources were commonly mounted film made from 
VYNS Resin obtained from the Chemistry Department. The films were very 
lightly gold-plated ensure electrical grounding. The spot size was defined 
drop insulin, later removed washing. The active drop applied 
this spot was dried with great care. occasions, when speed was essential, 
the active drop was dried thin piece niobium foil. 


IV. 180° 

30° trough 1/8 in. deep was cut across 1/2 in. wide piece teflon. Then 
fast rotating cutter made smooth sharp angled cut few thousands 
inch deeper the bottom this trough. 0.002 in. diameter tantalum wire, 
mounted the source holder, was placed this finer trough. The drop 
containing the activity pipetted into the center the trough. Owing 
the non-wetting the teflon, surface tension holds the drop the center. 
After drying with heat lamp about one-third the activity remains the 
wire surface and confined short region. 
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THE Mass SEPARATIONS 


The mass spectrograph was the same that described the previous 
paper except for alterations described below. 


Transmission 

The horizontal focusing plate system was replaced radial lens which 
provided some direct control the beam spread the vertical plane. Machined 
teflon spacers ensured exact and reproducible alignment lens components. 
produce line width the collector equal one-half the linear mass 
dispersion for praseodymium, viz. 0.020 in., required beam-defining slit 
0.015 in. wide. The transmission efficiency for the whole system was then 
about 15%. Owing, however, the low electric fields the hot filament 
surface and the consequent limitation the emission singly-charged ions 
praseodymium oxide, the over-all efficiency was 


Collection 

Initially, the separated isotopes were collected plane copper dural 
foil. means contact print X-ray film, the active isotopes were 
identified, cut out, and monitored. 

This procedure gave positive identifications mass numbers; but imperfect 
separations learned from the y-ray spectrometer. The main source such 
imperfections was not the resolution the instrument but either mechanical 
smearing deposits during manipulations the X-ray film errors 
cutting the foil guided only pin pricks. With little 
loss activity, these troubles were avoided system described the 
next paragraph. 

the central portion the collector foil, the individual ribbon collectors 
were precut along their sides with the known spacing the linear mass 
dispersion, 0.040 in.; but through their uncut ends, these remained 
integral part the collector foil. This foil was insulated from its grounded 
front supporting frame 0.030-in. teflon washer. Thus the collection area 
was connected, before, the amplifier circuit which indicated the ion 
current received. However, spot-welded the grounded frame immediately 
front was grid opaque vertical ribbons width 0.4d and spacing, 
Each member this grid was obviously centered front corresponding 
slit the collector foil. Thus control the accelerator voltage, the mass 
samples could brought through the 0.6d openings the grid the center 
precut ribbon collectors. For convenience, this voltage was selected that 
the stable (which alone gave detectable current) passed through the 
first opening. Continuous monitoring the accelerator voltage kept the 
collector yield maximum. Thus the mass number for each collector strip 
was predetermined. Thereafter the pure samples known masses were readily 
separated. 
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THE ABSOLUTE DISINTEGRATION RATES 


obtain the absolute disintegration rate for particular y-ray, the observed 
counting rate the full energy peak was divided efficiency factor. Using 
sources known strength, these efficiency factors for energies greater than 
100 kev had been calculated function distance from the crystal cover, 
(Bell and Skarsgard 1956). For energies below 100 kev, the applied efficiency 
factors were the computed total efficiencies McGowan (1954) also using 
in. diameter counter. low energies these total efficiencies are 
approximately equal the photoelectric efficiencies. However, the photo- 
electric effect produces both the full energy peak and the iodine escape peak, 
which kev lower energy. For X-rays about kev this escape peak 
too low energy detected and this unknown part the total intensity 
lost. 

obtain estimate the correction required for the escape peak, 
source lanthanum X-rays was sandwiched between two identical 
1.5 in.X1 in. Nal detector units placed end end. The combined 
outputs were 15% greater than the sum the two counter rates when 
separated, with the source contact with each turn. This indicated that 
15% the full energy peak intensity was escaping into the opposite crystal. 
This did not represent the entire escape peak, however, since some losses occur 
owing geometry and absorption the aluminum covers. 

similar experiment was performed with 57-kev tantalum X-rays where 
both the escape peak and the peak due 28-kev iodine X-rays from the 
opposite detector were observed. This showed that 60% the escaping 
iodine X-rays were detected the opposite counter. Assuming that the 
efficiencies were about the same both cases, 25% the lanthanum X-rays 
were not represented the full energy peak when the source was directly 
against the crystal cover. 

With single counter, the escape peak for tantalum X-rays was 22% the 
full energy peak with the source set directly the detector, and 16.5% 
Assuming the efficiency corrections were about the same both cases, 
this would indicate that about 19% the lanthanum X-rays were not repre- 
sented the full energy peak with the source distance for which 
the efficiency factors are tabulated for the detector). 

further complication arose from the fact that while cerium and praseo- 
dymium X-ray components were all above the iodine K-edge, all lanthanum 
components were not. The lanthanum Kaz component was below the iodine 
K-edge energy, and thus contributed nothing the escape peak. The 
ratios the intensities the are respectively about 2:1:1 
(Richtmyer and Kennard 1947). Thus only 75% the lanthanum X-rays 
contributed the escape peak and the efficiencies were lowered 
25% for cerium and praseodymium and 19% for lanthanum X-rays. The 
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limit accuracy the efficiency factors for y-rays energy greater than 
100 kev was about 10%. For the X-rays, considering both the escape peak 
corrections and the fact that the calculated total counting efficiencies used 
might somewhat large, the limit error might high 

addition the above factors, the observed X-ray intensity must 
divided the fluorescent yield give the total intensity A-shell events. 
While this factor not due the counting process itself, included 
this point for completeness. For the X-rays cerium and lanthanum the 
fluorescent yield about 0.88 (Bergstrom 1955). 


Studentships from the National Research Council were awarded one 
(G.T.D.) for the period the investigation. 
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THE STABILITY COUETTE FLOW AXIAL 
MAGNETIC FIELD! 


NIBLETT 


ABSTRACT 


theory the stability viscous flow electrically 
conducting fluid between coaxial rotating cylinders with perfectly conducting 
walls extended include the case non-conducting walls, and found 
that their effect reduce the critical Taylor numbers and increase the wave- 
length the instability patterns considerable amounts. experiment 
designed measure the values magnetic field and rotation speed the onset 
instability mercury between perspex cylinders described. The radioactive 
isotopes and were used trace the flow. The results support the 
theoretical prediction that the boundary conditions can have large effect 
the motion. 


INTRODUCTION 


The problem the stability viscous flow between coaxial rotating 
cylinders (Couette flow) has been investigated Taylor (1923), Lewis (1928), 
Meksyn (1946, c), and others. Taylor’s work this connection was the 
first successful correlation between theory and experiment problems 
hydrodynamic stability. Within the framework certain approximations, 
was able define the conditions under which instability would occur, 
and also predict the motion which would develop after its onset. Taylor’s 
criterion for the onset instability depends the value the dimensionless 
parameter 


where 


2 2 
0 


the kinematic viscosity, the radial separation the cylinders, the 
speed the inner cylinder, the ratio the speeds the two cylinders 
and the radii the inner and outer cylinders respectively, and 
called the Taylor number, and its critical value when the 
cylinders rotate the same direction 1708. 

Chandrasekhar (1953) has studied the effect axial magnetic field 
the stability Couette flow when the fluid considered electrical con- 
ductor. The critical Taylor numbers were found increase quite rapidly 
the magnetic field increases, indicating that the effect the field inhibit 
the onset instability. The vortex compartments which appear the transi- 
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tion were found elongated the direction the field. The inhibiting 
effect the field was expressed terms the Hartmann number 

(2) 

where the permeability the fluid, the conductivity, the density, and 
the intensity the magnetic field. 

The practical application Chandrasekhar’s theory limited his choice 
boundary conditions which depend the cylinder walls being perfectly 
conducting. this paper the characteristic value problem solved for the 
case non-conducting walls, and some experiments are described which 
mercury was used between cylinders made from perspex. 


THEORY 


(a) The Equations Motion and Marginal Stability 
The governing equations are the basic equations magnetohydrodynamics 
for homogeneous incompressible fluid. These are: 


where the fluid velocity, the pressure, and the 
Electromagnetic quantities have been expressed e.m.u. 

The problem hand differs from that solved Chandrasekhar (1953) 
only respect the boundary conditions, and the formulation will not 
repeated detail. The perturbations velocity and magnetic field intensity 
superimposed the steady-state solutions the foregoing equations are 
assumed the form 


where are functions only. Here the variables are expressed 
terms the cylindrical polar co-ordinates (r, and disturbances are 
assumed symmetrical about the z-axis and therefore independent 
For the case when the cylinders rotate the same direction and their separa- 
tion, much less than the mean their radii, Ro, the equation marginal 
stability derived Chandrasekhar 


(8) 


This equation valid provided Radial distances from Ro, and 
the wave number have been expressed dimensionless form setting 


jj 


— 
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(b) The Boundary Conditions 

Eight boundary conditions must prescribed, four each the two 
cylindrical walls between which the fluid contained. The assumption that 
there slipping the fluid the walls establishes the six conditions 


The remaining pair conditions must involve the magnetic field vector. 
When the cylinder walls are non-conducting current can cross the boun- 
daries, that the radial component current density must satisfy 


must, therefore, require that 


Equation (8), along with the boundary conditions (9) and (10), forms 
eigenvalue problem for the parameters and 


(c) The Solution 

Chandrasekhar has frequently used the calculus variations order 
circumvent the tremendous labor solving high order linear equations with 
constant coefficients conventional methods. implement the variational 
procedure necessary derive integral expression for one the para- 
meters (in this case which yields minimum value when the equation 
governing marginal stability satisfied. For the cylinder problem with 
perfectly conducting walls Chandrasekhar (1953) proved that the expression 


e+} 

can meet this requirement. Here 
can verified that the same expression valid for the case non- 

conducting walls, when the boundary conditions are equivalent 


The problem thus differs from that Chandrasekhar only that the boundary 


tegration. 
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Following Chandrasekhar (1953) express and terms 
variational parameter that 


and 


both which vanish conformity with the boundary conditions. 
Then must satisfy 


where 


and are constants integration. The boundary conditions 


From (17) and (18) have 


Both numerator and denominator the integral expression (11) can 
written terms the variational parameter using (14), (16), (19), and (20). 
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The numerator (11) identical with that derived Chandrasekhar 
which 


where 


and 

After integrating and arranging the terms, find for the denominator 

where 


iff3 


and 
= — 3”) 
FP 
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Thus (11) becomes 


(23) 


The value for which minimum the root the equation 


(24) 

Critical Taylor numbers have been evaluated for seven values ranging 
from zero 6400. this was calculated from (24) for each suitable 
series values The minimum value i.e. was found plotting 


TABLE 


TAYLOR NUMBERS, AND WAVE NUMBERS THE MODES 
WHICH APPEAR THE ONSET INSTABILITY FOR THE CASE WHEN THE 
TWO CYLINDERS ROTATE THE SAME SENSE AND 


T./To 
2.565+ .01 4.489 108 2.628 
100 1.69 +.01 1.082 « 104 6.335 
400 0.775+.01 4.273 
1600 0.39 +.01 1.70 X105 99.53 
6400 0.19 +.01 6.85 
*The limits indicate the accuracy to which the wave numbers could be determined 
graphically. 
7 
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Fic. instability curves. The critical Taylor numbers functions 
for the cases perfectly conducting and non-conducting cylinder walls. 

Fic. Wave number the mode which first develops the onset instability 
function 
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The dimensionless wave number equivalent 
The value which corresponds critical Taylor number defines the wave- 
length the mode which develops the onset instability. The values 
the critical Taylor numbers and wave numbers are given Table and 
curves for perfectly conducting walls are shown for comparison. Since the 
two cases represent extreme conditions far the boundaries are concerned, 
reasonable suppose that they provide upper and lower bounds the 
general solution which the boundaries have finite electrical conductivity. 
must remembered that these analyses are valid only when the separation 
the cylinders small compared their mean radius, and when both 
rotate the same direction. 

would expected, the present analysis gives the same result 1708) 
the field-free case was obtained Chandrasekhar (1953) and Taylor 
(1923). seen Figs. and that the effect non-conducting walls 
reduce the critical Taylor number factor nearly five when 
and increase the wavelengths the instability patterns considerable 
amounts. surprising feature analysis that the wave 
number first increases with the expected decrease not taking place till 
50. When the walls are non-conducting, however, the wave number 
decreases steadily increases from zero. 


EXPERIMENT 

(a) The Magnetic Field 

When the experiments were attempted, the solution the previous section 
for non-conducting walls had not been worked out. However, was evident 
from theory that field intensities about 5000 oersted 
would required demonstrate the inhibition instability mercury 
between two cylinders moderate dimensions. large rectangular water- 
cooled solenoid was wound from heavy copper conductor produce fields 
this size. The finished coil consisted five layers windings, each layer 
having turns. was approximately length with inside diameter 
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Fic. magnetic field along the axis the solenoid. 
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The solenoid was designed operate from d-c. motor-driven generator 
capable delivering 2100 amperes approximately volts. The power 
output was continuously variable from zero the maximum about 
25,000 watts. The curve plotted Fig. shows the measured field variation 
along the axis the coil. The field was reduced only distance 
either side the center, but fell off rapidly greater distances. 


The Cylinder Apparatus 

The inner and outer cylinders are shown separately and assembled Fig. 
Only the inner one was able rotate. They were made from perspex the 
following dimensions: 


inner cylinder 


radius central portion 1.740 cm, 

radius end portions 2.012 cm, 

over-all length cm, 

length central portion 22.0 cm, 
outer cylinder 

internal radius 2.554 cm, 

over-all length 
annular space between the cylinders 

width central portion 0.814 cm, 

width end portions 0.452 cm, 

volume 376 cc. 


The inner cylinder was made larger its ends effort counteract the 
effect the reduced magnetic field points remote from the center the 
solenoid. While the effect the reduced field was enhance the onset 
instability, that the reduced gap was delay it. Thus the cylinders were 
designed prevent, far possible, the premature development tur- 
bulence near the ends which might upset the laminar flow the remainder 
the fluid. 

The cylinders were mounted dural base plate that the assembly 
could rigidly supported inside the solenoid. The heavy brass flywheel 
(Fig. attached the lower spindle the inner cylinder served stabilize 
the drive, and also give the cylinder sufficient downward thrust overcome 
the buoyancy force when the apparatus was filled with mercury. small 
phonic motor bolted the floor was connected the flywheel shaft through 
reduction gear and flexible rubber coupling. The speed the inner cylinder 
could varied continuously from r.p.m. varying the frequency 
the audio input signal the motor. 


(c) The Tracer 

Since all liquid metals are opaque, their motions are not easily traced 
means dyes suspended particles. this experiment the use the radio- 
active isotopes and conjunction with suitably placed geiger 
tube was attempted. has half-life hours and emits gamma 


t 
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Fic. The perspex cylinder apparatus. 


radiation 0.077 and 0.191 Mev, while Hg?® has half-life days and 
emits 0.279-Mev gamma ray. small amount these isotopes carefully 
introduced the middle the inner cylinder would expected remain 
adjacent that cylinder long the motion was laminar. the onset 
instability, the development radial component fluid velocity would 
cause some the radioactive material carried across the annular space 
the wall the outer cylinder. Provided the radiation efficiently absorbed 
the intervening layer mercury while the motion laminar, the onset 
instability would accompanied large increase the counting rate 
geiger tube mounted near the outer wall. 

data were located the absorption coefficient mercury for gamma 
rays the required energy range (0.2 0.3 Mev). Radiation this energy 
should obey approximately the absorption law 


where the initial radiation intensity, the intensity after absorption, 


? 
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the mass absorption coefficient cm? per the thickness the 
absorbing medium per and its density. The value required 
reduce the intensity factor two 


Rough experimental determinations the absorption were made using 
the source, and ordinary mercury between the cylinders the 
absorbing medium. The source was placed sealed specimen bottle which 
was lowered down the hole the axis the inner cylinder (Fig. 4). The 
activity was measured means ring-shaped geiger tube which could 
mounted about the outer cylinder. Counting rates were registered scaler 
when the absorbing layer was present and again when was absent. The 
results three such experiments, after allowance had been made for the 
background counting rate, are given below. 


Experiment counts/minute counts/minute 
10055 0.43 0.12 
141 11059 0.39 0.13 
12147 0.44 0.12 


Since the layer mercury between the cylinders was 0.8 thick, was 
reasonable suppose that the shielding would sufficient enable one 
distinguish between laminar and unstable motion the radioactive tracer 
method. Furthermore, fresh isotope mixture most the activity arose 
from which emits softer gamma ray than does absorption 
the perspex parts the cylinder apparatus was tested and found 


negligible. 


(d) The Experimental Procedure 

The radioactive mercury, obtained from the Atomic Energy Research 
Establishment Harwell, usually consisted 2.4 millicuries and 
120 millicuries mixed with the natural element. The mixture was 
dispensed with micropipette capable transferring measured volumes 

The radioactive mercury was introduced through glass capillary tube 
mounted hole drilled along the axis the inner cylinder (Fig. 4). 
small radial outlet was drilled from the outside the cylinder the bottom 
the capillary housing. The capillary tube projected about inches above 
the top the cylinder, the bore the upper end being enlarged form 
small funnel. subsidiary wide-bore glass tube containing tap was cemented 
the upper cap the cylinders (Fig. 4). 

Mercury from external reservoir was run into the apparatus very slowly 
ensure that all the air would forced through the outlets the top. The 
mercury was allowed rise till filled the cylinders and the capillary tube 


0.693 
g/cm’. 
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the level the bottom the small funnel. The tap the subsidiary tube 
was left open that mercury rose slightly higher level. measured 
volume the isotope mixture was then transferred the funnel that 
lay top the column mercury. The original level was then restored 
that the upper part the capillary column consisted the radioactive mixture. 
Then sufficient natural mercury was added the funnel that when the level 
was restored second time the radioactive material was the bottom 
the capillary and ready injected through the radial duct into the fluid 
between the cylinders. 

Power was supplied the solenoid, and the field adjusted till its value 
near the middle the cylinders was considered well above the critical 
value required preserve laminar motion. After this had been done, the inner 
cylinder was set rotating speed which was held constant throughout the 
experiment. The geiger tube was positioned opposite the radial duct, and the 
counting rate was noted this stage. The mercury levels were reduced such 
rate that flow the glass columns was barely discernible. The isotopes were 
thus forced low speed into the main body mercury while the inner cylinder 
was rotating. Great care had taken prevent the mixture from travelling 
across the annular space the outer wall the beginning. The field was then 
reduced zero steps about 350 oersted, interval seconds being 
permitted elapse each step during which the counting rate was noted 
rate meter. This procedure was called the run’’. Immediately after 
this series observations the field was restored its initial value and the 
process repeated for “blank run” which further radioactive mercury 
was added. The purpose the second run was ensure that the variations 
counting rate, time when the isotopes were more less randomly 
mixed with the main body mercury, were appreciably smaller than those 
which occurred during the live run. this were not so, sometimes happened, 
significance could attached the live run, and the results were discarded. 

Contaminated mercury which was drained from the apparatus was usually 
fit for use again after week two owing the rapid decay More- 
over, after the isotopes had been thoroughly mixed with the rest the fluid 
the background count was considerably reduced because the high ab- 
sorption coefficient mercury. 


(e) Experimental Results 

Experiments which significant rise the counting rate was observed 
during the live run were completed for five different speeds the inner cylinder. 
The observations and relevant data one these experiments 11.9 
r.p.m.) are given Table II, and the rate-meter deflections are plotted against 
magnetic field strength Fig. The data from all five experiments are too 
numerous record. Different runs which were done the same rotation 
speed are labelled etc. the rate-meter readings indicate, the back- 
ground count most the experiments was quite high. Many the runs, 
expecially the earlier ones, were unsuccessful and the high cost mercury 
prohibited the use uncontaminated fluid every run. The background 


H 
i 
| 
i 


1520 CANADIAN JOURNAL PHYSICS. VOL. 36, 1958 


TABLE 


OBSERVATIONS IN THE FIRST CYLINDER EXPERIMENT 


Temperature 


6.37 


Rate-meter Rate-meter 
Field readings readings 
intensity, live run, blank run, 
oersted 
Run 
3640 
3295 
2770 
2420 
2080 
1725 
1380 
1035 
690 
345 
Run 
3990 
3640 
3295 
2770 
2420 
2080 
1725 
1380 
1035 
690 
Run 
3810 
3470 
3470 
3370 
2940 
2600 
2245 
1900 
1550 
1210 
960 
520 


RATEMETER DEFLECTION 


Observations the first experiment. Rate-meter deflections plotted 
against magnetic field intensity. 
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count was measured with scaler over several 1-minute intervals before 
experiment was begun make sure that was constant. The time constant 
the rate-meter tank circuit was about seconds, that the ammeter 
needle was not subject rapid fluctuations. Whenever the rate meter in- 
dicated significant change counting rate given stage, the reading 
recorded was that which corresponded the end the 30-second period. 

was assumed that the critical magnetic field, was higher than the 
value which increase the counting rate was first observed, but lower 
than that the previous stage. The values and estimated from 
the five experiments are given Table III. The experimental results are 
plotted Figs. and where they may compared with the theoretical 
curves. 


TABLE III 
oersted Adopted 

11.9 3030 3470 3420 3420 
2945 2695 2695 2780 3.43 
22.7 10° 3295 3210 3210 4.56 
26.7 3640 3640 3640 3640 5.88 
29.7 3985 3985 3985 7.03 
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Fic. Critical speeds the inner cylinder plotted against magnetic field intensity. 


(f) Discussion 

While the experiments were progress was not realized that non- 
conducting walls could have such large effect indicated the theory 
the preceding section. The theoretical curve for non-conducting walls 
plotted Fig. shows that the motion should have been unstable the 
outset all the experiments except that for which the speed the inner 
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Fic. instability curves and experimental results. 


cylinder was 11.9 r.p.m. The following table gives the initial field intensities 
the experiments, and the approximate critical values based the theories 
conducting walls and non-conducting walls. 


Rotation Initial field, oersted 
speed, Conducting Non-conducting 
Experiment r.p.m. walls walls 
11.9 3640 3990 3810 1900 3710 
17.9 3810 3740 3740 2775 5600 
22.7 4160 4160 3500 7100 
26.7 4850 4850 4850 4100 8500 
29.7 4850 4850 4550 9500 


the onset instability was expected that cellular patterns would 
develop manner similar that observed Taylor (1923), i.e. the motion 
would and down near the wall the inner cylinder, radially outwards 
the cell boundaries, and down and along the wall the outer cylinder. 
the outer wall the upward and downward streams would meet the middle 
where radially inward motion would develop complete the pattern. 
The ‘‘wavelength” cell the distance between two radially outward 
streams. Whereas this distance was found Taylor very nearly equal 
twice the separation the cylinders the field-free case, theory indicates 
that the cells are, general, elongated when the motion influenced 
axial magnetic field. The wavelengths the cell patterns are plotted against 
critical rotation speed for the cases perfectly conducting walls and non- 
conducting walls Fig. The curves apply only cylinders having the same 
radial dimensions those used this series experiments. The two curves 

cet when the critical speed 0.615 r.p.m., corresponding the field-free 
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case. seen that the wavelengths increase far more rapidly when the walls 
are non-conducting. 


NON-CONDUCTING WALLS 


2n 


PERFECTLY CONDUCTING 
! 


° 10 20 x” 
CRITICAL ROTATION SPEED (rem) 


Fic. wavelengths the cell patterns the onset instability. 


The fact that unstable motion was not detected the outset the last 
four experiments can explained the long wavelengths associated with 
the instability patterns. rotation speeds 17.9, 22.7, 26.7, and 29.7 r.p.m., 
the predicted wavelengths are and from 
theory, and approximately 40, 49, 58, and when the walls are non- 
conducting. Thus these experiments the upward and downward flow along 
the inner cyclinder would carry the isotopes into the region the reduced 
gap, and right the ends the cylinders. Evidently the speed flow 
the start the runs was too slow permit the isotopes travel this far 
the seconds during which the field was held constant. Whether not 
the convection could develop properly uncertain, because, except the 
second experiment 17.9 r.p.m.), there was not sufficient space 
contain single complete cell pattern. the second the pattern may have 
formed, but the motion was certainly disturbed end effects. The results 
indicate that some the isotopes eventually reached the vicinity the 
geiger tube round-about path, the average travel time the four experi- 
ments being about two minutes. 

The travel times (in seconds) observed the last four experiments are 


outlined below: 


Run 2nd expt. 3rd expt. 4th expt. 5th expt. 
120-150 120-150 


While the observations not allow very accurate determination the times, 
clear that all the runs attempted the interval between injection the 
isotope and the increase counting rate close two minutes. Since the 
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speed transport should not have varied greatly from one experiment 
the next, and since the path length was the same each (about cm), 
constant travel-time expected. 

Thus the rise counting rate observed the last four experiments 
attributed long travel-time the isotopes, and not the onset instability 
the corresponding value the magnetic field. The fact that the results 
lay fairly smooth curve (Fig. 6), and did not depart more than 15% 
from the critical fields estimated from Chandrasekhar’s theory now con- 
sidered fortuitous. Different results would almost certainly have been 
obtained different time intervals had been used between the stages, 
different initial fields had been applied. 

The first experiment 11.9 r.p.m. gave result which was incompatible 
with the other four, and which lay close the curve predicted for non- 
conducting walls (Figs. and 7). this speed the predicted wavelengths 
are 4.3 for perfectly conducting walls, and 25.7 for non-conducting 
walls. Thus single cell pattern would, theory, extend about beyond 
each end the central portion the inner cylinder. reasonable sup- 
position that the elongated cell pattern did form during these runs, because 
run when the motion should have been unstable from the start small 
margin, the increase counting rate was not observed till seconds after 
the injection the isotopes. both runs and the motion should have 
been stable the initial stage and unstable the second. the increase 
counting rate came seconds after the beginning the second stage, 
while came seconds afterward. view the reduced wavelength, 
the isotopes were transported roughly the same velocity (0.3 cm/sec) 
they were the other experiments. 

The first experiment supports the theory the foregoing section. The 
adopted experimental value the critical axial field was 3420 oersted whereas 
theory predicts 3710 oersted. Allowing for time lag about one 
minute between the onset instability and the increase counting rate, 
can inferred that lay between 3900 and 3640 oersted run and 
between 3810 and 3470 oersted run The agreement good can 
expected from the technique employed, but the evidence not conclusive. 
The data obtained the first experiment not, unfortunately, establish 
beyond all doubt that the motion was laminar the start runs and 
However, rotation speed 11.9 r.p.m. the critical field the case 
perfectly conducting walls 1900 oersted. The experiment does, therefore, 
clearly indicate that when the cylinders are made from non-conducting 
material the onset instability much less inhibited magnetic field. 

Taylor’s experiments, when the wavelengths were close cm, the time 
taken for the colored liquid trace out complete cell pattern was 
seconds. Thus was originally considered that constant, nearly constant, 
counting rate over the initial seconds run would conclusive evidence 
stable motion. However, the sevenfold increase wavelength the in- 
stability patterns predicted the theory for non-conducting walls makes this 
assumption invalid. the first experiment, the initial field exceeded the 
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theoretical value 280 oersted run and 100 oersted run 
The margin was very small, and the observations should have started 
higher fields and extended over time intervals about minutes establish 
that the motion was initially stable. 

When the discrepancies between theory and experiment were realized, 
was longer possible run another series experiments. The apparatus 
described herein would probably satisfactory rotation speeds r.p.m. 
less longer time intervals were allowed between reductions the field. 
The use longer cylinders and larger magnetic fields indicated higher 
speeds. Experiments would then difficult because the greatly increased 
requirement for d-c. power. 
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ROTATIONAL ANALYSIS THE SYSTEM THE 
MOLECULE! 


SURYANARAYANA RAO? 


ABSTRACT 

The bands the system the molecule have been photographed under 
high dispersion (0.35 A/mm). rotational analysis the 0-0, 0-1, and 1-0 
bands given, which differs from the one previously given Sen Gupta. ad- 
dition, four more bands, namely, the and bands, have been 
analyzed. The bands are attributed the electronic transition, 
the lower state being the ground state the molecule. The new rotational con- 
stants for the ground state are the following: 


The spin doubling the upper state small. Perturbations the level 
the upper state, which were not reported previously, are observed and dis- 
cussed. They supply welcome confirmation the ccrrectness the analysis here 
presented. 


INTRODUCTION 


Two band systems the molecule are known, the system extending 
from 3200 3600 and the system extending from 2800 2300 These 
bands have been described Petrikaln (1928), Curry, Herzberg, and Herzberg 
(1933), Ghosh and Ball (1931), and Dressler (1955) and rotational analysis 
the system has been given Sen Gupta (1935). Recent work Singh 
(1957) the system has indirectly cast considerable doubt this rotational 
analysis, and the present work was undertaken check directly the analysis 
the system. The 0-0, 1-2, 2-1, 2-3, and bands have been 
photographed high dispersion and analyzed. Aided perturbations not 
previously found Sen Gupta, the new analysis shows quite definitely that 
Gupta’s analysis incorrect. The constants the upper and lower states 
have been determined from the new spectra. 


EXPERIMENTAL 
The spectrum was excited high frequency discharge tube through 
which helium containing trace oxygen and phosphorus vapor was 
pumped. 100-watt, 2450-Mc/sec oscillator was used produce the dis- 
charge and the spectra were photographed the third order 10-meter 
concave grating spectrograph. Iron lines obtained from iron hollow cathode 
discharge tube were used wave-length standards. The error the relative 
wave-numbers the lines estimated be, general, less than 0.03 
Under the conditions described, the band extremely strong 
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and could photographed with very short exposures. One the subbands 
the band and the corresponding subband the 2-1 band 
are shown Fig. Bands involving higher vibrational quanta are weaker 
and overlap each other. Also, with long exposures, bands and 
appear and become troublesome. result this overlapping only 
relatively small number the bands could analyzed and even these 
bands (except many lines are overlapped. 


ANALYSIS AND DETERMINATION CONSTANTS 


Each band consists two subbands, with each subband showing two heads. 
Eight branches could picked out each band. Their J-numbering has been 
determined the usual way. These features suggest that the electronic 
transition the bands The nature the electronic transition these 
bands will discussed detail later. The number expected branches for 
this transition 12, six them being the main branches 
and and the other six being satellite ones, namely, 
and The analysis shows that the spin-doubling the upper state 
extremely small and view this, four the expected satellite branches, 
and are superposed the main branches. The wave- 
numbers the identified lines the seven bands are given Table 


TABLE 


band 

62.42 71.62 82.62 94.50 17.78 
60.74 71.62 84.04 94.50 20.77 

11.92 
11} 56.60 71.96 11.00 79.85 14.30 34.11 
53.71 76.93 
12} 54.23 72.27 92.47 14.3 77.59 16.23 
14} 51.82 73.52 75.50 44.54 
25.85 
72.58 

234 19.52 58.28 07.50 


| 
5 
3 
t 
t 
| 


90.22 

92.12 

94.00 

39149.72 
44.07 
40.77 

32.72 

31.92 39154.85 
31.18 

30.02 
29.28 
29.28 

76.68 
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TABLE (continued) 


Q2(J) 


40449. 


15. 


72.07 


0-1 band 


74.73 
74.73 
62.21 74.73 
61.08 74.73 


57.44 76.46 


1-0 band 


80.93 


59.14 


84.07 
85.54 
87.09 
88.78 
89.40 
92.04 


94.90 


93.13 
98.22 


94.90 
97.91 

401.04 
04.01 


14.89 


18.39 
22.19 


31.01 
34.77 


52.34 


1528 
254 25.49 10.88 51.42 
31.86 14.71 58.34 
284 16.77 61.91 
294 38.58 18.90 65.61 
42.06 21.15 69.42 
53.09 28.61 81.40 
31.28 85.53 
60.91 34.04 89.91 
64.94 
37} 98.75 
393 708.09 
39161.14 
62.55 
64.05 
67.15 84.78 
103 
13} 75.33 98.52 57.44 77.47 
14} 56.44 78.26 
15} 79.35 55.93 79.14 03.71 
16} 83.08 10.70 55.38 80.05 06.31 
82.00 
87.78 18.84 54.74 82.72 12.33 
55.38 
90.47 23.19 54.74 84.07 14.89 
24} 05.55 92.47 
12.32 
293 
78.60 
79.73 
80.97 
82.33 
83.76 
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TABLE (continued) 


Q2(J) 


65.37 56. 88.63 

13} 70.64 93.77 52.64 72.54 94.25 
14} 27.79 72.54 51.52 73.05 
41776.59 41804.39 41949.76 41974.36 42000.70 
24.71 50.93 78.72 07.99 49.03 75.11 02.99 
18} 51.61 80.93 11.85 75.98 05.41 

52.38 83.24 15.61 47.74 76.93 

203 53.23 85.63 19.57 

21.7 54.16 23.61 79.08 

233 56.21 31.35 81.60 

28} 62.72 07.34 

91.04 

303 
94.60 
i 1-2 band 

50.30 
29.79 52.74 
31.09 54.09 
12} 04.63 23.12 43.19 46.78 
13} 23.71 45.30 69.16 
02.04 25.18 49.88 74.01 
01.38 52.28 25.71 50.30 76.62 
00.86 .06 54.74 51.41 79.33 
18} 00.39 28.13 24.96 52.74 

35.08 72.11 60.76 97.79 

64.80 04.86 

273 42.44 12.30 

71.29 
29} 73.78 
2-1 band 

4 86.12 99 93 22.82 


TABLE (continued) 
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79.12 05.18 02.18 28.38 43.70 
77.57 06.70 41923.61 00.85 29.99 
103 76.10 26.65 099.49 42114.65 
74.68 09.95 29.88 14.96 
73.38 41891.73 11.71 33.20 15.33 
72.15 92.03 13.54 36.53 15.79 60.24 
71.00 92.43 15.45 39.96 94.93 16.29 
69.93 92.91 17.46 94.00 16.89 
67.96 94.05 21.64 18.36 
64.62 98.31 33.54 23.59 58.94 
233 64.18 99.40 35.74 24.91 61.79 
243 63.81 900.58 64.73 
63.48 01.81 41.61 67.74 
63.27 03.14 44.46 29.36 70.84 
06.04 50.41 32.80 
293 07.60 53.49 34.68 80.71 . 
56.64 36.57 84.17 
10.96 59.90 38.61 
14.64 42.91 
16.61 70.15 45.18 
18.65 
363 20.74 50.02 
373 52.53 
383 25.14 55.08 

2-3 band 
75.13 92.01 10.39 30.30 98.81 15.65 34.18 
74.09 92.50 12.41 33.89 16.27 36.31 
73.15 14.52 37.54 97.11 16.98 38.57 
143 72.30 93.76 16.76 41.32 17.81 40.94 
71.60 94.57 19.09 18.74 
163 70.94 95.46 21.51 
173 70.37 96.47 24.04 95.30 20.95 
19} 69.58 98.73 54.48 
69.58 32.23 25.07 57.43 
214 01.48 35.16 60.53 
02.98 38.19 28.38 
233 04.59 41.32 30.20 67.14 
06.30 32.13 70.62 
10.02 36.31 
273 12.04 38.57 81.68 
14.12 40.94 85.56 
293 16.36 89.56 
303 39518.67 39745.99 
48.68 
23.59 


$ 
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TABLE (concluded) 


SOMO 


~ 


299 
96.11 
95.10 
94.19 
93.40 
92.73 
91.71 
91.32 
91.16 
91.16 
91.16 


11.02 
11.42 
11.95 
12.59 
13.34 
14.19 

16.23 
18.73 


21.68 


2-4 band 
18.36 
22.81 23.27 
24.53 21.97 
20.86 
19.97 
30.36 19.07 
32.49 18.19 
34.77 
42.28 16.04 
47.81 
50.76 
53.85 
56.99 
60.28 
63.68 


54.58 
57.43 
60.53 
41.57 
44.57 
46.24 38560.01 
34.56 66.77 
35.05 52.00 70.19 
35.69 54.19 73.98 
77.79 
37.31 58.84 81.72 
85.43 
38589.98 
66.77 94.35 
41.97 69.49 98.68 
45.20 75.79 
46.77 79.04 
82.41 
50.55 
52.61 89.51 
54.81 
59.60 
05.15 
64.83 09.38 
70.54 
73.63 
76.78 


The following are the wave-numbers the satellite branch lines the band separated 


from the main branch ones 


40385 


85.97 153 


401.18 


97.61 


12.34 


1531 
343 
28.89 
30.99 
33.21 
35.50 
40.50 
303 
333 80.05 
343 
354 87.01 
90.64 
94.35 
98.18 
73.52 
10} 
76.23 
90.83 
123 92.67 
94.78 
144 
154 
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The rotational constant, the level the lower state has been 
determined from combination differences plotting the mean value 
other vibrational levels the lower and the upper states have been determined 
from and from the values the bands obtained from their 
branches according the equation given below for the determination the 
band origins. Thus, using the values determined indicated above, 
the and values have been determined from the band 
and the 1-0 band respectively. Using the thus determined, 
the has been obtained from (Bj—BY) the band. Similarly all the 
other values have been obtained. This method has been used makes 
possible get higher relative accuracy the values and also some the 
bands are not developed well enough give sufficient number com- 
bination differences for direct determination the values. 

The origins the bands have been obtained using the relation 
yields straight line whose intercept with the ordinate axis gives 
2vo+2B"’. 

The spin-doubling the lower state, has been evaluated from the 
relation 


where A/B,. Experimental values (J) can found from 
against straight line obtained whose intercept the ordinate 
axis gives (Y—4). 

the transition the A-type doubling the lower state, pro- 
duces doubling the rotational lines but only combination defect, that 
is, the difference between and The mag- 
nitude the combination defect each band comparatively large one 
the subbands while very small the other. Since state the 
A-type doubling larger the than the component for low 
values, the subband which the combination defect larger has associated 
with the and the other, with the state the A-type 
doubling is, according Mulliken and Christy (1931), 


(J+1/2) 


positive. From the observed branches, 


Since the A-type doubling varies very little with the vibrational quantum 
number, the experimental values all the bands have been plotted 
the same graph, against The slope this straight line gives the value 


| 
| 
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The spin-doubling the upper state, could not determined the 
resolution not adequate. 

All the molecular constants determined indicated above are given 
Tables III, and IV. 


TABLE 
Molecular constant Lower state Upper state 
0.7303 0.7775 
0.7248 0.7719 
0.7192 0.7666 
0.7140 
0.7084 
224.03 
224.48 
0.0073 
TABLE III 
Band origins (in 
41862.71 
TABLE 
1391.16 6.99 0.7801 0.0054 
DISCUSSION 


The fact that there measurable combination defect proves that the 
transition either But analysis shows that the doublet- 
splitting the lower state large (224 and that negligibly small 
the upper state. Thus must conclude, agreement with previous in- 
vestigations, that the bands arise transition. Since the 
subband lies the shorter wave-length side the the lower state 
regular. Thus the bands are due transition, where the 
lower state the ground state the molecule. 

the two lowest excited states are the state 43966 and 
the state 45918 and transitions from these the ground state 
give rise the and systems respectively. the two known 
excited states are 30846 and 40808 The recent work Singh has 
shown that the upper state the bands state. The present work has 
confirmed the fact that the upper state the system also state. 
Furthermore comparison the combination defects the two band systems 


4 
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shows that both states have the same symmetry and presumably are both 
states. The two known excited states therefore not correspond 
the and states NO. 

The levels the vibrational level the upper state the 
bands are perturbed and 15}. The levels are perturbed 
and 17}. These perturbations consist abrupt increase the 
separation successive members the branches and the appearance 
extra lines. They can clearly seen the spectrogram reproduced Fig. (a) 
the subband the band comparing with Fig. (6), which 
reproduction the unperturbed subband the band. Where 
the perturbations occur, the satellite lines are separated from the lines the 
main branches. Similar perturbations are observed the band. Fig. 
2(a) and are plotted the differences between the observed and the cal- 
culated term-values for the two component levels (J) and (J), the 


level the state. has not been possible from the observed per- 
turbations determine the nature the perturbing state. But may 
remarked that, calculation, has been found that the level the 
upper state the bands almost the positions perturbations, there 
are rotational levels close energy the corresponding levels the 
level the upper state bands. This suggests that the observed 
perturbations might due the upper state the system. this were 
state might possible account for all the perturbations observed, 
two the and two the levels. But Singh’s analysis the bands 


— 


+1.6 
(a) 
° 
Or 
° ° 
| 
-0.4 
-0.8 
° 
Sz 105 205 255 
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shows that their upper state state and not state. mutual 
perturbation two states, not possible, because the selection rules, 
that each component level can perturbed more than one value, has 
been found here. possible, however, that some the perturbations are 
due the upper state the system and the others another state. 

The ground state, the molecule separated from the perturbed state 
about 41,000 unlikely that the ground state causes the per- 
turbations since the dissociation energy the molecule the ground state 
about 44,000 Further, found that the shift energy caused 
the perturbations goes back zero higher values. Thus does not seem 
very likely that the observed perturbations can heterogeneous character. 
The possibility quartet state the perturbing state cannot ruled 
out, since such intercombination perturbation would give rise pertur- 
bations similar those observed. 

The rotational constants derived for the ground state the present 
analysis agree fairly well with the preliminary values obtained Singh from 
his analysis the system bands but not with those given Sen Gupta. 
Further, the perturbations observed the level the upper state give 
unambiguous numbering the lines. Moreover, many wave-numbers 
the individual lines given Sen Gupta not agree with those the present 
investigation. Therefore appears that the analysis given Sen Gupta 
incorrect. 
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PROPERTIES CADMIUM SULPHIDE PHOTOCONDUCTIVE 


WILLIAMS 


ABSTRACT 


Measurement the noise spectra cadmium sulphide photoconductive cells 
from has failed reveal transition from the 1/f the 1/f? 
noise frequency dependence. This clearly demonstrates complete lack de- 
pendence the 1/f noise spectrum the photoconductive time constant, which 
approximately milliseconds. 

noise correlation has been observed, indicating minority carrier lifetime 
less than seconds. 

The extreme difference between the minority carrier lifetime and the photo- 
conductive time constant explained assuming that most the recombina- 
tion takes place the surface the film, but that the rate recombination 
controlled diffusion. Slow diffusion the minority carriers consequence 
their very short lifetimes. 

The noise power and time constants CdS have been observed depend 
the value the series resistance used with the cell. The two effects are believed 
related and arise from the presence intercrystalline barriers. 


INTRODUCTION 


Noise Spectra 

Two types excess noise spectra have been observed for semiconductors. 
These are and 1/f noise. 
noise has spectral form identical with the frequency response the 
cell, and predicted the basis the random generation and recombination 
electrons and holes semiconductors. This type noise will always 
present even ideal semiconducting materials (Petritz 1956). 

Noise the 1/f form has been observed for most semiconducting cells and 
many cases completely dominates the g-r. form. appears that the 
1/f noise related the surface properties semiconductors (Kingston 
1956; Morrison 1956). However, Hill and Van Vliet (1958) have shown that 
recombination controlled surface states does not necessarily result 
dominance 1/f noise. Van Vliet and Van der Ziel (1958) have also shown that 
any process that depends diffusion will not give rise 1/f noise. Possible 
sources 1/f noise are discussed McWhorter (1956). 

Since cadmium sulphide cells are polycrystalline, 1/f noise could result 
from the passage current through barriers junctions controlled 
manner. Gianola (1956) observed 1/f spectrum when current was passed 
through crystal diode. This follows the type argument suggested 
Petritz (1957). 

Hyde (1956) reported transition from spectrum, and 
occurred near frequency equal the reciprocal times the photo- 
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conductive time constant, The noise measurements were not accurate 
enough ensure that the transition from spectrum was real. 
lack any dependence the 1/f noise spectrum has been demonstrated 
the present work. fact, transition was observed the 
frequency region cycles/sec Mc/sec. The photoconductive time 
constant the CdS studied approximately milliseconds. 

the course the investigations, was noticed that changes the magni- 
tude the load resistor affected the noise current unexpected manner 
well the photoconductive response time. This effect has been investigated 
detail. 

means correlation experiments, »mery (1952) appeared 
establish that could related the The noise voltage observed 
one point along germanium transferred part second 
point near with time delay predicted from the known minority carrier 
drift mobility. similar type measurement has been performed with CdS 
ceils the course the present investigation. correlation was observed 
with these cells. 

The resistance the CdS cells investigated was usually too great permit 
accurate measurements over the whole frequency range from cycles/sec 
Mc/sec when the cell was darkness. Practically all the noise measure- 
ments were therefore made with the cell resistance reduced radiation from 
small tungsten lamp. series noise and time constant measurements were 
performed elevated temperatures see the thermal and radiation 
effects are different the cells being tested. Van Vliet (1956) has shown 
that with single crystal cadmium sulphide the wavelength the radiation 
falling single crystals affects the time constant and the noise, but only 
small degree. 


Cadmium Sulphide Cells 

The cadmium sulphide cells used were manufactured the Canadian 
Marconi Co. The electrical contacts were indium and proved give good 
ohmic behavior. review the characteristics similar cells given 
Bube (1956) and Kroger al. (1954). The cell characteristics given the 
manufacturer are: electron mobility, approximately 300 and time 
constant from 100 milliseconds. 

The cells are multisectioned cells, shown Fig. For the present 
experiments, the individual sections cell were separated scraping off 
the indium contact between successive cell sections. Electrical contact was 
made the indium placing fine wire against and then covering the 
wire with silver paint. The cell and contacts had linear 
characteristic least the experimental accuracy 4%. photovoltaic 
effects could observed with these cells, even for light signals which changed 
the average resistance factor 10. Figure illustrates the form used 
for noise measurements, while Fig. shows the arrangement used when 
correlation effects were studied. 
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area 


ELvectaicaL LEAD 


BAKE AREA 


cds 


(c) 


Fic. CdS cells showing their structure and method making contact. 


NOISE SPECTRA 


Technique 

The noise spectra CdS cells have been measured from cycles/sec 
Mc/sec. The equipment used consisted units. the region from 
cycles/sec kc/sec General Radio wave analyzer was employed. The 
signal was tapped from the stage preceding the meter and passed Ballan- 
tine r.m.s. voltmeter. The average value least seconds the recorded 
d-c. mean square output the meter was used for each reading. The zero 
current noise level, the Johnson noise, was recorded for calibration the 
apparatus and enable zero current noise corrections made. the 
square the noise voltage was being recorded, only simple subtraction 
the zero current noise from the recorded value was necessary. 

The high frequency region was covered two specially constructed wave 
analyzers whose frequency ranges, kc/sec 200 kc/sec and kc/sec 
Mc/sec, were chosen overlap appreciably. The sensitivity these 


(b) 
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laboratory analyzers was least good that commercial models. all 
cases the noise was measured with constant band width. 

With these three analyzers, there was sufficient frequency overlap for the 
relative noise spectrum from cycles/sec Mc/sec drawn 
adjusting the ordinates the curves obtained with the different instruments, 
that the curves coincided the frequency overlap regions. This eliminated 
the necessity for calibration the instruments. 

the high frequency region the noise voltages all were small, and care had 
taken reduce amplifier noise. Using cascode circuit and 5’s 
input tubes, equivalent noise resistance approximately 1500 ohms was 
realized. The minimum cell circuit resistance used was 1000 ohms that 
amplifier noise was not serious handicap the measurements. 

The cells were prepared shown Fig. and only one section was in- 
vestigated time. 


Results 

Detailed frequency spectra three cells were taken along with limited 
noise measurements four others. The magnitude the noise voltage was 
comparable all cases and the noise spectra were all very similar. Cell 
B-68-6 was investigated most thoroughly and will discussed representa- 
tive all the cells. 

Figure graph the spectra obtained room temperature with various 
currents J,. The series resistance was 1000 ohms and the illumination 
falling the cell was such reduce the cell resistance from approxi- 
mately megohms 8000 ohms. The other symbols used the graphs are 
the cell voltage; the frequency operation; and the mean square 
noise voltage observed. 

Figures and show the same type measurement. Fig. the series 
resistance was increased 100 kilohms and Fig. the cell resistance 
was raised from 100 kQ, the series resistance now being megohm. 

Figures and not cover the full frequency range, since the high value 
made impossible take accurate measurements the highest fre- 
quencies. 

The curves Fig. are nearly the megacycle region, but are closer 
the kilocycle region. 

transition from frequency spectrum would indicate short 
time limit the surface trapping times (Van der Ziel 1950). Fig. 
transition occurs. Using the fact that the noise spectra not change 
such transition, some surface trapping times must less than 
seconds. 


Analysis 

analyze the spectra, data were taken from curves the type shown 
Figs. 2-4 the kc/sec region, where noise approximating the 1/f 
spectrum was found all curves. form was assumed for the 
noise current The voltage observed was related the equation 


F 
i 
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ARBITRARY UNITS 


2 


FREQUENCY 


Cell current indicated graph expressed milliamperes. 


Table for cells B-68-6 and B-68-8. 


TABLE 
ANALYSIS NOISE DATA ASSUMING FORM 
Re, Rs, To, f, A T, 
B-68-6 section 
100 100 0.9 1.6 
100 1000 2.0 0.38 
B-68-8 section 


1000 1000 


‘+, 
* 
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2 


VOLTAGE) 


wound. Cell current indicated graph expressed milliamperes. 


Figure curve for fixed values and but for variable and 
R,. The curves fact are but under the conditions measurements 
was expected that would independent but this proved 
not the case. 

For fixed cell resistance, almost independent R,, but for 
with most marked for values less than kQ. For the converse 
value which shifts higher values increased. 

When small values series resistance were used, readings were often 
difficult take many were observed. These were large 
transients which sometimes made readings impossible. the whole band-pass 
the preamplifiers were observed oscilloscope, sudden changes the 
mean noise level would observed with slow drift back the mean level. 
and fluctuations could seen voltmeter used observe the voltage 
across the cell. This suggests AR/R changes the order 1%. This very 
large comparison the AR/R changes associated with normal noise when 
AR/R the order (Brophy 1957). 


«x 
* 
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Cell current indicated the graph expressed milliamperes. 


TIME CONSTANT EFFECTS 

The photoconductive response times the cells were determined using 
light source neon bulb triggered manually operated mercury switch. 
The response time the neon was known less than microsecond, and 
that the circuit used measure time constants was approximately 
microseconds. The measurements were taken visually from oscilloscope 
with estimated errors about 20%. 

Table lists the measured response times four different cells with 
various values and R,. For each cell resistance, the series resistance was 
varied from 1000 The response time depends the series resistance 
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X=>Re VARIABLE 
MA. 


O:Rs VARIABLE 
Rez 100K 


0.2 


(ARBITRARY) 


Fic. Variation with (X) and Cell B-68-5 section 


and, within the accuracy the measurement, always increases de- 
creases, approaching approximately constant value for small The 
variation response time with much more marked for low values cell 
resistance. 


TABLE 
DEPENDENCE THE RESPONSE TIME CADMIUM SULPHIDE CELLS THE SERIES RESISTANCE 


T T 
B-68-6-8, B-68-5-3, B-68-4-3, B-68-8-3, 
sec sec sec sec 
100 kQ 1 MQ 009 O19 .032 Oll 


noted that the conditions under which both the time constant 
and vary with the value one high current. 
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CORRELATION 


Experimental Technique 

The method used measure correlation effects followed essentially that 
Montgomery (1952). The electrical circuitry associated with the cell 
indicated Fig. Two separate amplifiers were connected that the mean 
square value the noise voltage could determined for contacts (1) and 
(2), and their sum and difference (addition 180° out phase). Ballantine 
meter was used convert the signal d-c. voltage proportional the mean 
square the noise voltage. The mean square output was recorded mentioned 
earlier. 


Results 

Four cells were tested all, and case was there evidence correlation 
excess two per cent. This limit set the accuracy the apparatus. 

For low values the correlation coefficient the relationship developed 
minority carrier lifetime and the time for electron travel the distance 
between the probes, Fig. 1(c). 

Before estimating using the above equation, the physical arrangements 
the electrodes will considered. Figure indicates cross section the 
cell. 


Fic. Constant-voltage method operating cell for correlation studies. and are 
used for filtering the battery supply. 
Fic. Cross section CdS cell. 


The feature the cell considered that the electrodes are metallic 
and short out the field the film below this region. Montgomery suggested 
reduced value because the reduced field the region the contact. 
the case under consideration neglect completely the region under the 
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electrode. Since the contacts were ohmic and photovoltaic effects could 
observed with the cells, the contacts were not believed noisy. this 
the case and the electrodes are completely shorting, then electrons coming 
one the electrodes cause equal number electrons leave the other 
side the electrode. 

Assuming electron mobility 300 cm?/v-sec and using the experimental 
0.03, must less than seconds. 

detail model for the production 1/f noise lacking, the significance 
must considered. With germanium there little trapping 
carriers that the drift mobility essentially the Hall mobility. 1/f noise 
results from surface trapping carriers any disturbance resulting from 
trapping, emission, the surface can propagated down the filament 
with velocity equal the mobility times the electric field. With CdS trapping 
minority carriers the crystal prevents the drift charge density dis- 
turbance arising from surface traps. The value obtained, seconds 
less, can taken evidence very small minority carrier lifetime. The 
long photoconductive time constant then represents the majority carrier life- 
time, the time between creation electron hole pair and their eventual 
recombination, the minority carrier having been trapped very quickly. Woods 
(1957) has demonstrated that photoconductivity lead sulphide results 
from majority carrier conduction. 


EQUIVALENCE THERMAL AND OPTICAL EXCITATION 
Double time constant effects and quenching are often observed for par- 
ticular wavelength light. Van Vliet a/. (1956) have fact shown that 
light wavelength less than 5900 produces effects CdS single crystals 
different from those produced light wavelength greater than 5900 
Thermal excitation was therefore used establish that the light source was 
not responsible for some the effects observed. 


Experimental Technique 

The cells were mounted brass block whose temperature could 
controlled. The cell mount was suspended vacuum chamber whose outer 
wall was contact with constant temperature bath, normally The 
cell was heated various temperatures 170° and the noise spectrum 
and time constant were determined. The temperature used varied between 
and After these measurements had been completed the cell was 
cooled and air admitted the system maintain the cell the tem- 
perature the bath. filament lamp was used illuminate the cell, the 
radiation being admitted through sapphire window the end 
nickel tube extending outside the bath. The illumination was varied until the 
cell resistance reached one the values used the high temperature measure- 
ments: then the noise spectrum was measured again. the cell resistance 
was brought the same value both cases, the same voltages and current 
were immediately available for both noise spectra measurements. Any changes 
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the amplitude the noise voltage then indicated changes the noise 
power under the two operating conditions. Calibration the amplifier was 
not necessary the zero current noise level was used reference for com- 
paring the sets data. Under conditions high illumination, the zero current 
noise was found the same that obtained from wirewound resistor, 
equal resistance the cell. 


Results 

Figure representative the data which were obtained. general, the 
actual magnitude the noise was not very different from the value obtained 
optical excitation. However, for given resistance value, the noise voltages 
were usually slightly higher the resistance was determined heating the 
cell, rather than illuminating it. The form the noise spectra was not 
apparently changed the optical excitation. 


1.6 sec 


spectra taken with the resistance reduced thermally and illumination 


The photoconductive time for thermal excitation about one-third 
the value obtained under optical excitation conditions, 

The significance these experiments that the peculiarities the noise 
spectra cannot attributed the excitation 
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INTERPRETATION 

1/f 

the discussion which follows will assumed that 1/f noise results 
from surface effects. The points that have been demonstrated are that the 
noise and the photoconductive time constant can not related and that 
minority carrier lifetime less than seconds. addition, other infor- 
mation can inferred from the time constant milliseconds. surface 
recombination was active, for film the order thick, using the 
equation for bulk recombination dominate d/2s, 
must less than 1/50 cm/sec. The surface recombination velocity can 
much larger diffusion controls the rate recombination. will 
shown subsequently that the lack correlation, i.e. minority carrier lifetimes 
less than seconds, consistent with this picture. Assuming that trapping 
carriers the surface responsible for noise, clear that the minority 
carriers must able get the surface that trapping can take place. 

The diffusion carriers the surface depends the ability both 
carriers diffuse. The diffusion constant for both the electrons and holes 


(Rittner 1956): 
(2) 


where and are the electron and hole concentrations, and and the 
diffusion coefficients for electrons and holes. For n-type CdS which holes 
spend large part their time trapped D,. Defining the difference 
between the average density electrons and the instantaneous value 
the one dimensional continuity equation is: 

Shockley (1950) has solved this equation while considering the effects 
the surface filament the recombination holes. One the results 
his analysis was that carriers diffuse the surface piece material 
and recombine there, the decay with time density disturbance will 
exponential. For rectangular bar cross-sectional area with surface 
recombination velocity and which the diffusion coetticient when 
the exponential time constant 


Tp Tr a 


the discussion follow will assumed that the bulk lifetime 
large small the condition governing the validity the above 
The model proposed for sensitized CdS photoconductive films 
follows. the sensitization process, states which normally 
combination paths for electrons and holes become charged electrically such 
way that the capture cross section for holes increases, while that tor electrons 
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decreases. When the CdS film sensitized the hole lifetime very short and 
the electron lifetime long. The lifetime electrons, fact, becomes long 
that recombination takes place primarily surface states. The very fast 
surface states, trapping some the carriers, give rise the noise. 

the sensitization process, well the bulk recombination lifetime 
becoming long, the diffusion coefficient becomes small. This arises the 
following way. the hole lifetime short free hole will not diffuse 
appreciable distance before being trapped. Actually, the holes spend very 
high percentage their times traps, during which time they are immobile. 
The effective diffusion coefficient with traps, then reduced from the 
trap free value, the ratio the time spent free, the time until 
recombination takes place, 


(5) D, = 


For the minority carriers the sum the time spent traps and the 
time spent the valence band, ry. 

Because small, small and the time constant 
long. 

predicts value cm?/sec for D,. This implies that for the condition 
the case, since for germanium, the smallest value that can readily 
obtained 200 (Shockley 1950, 324). 

The value for very small. Using this value for and equal 
the photoconductive time constant, equation (5) predicts value 
seconds for ry. This value less than the lower limit set the correlation 
experiment but Rose (1957) has pointed out that from photovoltaic 
effects the value must less than seconds. The small value 
thus not unreasonable. 

The very short hole lifetime and long electron lifetime CdS can ex- 
plained. When chlorine added CdS the sensitization process, electrons 
from these levels attach themselves impurity sites with energy levels near 
the valence band. large number these states become negatively charged. 
With this charge the capture cross section for holes increased and that for 
electrons decreased. The charging the sites results not only very short 
hole lifetimes but also long electron lifetimes. 

Such situation has been shown exist for manganese doped germanium 
(Tyler and Woodbury 1956). 

The role diffusion the manner discussed above has been used the 
author (Williams 1958) explain the production 1/f noise when cells 
are irradiated with short wavelength radiation. 


Breakdown 


variation the time constant and current noise with the value the 


— 
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series resistance attributed the presence intercrystalline barriers. 

noted that for small when the discussed above was 
function R,. The time constants are also affected the series resistance. 
The following process suggested explain these effects. proposed that 
there are barriers present the film and that sufficient field applied 
break some these down. When R,, the breakdown barriers corre- 
sponds reduction and consequence reduces the electric field. 
The constant current condition, would thus tend quench any 
such runaway process reducing the field. Under constant voltage conditions, 
breakdown would not quenched readily. 

The form the spectra obtained Fig. not precisely 1/f. The deviation 
from 1/f ascribed the peculiarities the breakdown. This suggested 
since, high series resistances (Fig. 4), the spectrum very and, 
under these conditions, none the mentioned earlier were observed. 

The variation time constant with accord with the model suggested 
above. The time constants the short-circuit cells are longer than those the 
open-circuit cells. This can explained saying that the high external im- 
pedance quenches the runaway processes giving short response times. When 
the external impedance low the breakdown runs itself out and the time 
constant longer. 

emphasized that the effects mentioned above were general 
small but reproducible. The role intercrystalline barriers has been considered 
detail Mahlman a/. (1956). The present work indicates that the 
barriers probably exist but play only secondary part photoconductive 
films. The paper (1957) indicates that such the case for PbS. 

The results reported suggest that trapping the minority 
carrier large degree, the characteristics CdS cells. addition, 
the series resistance effects are taken evidence small internal barriers. 
The illumination used vary the cell resistance was not responsible for any 
the effects observed. 
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SPACE CHARGE RELATIONS FOR THE LOW-FREQUENCY 
SILENT ELECTRIC DISCHARGE 


JATAR 


ABSTRACT 


The characteristics hydrogen under the low-frequency 
silent electric discharge using all-glass ozonizer are examined from the stand- 
point the space charge relation and the modified Townsend 
Werner equation. concluded that the three distinct regions the charac- 
teristics correspond the space-charge-free, the space-charge-limited, and the 
steady corona regimes. 


INTRODUCTION 

Studies (Jatar 1951, the influence ionic space charge the 
Joshi effect have shown that even chlorine the conductivity and the Joshi 
effect are determined chiefly the space charge positive ions rather than 
that negative ions, the latter which may not all necessary for the 
Joshi effect. This was confirmed the observation (Jatar, press) 
large Joshi effect over wide range pressure and potential the so-called 
‘electron-free’ hydrogen, which has very low electron affinity. During these 
latter studies, the current appeared rise parabolic function the 
applied potential (Jatar 1956a), and the results could explained assuming 
space-charge-limited regime over certain range potential. Other workers 
(Ahmed and Gill 1950; Deb and Ghosh 1948; Haries and Engel 1951) have also 
assumed the existence positive ion space charge explain the production 
high frequency pulses under the silent electric discharge and the Joshi effect. 
Recent studies (Jatar Jatar and Sharma 1956) have, however, shown 
that for potentials near the threshold potential, space charges are not 
active. was therefore interest investigate the time which the space- 
charge-free regime changed over space-charge-limited regime 
potential increased. 


EXPERIMENTAL 
The general experimental arrangement was essentially similar that 
employed earlier (Jatar 1957) and shown Fig. The hydrogen used 
was prepared electrolysis barium hydroxide solution 
finally over all-glass ozonizer the was used the 
discharge tube and was excited h.t. using The 
current flowing through the discharge tube was measured galvanometer 
the plate circuit the detector (RCA connected the line across 
The wavetorm the discharge current was investigated 
feeding the potential drop across series resistance the line the 

vertical input the cathode-ray oscillograph. 
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given pressure, the potential applied the tube was increased 
gradually from zero, the galvanometer registered practically zero current till 
critical potential was reached which there was sudden rise the current 


h.t. 


RCA 
30, 
OSCILLOGRAP, 


electrical circuit diagram. 


flowing through the discharge tube. This was accompanied the initiation 
glow characteristic the gas the system and also the appearance 
high frequency pulses the sinusoidal current waveform. This critical 
potential called the threshold potential, Beyond this, the potential, 
applied the tube was increased regular steps and the current was 
measured both dark, and under irradiation 200-watt bulb, 


RESULTS 


Figure shows the variation and with for few typical pressures 
hydrogen. While these are discussed detail elsewhere (Jatar 1957) 
seen that the characteristics exhibit three distinct 
regions: (a) the initial region rapid rise near the inter- 
mediate sag region corresponding slow rise the conductivity 
followed (c) final fast rise conductivity. was pointed out (Jatar 
1957) that the three distinct regions the characteristics 
correspond the space-charge-free, the space-charge-limited, and the steady 
corona regimes. 

The suggestion space-charge-limited current above investigated 
from the standpoint the 3/2 power law, viz. 


the electrode distance kept constant. Here and refer the charge 
and mass the charged particle. This equation was first deduced Child 
(1911) for ions arcs and independently Langmuir (1913) for electron 
currents vacuum. While this obtained the supposition that there 
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Fic. Variation the discharge current dark, ip, and that under light, with the 
applied potential, Current arbitrary units, (r.m.s.). Numbers curves indicate 
pressure hydrogen Hg. 


ionization between the electrodes, Thomson (Thomson and Thomson 1934) 
has shown that this holds good when there uniform ionization. This equation 
has been further extended Langmuir and Compton (1931) for concentric 
cylinders, the current assuming negligible velocities emission, being 
given 


where the length the cylinder, the radius the collector (anode for 
electrons) and coefficient taken from the data given Langmuir 
and Compton (1931). This equation also applicable the case discharge 
coaxial cylinders when the collector surrounded space charge sheath 
radius first approximation, assumed that the change ro/r 
and hence due change negligible because the large value 
the inner cylinder, expected that the double logarithmic plots will 
linear. These are shown Figs. and which show that about 200 
Hg, the curves are linear beyond For pressures between 
and 200 these curves also indicate second linear region corresponding 
the region final fast rise conductivity beyond the transition potential 
these high pressures and potentials, not unlikely that the assump- 
tion not valid that does not change appreciably result change 


examination the slope these double log plots, however, indicates 
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Fic. plots for the relation. 


3.2 Log 


Fic. plots for higher pressures hydrogen. 


that the slope more nearly instead 3/2, expected from the 3/2 
power law. The data therefore further examined from the standpoint the 
Werner’s equation (1934), which reads 


where refers the applied potential, the starting potential, and 
are, respectively, the radii the outer and inner electrodes, the radius 
the zone intense ionization surrounding the gas pressure, the 
proportionality and equal This equation was derived 
Werner study Geiger counter action. Neglecting compared 
and compared the equation becomes 
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which was derived Townsend (1914) early 1914, the same reasoning 
but with simplified assumptions. According these equations the current 
While the general shape the curves the present case does correspond 
approximately parabola least above further examination revealed 
that the curves show appreciable deviation especially for higher pressures. 
This shown Figs. and which plotted against Ac- 
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Variation and with V(V— for the modified Townsend—Werner equation. 
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Fic. and with V(V— Vm) for higher pressures 
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cording Townsend’s equation, these curves should linear. The same 
holds good according Werner’s equation the assumption that does 
not change with increasing Figs. and show that varies linearly with 
only for the intermediate pressure range, from 100 Hg. 
interest note that for lower pressures, e.g. and Hg, while 
the curves for are linear, those for are not (Fig. 5). The curves for 
higher pressures (Fig. show three almost linear regions corresponding 
the three regions the characteristics. suggested earlier 
(Jatar 1957; corresponds the potential which space charges 
set in. change the slope the curves seen Figs. and therefore 
anticipated. Werner’s equation, given above, may modified 
where has the dimensions resistance and expresses the 
corona gap resistance. When becomes zero because the term 
brackets Werner’s equation becomes zero. Under these conditions, theo- 
retically space charges are active. This condition holds for the region 
above change slope seen Fig. lends support the above 
view. significant note, however, that Vg, does not become infinite, 
Fig. indicates decrease the corona resistance, the latter does not 
become zero any potential the range investigated here. This suggests that 
even these high potentials, while they not play prominent role, meager 
space charges are present. view the simplifying assumptions made 
deriving the equation and also applying the present case, the agree- 
ment more than satisfactory. concluded that the three distinct 
regions the characteristics correspond the space-charge-free, the 
space-charge-limited, and steady corona regimes. 
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THE VISIBLE EMISSION SPECTRUM 


VENKATESWARA RAO AND TIRUVENGANNA RAO 


ABSTRACT 


Two visible emission band systems attributed the ionized molecule 
were studied using the second order 21-ft concave grating spectrograph (dis- 
persion 1.25 Rotational analyses have been carried out for six bands 
system and six system III for which new vibrational analyses have recently 
been proposed Haranath and Rao. The rotational analyses the bands 
these two systems are consistent with the vibrational assignments. Rotational 
constants for the upper and lower states the two systems have been derived. 
Both the vibrational and rotational analyses systems and III indicate case 
(c) 1/2-1/2 and 3/2-3/2 transitions respectively for the two systems. Potential 
energy curves are drawn for the upper and the lower states each the two 
systems. The lower states dissociate into normal while the 
upper states dissociate into normal excited sug- 
gested that the lower state system III probably the normal state 
theoretical and experimental evidence. 


INTRODUCTION 


The emission spectrum chlorine was investigated early workers 
(Uchida and Ota 1928; Elliot and Cameron 1937 and 1938) and well known 
consist series bands the visible region degraded towards the red. 
These bands not correspond the well-known absorption bands chlorine 
the visible region and were attributed Elliot and Cameron (hereafter 
referred EC) the molecule from vibrational analysis and 
rotational analysis some the bands. proposed vibrational analysis 
the bands the basis two systems designated systems and 
and suggested that the observed features the two band systems indicate 
transition the type The following vibrational formulae were derived 
them. 


The rotational analysis some the bands belonging both the systems 
showed that the quantum number half integral, which confirmed that 
the two systems arise from the ionized molecule Although the vibrational 
analysis was supported the observed chlorine isotope effect, there were 
several inconsistencies, some which were pointed out recently Howell 
(1953). (i) The intensity distribution both the systems inconsistent 
with the values and derived from the analysis. (ii) The values the 
rotational constants and not run parallel the derived values 
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that the ratio not constant for each the pairs elec- 
tronic states systems and II. mentioned Howell for the neutral 
halogen molecules the ratio approximately constant for pair 
electronic states the same molecule. Using the rotational and vibrational 
data for found that this rule fails badly. (iii) The value 
the energy dissociation the normal state obtained from 
Sponer extrapolation 4.4 against 2.23 ev, which can indirectly 
determined from the accepted values the ionization potentials and 
and the dissociation energy (iv) detailed electron configurations 
and predicted states for indicate (Mulliken 1934) the probable transition 
expected for the two systems corresponds the transition 
electron from bonding orbital antibonding orbital and should 
accompanied considerable difference the relative values and 
(v) The large number gaps the vibrational schemes suggests funda- 
mental weakness the analysis which leaves large number prominent 
bands unclassified. 

view the above weak points EC’s analyses, Haranath and Rao 
(1956) re-investigated the emission spectrum and proposed vibra- 
tional analysis the basis three systems designated systems II, and 
III. The possibility another system (IV) was also indicated. detailed 
paper the analysis now the course publication the Indian Journal 
Physics. The vibrational constants systems and III are summarized 
below: 

System 20448 350.0 2.0 4.6, 
System III 347.5 2.0 5.5. 

can seen from Figs. and most the prominent bands belong 
these systems. The new vibrational assignments were supported the 
agreement between the observed and calculated chlorine isotopic separations. 
The proposed vibrational analysis these two systems entirely different 
from that EC. The vibrational constants seem indicate that these two 
systems may form two subsystems. 

The present investigations the visible emission bands were 
carried out 21-ft concave grating spectrograph with view checking 
the correctness the new vibrational assignments Haranath and Rao and 
derive the rotational constants the upper and lower states systems 
and 


EXPERIMENTAL 


The spectrum was excited high frequency discharge from 100-watt 
oscillator Mc/sec frequency. The discharge tube was made pyrex 
glass length and diameter. had two side limbs which 
served inlet and outlets for chlorine vapor. Chlorine was generated elec- 
trical heating cupric chloride separate tube and was allowed pass 
through absorbing tower containing before entering the discharge 
tube. The main tube was also provided with central tube length 
and was dipped liquid air. The pressure chlorine was regulated 


RAO AND RAO: EMISSION SPECTRUM 1559 


adjusting the current heating coil. was seen that chlorine vapor condensed 
both the central tube and the liquid-air trap which was introduced 
between the discharge tube and Cenco Hyvac pump. Occasional adjustment 
the liquid-air level the central flask was found necessary for obtaining 
pure spectrum. 

The spectrum was photographed the region 5000 4000 the second 
order 21-ft concave grating spectrograph which was recently set this 
laboratory one the authors (V.V.). The grating has dispersion 
1.25 A/mm and resolving power 1,00,000 the second order. Exposures 
varying from hours were found suitable for obtaining good spectra using 
Ilford Special Rapid plates. this region most the prominent bands 
belonging systems and III were clearly obtained. Plates were measured 
using iron arc standards and the vacuum wavenumbers the rotational lines 
were calculated. The relative accuracy the measurements the sharp 
lines 0.07 however, for few the broad lines the error measure- 
ment may great 0.15 


ROTATIONAL ANALYSIS 


the basis the predicted electronic states from electron configurations 
the possible electronic transition for systems and III expected 
the type The vibrational constants the upper and lower states 
these two systems indicate that they may the two components the 
above transition. Further the doublet splittings the two states are 
expected large that the transition may with both the 
states approximating Hund’s case (a). Since large multiplet splittings the 
order 610 and 690 are known the ground states and atoms, 
one may even expect tendency towards case (c). With case each 
band splits into two subbands and and each subband 
consists single and branch and very weak branch. Hund’s 
case (c) the two systems may better represented 3/2-3/2 and 
transitions with the over-all rotational structure remaining the same 
case (a). 

The wavenumbers and lines follow the equation (Herzberg 1950) 


where J+1 for the branch and for the branch. The effective 
rotational constants the upper and lower states were determined from the 
equation 


A:F(J) 


which the combination differences the upper and lower states are obtained 
from 


The detailed rotational analysis was carried out for six the bands belonging 
system and six the bands belonging system III. The vacuum wave- 
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numbers and rotational assignments for the bands 11-3 and 11-4 system 
and and system III are given Tables and respectively. 


35. 
36. 
37. 
38. 
39. 
40. 


System 
this system the 10-3, 11-8, 11-4, 11-2, and 12-3 bands 21669, 


VACUUM WAVENUMBERS AND ROTATIONAL ASSIGNMENTS FOR SYSTEM I 


P(J) 


11-3 


TABLE 


R(J) 


975.23 
974.67 


P(J) 


21347 


346. 
345. 
342. 


341 


339. 
337. 
335. 
333. 
330. 
328. 
325. 
322. 
319. 
316. 
312. 


305. 
301. 
297. 
21292. 


278 
273 
268 
262 
256 
251 
244 
238 
232 
225 
219 
212 


11-4 


R(J) 


21347. 
346. 
345. 
342. 
339. 
337. 
335. 
328. 
322. 
319. 
315. 


21304. 
300. 
296. 
292. 
282. 
277. 
272. 
267. 
261. 
255. 
250. 
243 
231. 
224. 
217. 


21041, 21976, 21349, 22610, 22280 respectively were found suitable for 
carrying out the rotational analysis. All these are accompanied weaker 
isotopic heads the longer wavelength side belonging species. 
However, for the 10-3, 11-3, and 11—4 bands the rotational lines the 
more abundant isotope can clearly distinguished from those the other. 
Hence these four bands are best suited for accurate determination upper 
and lower state rotational constants. For the and 12-3 bands the rota- 
tional lines the two species are partially superposed. 


21974.67 
972.85 
971.55 
970.36 
971.55 
970.36 
950.14 
928.21 
914.91 931.53 
918.57 
894.35 913.88 
864.72 887 .65 
881.76 
851.41 875.79 
844.73 
863 
842.85 
210.80 
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TABLE 
VACUUM WAVENUMBERS AND ROTATIONAL ASSIGNMENTS FOR SYSTEM III 
12-1 12-2 

P(J) R(J) P(J) R(J) 
5.5 
6.5 23616.38 982.71 22987 
7.5 614.66 
8.5 612.59 985.14 
9.5 610.38 23617 976.86 983 
10.5 608 615.55 974.60 
605.41 613.62 972.06 980.17 
12.5 611.53 969 978.09 
16.5 601.24 968.11 
17.5 952.84 965.12 
18.5 581.96 594.94 949 .04 961.91 
19.5 591.48 944.87 958.55 
20.5 940.78 954.95 
21.5 936.18 951.08 
564.39 580.00 931.60 947 
23.5 575.70 926.76 
24.5 554.40 571.24 921.79 938 
26.5 61.84 911.10 929 .39 
28.5 532.11 551.69 899.76 919.43 
29.5 546.40 893.71 914.06 
30.5 519.68 540.72 887 .55 
31.5 513.29 534.89 881.21 
32.5 506.61 528.90 874.76 897 .05 
34.5 
35.5 509.79 
37.5 


The numbering fixed for 10-3, 11-3, and 12-3 bands comparing the 
ground state combination differences. These combination differences are 
found agree well within the accuracy the measurements, thus confirming 
their vibrational assignments. Similarly for the and 11-4 bands the 
assignments are made. With numbering derived above for the 10-3, 
and 11-3, 11-4 bands found that the combination differences the upper 
states each these two pairs agree well. The combination differences for 
upper state turn agree with those Thus the rotational analysis 
these six bands consistent with the vibrational assignments. in- 
teresting note that although four the bands 10-3, and 11-4 
belong two different systems EC’s vibrational analysis the absolute 
numbering for these bands the same both the cases. The rotational lines 
some the bands are extended high values our present work that 
has become possible determine the values. The rotational analysis 
also shows that half integral thus confirming EC’s earlier conclusion that 
the emitter these bands 


Rotational Constants 
After fixing the numbering the lines the manner described, the 
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constants and are determined plotting against 
(J+1/2)? both for the upper and lower states. The resulting and values 
are shown Table deriving and the average lower state com- 
bination differences the 11-3, and 12-3 bands were these 
have well resolved and branches. Figure reproduction the 11-3 
band. For and the band taken the band has superposed 
and branches throughout. The band reproduced Fig. 


TABLE III 
THE BAND ORIGINS, AND VALUES FOR SYSTEM 


*Determined from the averages of combination differences. 


The constants and are derived from the 10-3 band and and 
from the average upper state combination differences 11-3 and 11-4 
bands. only few the rotational lines are identified and 12-3 
bands because superposition, rotational constants and are not 
derived. The band origins shown Table III are determined plotting 
against according the equation (Herzberg 1950) 


examination the 11-3 and 11-4 bands shown Figs. and shows 
the rotational lines the branch down 1.5. Although the and 
lines are superposed low values, important conclusion can drawn 
regarding the nature the transition system both these bands. the 
line clearly observed. the short wavelength side this 
line the sudden drop intensity the rotational lines may taken 
indication that the first line has 1.5 and 0.5. Thus the upper 
case (c). 

Since the and are accurately determined the difference these 
two taken approximate value for and was then deduced from the 
obtained from and The and thus derived are given Table IV. 
The and values for the upper and lower states are determined the 
usual manner (Herzberg 1950) and are shown Table IV. 


Fic. Enlargement the 11-3 band system showing the assignments the 


rotational lines. 
Fic. Enlargement the band system showing the assignments the 


rotational lines. 
Fic. the band system III showing the assignments the 


rotational lines. 
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TABLE 
CONSTANTS FOR SYSTEM 


Lower state Upper state 
0.2716 0.193; 
0.0017 0.0017 
2.233 


System 

this system the bands 12-1, 11-2, and 12-2 23327, 
23624, 23925, 21476, 22694, and 22990 respectively, are found suitable 
for rotational analysis. the first three bands which have clearly separate 
and branches the numbering fixed comparing the lower state 
combination differences. These are found agree exactly within the accuracy 
the present measurements, thus confirming that all three have common 
lower level. Similarly the numbering was fixed for and bands 
and again good agreement was found. 

Using the numbering derived above found that the upper state 
combination differences for the pairs bands 11-2, 12-1, and 12-2 
agree with each other. the six bands analyzed here the and bands 
belong system and 12-1, 12-2 belong system whereas the and 
7-2 bands are unclassified EC’s analysis. 

Rotational Constants 

The rotational constants for the upper and lower states this sytem are 
summarized Tables and VI. the band has superposed and 


TABLE 
THE BAND ORIGINS, AND VALUES FOR SYSTEM III 


7-2 21475.40 


*Determined from the averages of combination differences, D’ is negligibly small. 
branches throughout the and were determined from the and 
12-2 bands only. The constants and By; can fitted into the equation 


Fic. the band system III showing the assignments the 


rotational lines. 
Fics. and Fuess spectrograms showing the prominent bands systems and 
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showing that linear The values shown Table calls for 
comment. almost equal the This irregularity may probably 
done because perturbations the various excited electronic states which 
cross the upper potential curve. 


TABLE 
CONSTANTS FOR SYSTEM III 


Lower state Upper state 


266 
0.0017 0.0016 
1.891 2.266 


The 7-2 band reproduced Fig. shows line the long 
wavelength side vo. There clear gap between this line and the band 
origin. This may taken indication that the first line has 2.5 
and 1.5. Thus the upper state has 3/2 and the transition 


POTENTIAL ENERGY CURVES AND DISSOCIATION PRODUCTS 


Using the constants B,, D,, and for the upper and lower states derived 
from vibrational and rotational analyses systems and III the potential 
energy curves are drawn employing the Morse function. These curves are 
shown separately for systems and III Figs. and For deriving the low 
lying electronic states from the states the separate atoms and 
may limit ourselves the consideration the three combinations 
and These combinations 
will give rise the following states according the cor- 
relation rules 


the basis electron configurations discussed later, the lower states 
systems and III are regarded and arising from the con- 
reasonable assume that the state, the lower state system III 
dissociates into while the state, the lower state 

The lower state systems and III can correlated unambiguously 
with one the states arising from The upper states 
these systems may identified and configuration 
For the dissociation products the upper states systems and III 


Qe 
le 
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only two possible combinations need considered, namely 
and The latter possibility can easily ruled out each 
case will lead values dissociation energies the lower states which 
are smaller than some the observed vibrational levels these states. The 
excited and states systems and III dissociate therefore into 
This clearly illustrated Figs. and for systems 
and III respectively. Experimentally will not possible reach 
unambiguous decision whether the atom produced the 
state. However, there some evidence from consideration the 


20 2-5 3:0 3.5 


Fic. Potential energy curve, system 


Fic. Potential energy curve, system III. 
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bands and shown Table VII that the interval the 
upper state shows decrease with increasing vibrational quantum number. 
this interval decreases zero the dissociation limit approached then 
the two components dissociate into Theoretical evidence 
outlined Mulliken (1928a, 1930a, and series papers and 
discussed Durie and Ramsay (1958) for analogous upper state ClO 
shows that the halogen atom probably produced the state. 


TABLE VII 


Classification, System System Ap, 
10-2 22390 22301 
11-2 22694 22610 
12-2 22990 22914 
13-0 24562 24485 


From the dissociation products indicated above for system III possible 
derive the dissociation energy the lower state which probably 
the normal state estimate the dissociation energy this state 
can made one the two ways. (i) means extra- 
polation the observed vibrational levels the upper state. Applying the 
equation for the upper state obtain value 15094 
assume that this value 20% high the true value upper state 
obtained 12075 now add this value the value system 
III and subtract excitation energy the atom (from Moore’s 
tables, 1949) obtain D’? 20992 2.60 ev. (ii) The value 
can also determined directly from Birge-Sponer extrapolation the 
observed vibrational levels the lower state. This obtained 2.42 ev. 
assume that this value 10% high, seen agree closely with the in- 
directly determined value (2.23 ev) the energy dissociation the 
normal state 


ELECTRONIC TRANSITIONS 


Since isoelectronic with reasonable expect the same 
electron configuration and term type for the ground state. Following 
Mulliken’s notation (1932) this may written 


> 2 2 4 3 9 
The first excited state may also derived from 
2 2 2 3 4 9 
The transition from the excited state the ground state the 
type giving rise and components representing systems 
and III respectively. The transition electron from the bonding 
orbital the corresponding antibonding orbital should accompanied 


considerable reduction the vibrational frequency. This consistent 
with the observed fact that value decreases from 655 348 similar 
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decrease the vibrational frequency also observed the case similar 
transition ClO (viz. 868 557 The increase vibrational frequency 
from 565 655 the normal states and also consistent with 
the fact that the electron removed the process ionization molecule 
belongs the antibonding vz, orbital. 

The coupling constant ‘A’ for each the above states expected 
intermediate between the atomic coupling coefficients ‘a’ and atoms 
i.e., between 588 and 645 and may taken ~610 
interesting note that the pairs bands systems and III shown Table 
VII are separated wavenumber interval varying between and 76. 
This small separation indicates that the interval both the states 
the same order magnitude predicted theoretically. the interval 
the lower state greater than the upper state then the 
system expected the short wavelength side the component 
system observed. However, for the two component systems CIO the 
situation slightly different the component occurs the 
short wavelength side the component. This implies that 
the interval the lower state smaller than that the upper state. 

The values for upper and lower states the two components systems 
obtained from the present analysis are somewhat different from each other. 
Hund’s case (a) strictly holds both for the upper and lower states they 
should equal. If, however, there departure towards case because 
spin uncoupling, then the derived values should regarded 
values where 


where the true value and the sign refers the component. The 
values differ the amount For the lower state this difference 
two values. The observed values for the lower state differ 
much 0.0021. This discrepancy can explained only assume that 
the coupling the lower state shows tendency towards case (c) rather than 
case (6). similar discrepancy for the upper state also points the con- 
clusion that the coupling this state also shows tendency towards case 
(c). not unreasonable expect case (c) coupling both the states when 
doublet intervals the order 610 are expected. Hence systems and 
III may better represented case (c) and 3/2-3/2 transitions 
respectively. 


conclusion the authors would like express their deep indebtedness 
Dr. Herzberg for his interest the work and for the contribution some 
castings and blue print his mount Prof. Rao, which were used 
the setup the 21-ft concave grating spectrograph. One the authors 
(V.V.) grateful Krishnaswamy Chetty for his valuable guidance 
and advice setting the grating. The authors’ thanks are due Prof. 
Mulliken and Prof. Rao for their interest this work. 
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POLARIZATION FLUCTUATIONS FERROELECTRIC 


ABSTRACT 


The theory thermal fluctuations electrical polarization ferroelectric 
crystal considered with special attention temperatures the neighborhood 
the Curie point. Models giving second-order transition are discussed thermo- 
dynamically, and the free energy with its derivatives yield the mean value and 
variance dipole moment. The fluctuations and their spectral density are 
related respectively the real and imaginary parts the susceptibility the 
crystal conformity with theorem, which not violated the 
presence constant applied electric field. The theorem shown 
applicable equilibrium and non-equilibrium conditions, and its extension 
lambda-point phenomena generally discussed. 


INTRODUCTION 


ideal gas the variance the number particles small sub- 
volume large enclosure equal the average number particles. This 
result, which essentially consequence probabilistic considerations, re- 
quires modification real gas take into account the molecular inter- 
actions. this done thermodynamically (Fowler 1936) then high 
approximation the variance may expressed terms the isothermal 
compressibility the gas which appropriately represents these interactions 
for the purposes fluctuation analysis: 


Under normal conditions this variance still the order and the 
distribution gaussian. the Appendix the variance evaluated for 
various forms the equation state. However, near critical point (e.g. 
the gas-liquid critical point) the compressibility becomes very large and 
fact the first and second derivatives and vanish simultaneously 
this point and determines the fluctuations. these conditions the 
fluctuations are greatly enhanced and the effect visible opales- 
when light scattered from the fluid. 

The formula given above longer valid the critical point where the 
third derivative determines the fluctuations which, instead 
being order are now order Furthermore the distribution law for 
longer gaussian but the sharper form 

relevant consider what form critical fluctuations should ob- 
servable solids. The lambda-point transitions which occur ferroelectric 
crystals are particularly interesting since the macroscopic phenomenon 
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readily observed and the fluctuations may detected electrically spon- 
taneous electromotive forces electrodes attached the crystal. Various 
classes crystal exhibit spontaneous polarization over certain range 
temperature whose upper limit there occurs transition with very large 
peak permittivity and associated specific heat anomaly. Since 
least some cases these transitions are second order (and also because lambda- 
point phenomena are general interest physics) shall discuss the 
fluctuations which occur the vicinity second order transitions. 

specific model dipole polarization will considered indicate the 
statistical background the problem and general thermodynamical dis- 
cussion will given which independent any specific model. Compliance 
the fluctuations thermal equilibrium with Nyquist’s theorem established 
implying, expected, that enhancement the dipole moment fluctuations 
around the Curie point must correspond peak the dissipative part 
the permittivity the same temperature range. The application steady 
electric field modifies the average dipole moment and its fluctuations. 
independent derivation readily obtained from the author’s 
(Burgess 1955, 1956) identifying the process with dipole 
aligning itself parallel the field and the process with anti- 
parallel alignment; the kinetics the change alignment lead the variance, 
relaxation time, and spectral density the net dipole moment. 

Although the relation between the macroscopic and fluctuation aspects 
second order transition discussed the specific terms ferroelectric 
crystal, can extended other cases where such phase transitions occur, 
e.g. alloy (such CuAu) with long-range order, liquid helium 
the lambda point, superconductor the absence magnetic field. 
all these processes co-operative effects are strongly acting produce 
ordering below the critical temperature, and the degree order decreases 
rapidly this temperature approached. 


DIPOLE MODEL 


Consider crystal volume containing identical dipoles all situated 
equivalent sites with moment which may oriented parallel 
antiparallel with the applied uniform external field any instant, 
the dipoles are the parallel position can introduce order parameter 
net dipole polarization. 


(2) 
where N/V. The net dipole moment Nmx. 

real crystal, other sources polarization may present but these 
will ignored order clarify the situation; this will done justifiably 
since near the transition the dipole polarization will dominate over other forms 
the present model. 

The local electric field acting each dipole the Weiss form 


(3) 
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where the Lorentz factor, however, does not need given the 
value but simply expresses the linear dependence the local field 
the polarization. 

The internal energy 


P 
0 


which represents the work done orienting the dipoles the presence the 


local field. 


The entropy the configuration 


and the free energy 
(4) F=U-TS 


where Stirling’s theorem has been applied the factorials. The condition 
equilibrium obtained minimizing the free energy with respect the 
dipole moment. 


giving 


This equilibrium condition equivalent the Boltzmann ratio for the 
populations the parallel and antiparallel states: 


readily deduced that the critical temperature, which the upper 
bound for non-zero root this equation exist when given 


(8) 


Just below the Curie point 


1/2 


Im 
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The susceptibility appropriately defined dP/dE rather than P/E 
since this removes the effect the spontaneous polarization below the Curie 
point and furthermore the differential coefficient required fluctuation 
problems relate the small variations and around the average values. 

Now the reciprocal the susceptibility 


which positive above and below the Curie point. the region there 
are two forms for the zero-field values: 


This model gives polarization and entropy which are continuous through 
the Curie temperature although their derivatives are discontinuous, charac- 
teristic second order transition. 

discussing the thermodynamical aspects the transition and the thermal 
fluctuations polarization necessary have the first four derivatives 
free energy with respect dipole moment the equilibrium value and for 


(14a) 


the absence applied field the transition occurs with 
which implies that the critical point 


FLUCTUATIONS 


Total Fluctuations 
The total mean square fluctuation the dipole moment can determined 


from the free energy function. expand about the equilibrium value 


where the suffix zero denotes the mean value equilibrium. see that 
long does not vanish have for small fluctuations 


(16) 


E 
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Now 


the form the probability distribution function which implies that 
has gaussian distribution about with variance 


This fluctuation polarization will manifest itself fluctuation 
potential difference between electrodes applied the crystal. the next 
section these observable fluctuations will related the dipole moment 
fluctuations. 

The expressions above for fluctuations are valid whether the temperature 
the crystal above below the Curie temperature (but not too near it) 
and also the presence steady electric field The potential difference 
fluctuation corresponds the classical equipartition formulation applied 
the electrostatic energy the crystal, might expected from equation (7). 

generalization the preceding arguments readily achieved such that 
quite general dependences the internal energy and the entropy the 
order parameter will still lead the equipartition result for the total fluctua- 
tions. 

Thus, writing the free energy its most general functional dependence 
the applied electric field and the dipole moment 


U(E, 


gives, for the equilibrium condition Mo, 
U'(E, Mo) = kTS'(M_), 


where prime denotes partial differentiation with respect 
The susceptibility found from the variation with infinitesimal 
change while equilibrium maintained: 


or 


Thus var demonstrating the general validity 
equation (18) under non-critical conditions where does not vanish. 
Just higher derivatives need taken into account considering 
critical point fluctuations gases, the ferroelectric need consider 
additional terms the development the free energy electric field 
terms first approximation, near but not the Curie point, 
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have from (146), (12), and (13) for the dipole moment fluctuations the 
absence applied electric field: 


These expressions indicate that the dipole fluctuations have pronounced 
maximum (theoretically infinite the first approximation) when 
This infinity removed going the fourth derivative, the first 


three all vanish; since tind that 


This not gaussian distribution and gives dipole moment variance 
(22) 


which implies that crystal having the fluctuations dipole 
moment given (19) can used temperatures quite close the Curie 
point. 

‘The mean internal energy due fluctuations the total dipole moment 
about zero above the Curie temperature 


terms number fluctuations, the variance the number dipoles 
oriented parallel the field found from 


(24) var var M/4m? 0.29 


which illustrates the tendency critical fluctuations order rather 
than order 

The mean fourth power fluctuations and are also interest for the 
type distribution given equation (21): 


which gives for the increment free energy due fluctuations around the 
critical point: 
(26) kT,/4 


contrast with the value when normal quadratic fluctuations occur. 

Some reservations about the results obtained the Curie point should 
mentioned. The representation the long-range co-operative action the 
dipoles parameter requires reexamination the Curie point since 
localized interactions will tend predominate. This evident the macro- 
scopic properties ferroelectrics and far there seems theory which 


{| 

| 
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predicts finite value for the vermittivity the Curie point, and furthermore 
the occurrence domains alternating polarities represents complexity 
the interaction. 

When the crystal not near the critical temperature the proportionality 
the variance the dipole with the volume (equation 18) 
implies that this formula vatid for each subvolume the crystal and that 
the fluctuations all elements are independent and additive. Near the Curie 
point covariance becomes important and the model used inadequate because 
the need include distance-dependent correlation between the fluctua- 
tions any two elements. Landau has suggested (Frenkel 1946) that this 
may represented including term |grad the energy density 
associated with the polarization. Since there satisfactory model for this 
case our remarks will remain qualitative; will noted likewise that the 
finite maximum the susceptibility the Curie temperature not predicted 
existing theory. 


3.2 Spectral Density Fluctuations 

Since thermal equilibrium the fluctuations obey Nyquist’s theorem 
may readily deduce their spectral density terms the dissipative part 
the susceptibility. frequency the admittance parallel plane crystal 


where X(f) the complex susceptibility Hence the spectral 
density the fluctuation potential difference 


and the spectral density the short-circuit current fluctuations 
(28) Si(f) 
and thus the fluctuations dipole moment are 


When integrated over all frequencies this correctly leads the earlier ex- 
pression (18) for the total moment fluctuations: 


using the relation and noting that the susceptibility 
defined earlier dP/dE synonymous with the low frequency limit the 
real part X’(0) the complex susceptibility. 

These results indicate that measurements fluctuations polarization 
ferroelectric crystal give more information about the internal co-operative 
and relaxation mechanisms than would measurements the complex sus- 


| 
7 
7 
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ceptibility the crystal. The peak near the Curie temperature will 
give rise proportional peak the fluctuation current measured 
narrow band frequency the dissipative mechanism can characterized 
polarization relaxation time then under conditions where the order 
parameter small compared with unity (so that linear function 
and can regarded the driving force) the complex susceptibility 


(30) 


Thus the peak the temperature taken through should most 
pronounced low frequencies where 
For this model the spectral density the dipole moment fluctuations 


The fluctuations observed electrically will depend whether the measurement 
circuit which the crystal connected has lower higher impedance than 
the crystal. The variation with temperature will also determined the 
variation which related the tangent the loss angle the dielectric: 


c 


where the spectral relaxation time appropriate the representing the 

The equivalent circuit which represents the admittance the ferroelectric 
crystal and its fluctuation properties consists the vacuum capacitance 
parallel with series combination the dipole capacitance 
AX(0)/L and the loss resistance The spectral density 
the potential difference fluctuations the terminals the system 


4kTR 


(33) 
that the major change S,(f) passes through will due the 
change with temperature. 


The total fluctuations potential difference integrated over all frequencies 
can related the dipole moment fluctuations: 


where the zero frequency permittivity the crystal. 
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3.3 Application the Theorem Polarization Fluctuations 

the most general sense proceeding such that the state the system can 
specified number variable (e.g. number particles state then 
g(y) the rate which the reaction tends increase while r(y) the rate 
which tends decrease such that dy/dt g—r, the steady state 
given 


and the mean square fluctuation about the mean value given 


where the prime denotes differentiation with respect 
The relaxation time for the reaction, defined terms the return the 
steady state after slight perturbation, given by: 


corresponding spectral density the fluctuations given 
(37) 


These results have been applied generation and recombination processes 
involving mobile carriers and ions semiconductors. However, they are 
much wider applicability and can used discuss diffusional processes, 
photon fluctuations, and phenomena. will now apply 
the dipole polarization problem identifying the variable with the number 
dipoles oriented parallel the field, and consider the 


antiparallel dipole parallel dipole +2mE’. 


would then write and the forms 


where constant and the height the potential barrier between the 
two positions. 

evident from earlier work the representation that only the 
ratio g/r important determining the moments the variable under 
consideration. Thus the parameters and are not relevant for the moments 
although they are important determining the rate which the 
process proceeds and hence the relaxation time and spectral density. 

seen that and depend directly and also indirectly through 
which linearly related N;. This feature the consequence the 
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co-operative manner which the state the crystal whole determines 
the transition probabilities any one dipole between the parallel and anti- 
parallel states. The logarithm related the change free energy 
when one dipole reversed from the the parallel position: 


and this quantity zero g). 
The derivative this expression 


and this must positive equilibrium state corresponding minimum 
the free energy. This seen the case since 


and E/x taken positive corresponding the crystal polarization having 
the same polarity the applied electric field. 

The spectral time constant given equation (36) now 


exp 


To = 172° 


Thus comparison with (32) seen that the introduced there (for 


which agreement with the earlier equation (18) for the dipole moment 
fluctuations. 

For comparison, had been naively assumed that the dipole statistics 
obeyed the binomial distribution corresponding Bernoulli ‘‘trials’’ with 
seen that approaches this value for but that when 
only few times the thermal energies are not sufficient override 
the co-operative effect (represented which enhances the fluctuations 

The theorem based statistical considerations and does not require 
thermodynamical equilibrium obtain and consequently can used 
determine the validity free energy analysis when uncertain. 
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example the verification (18) shows that whether not steady electric 
field being applied the crystal the variance correctly given this 
formula. 

When the critical condition approached, equation (35) longer correctly 
gives the variance and necessary higher derivatives and 
just for the free energy. fact may shown quite generally that 
the free energy related and 


whence 

(41) 
which vanishes the steady state where go. 

Thus 


kT 2 72 ” 


seen that must also hold because the necessity for dy/dt 
odd function requiring the even derivatives (r—g) vanish. 
Thus the critical point corresponds contact the second order the two 
curves g(y) and r(y) occurring common point inflection where the 
slope both curves zero. This demonstrates that the point region 
relatively large fluctuations because the restoring force when departs 
slightly from small since the rate the reaction depends the cube 
the disturbance: 


(45) dy/dt 


6 — (go yt... 


From stability considerations clear then must positive- 
also seen that because the cubic behavior impossible accord 
relaxation time the reaction the critical point. this condition the 


exp 


which the analogue the distribution terms free energy critical 
point (cf. equation 21). 


“= 
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These considerations relating the higher derivatives (r—g) are directly 
related the derivatives the free energy and correspond equations (14) 
critical point. 
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APPENDIX 


DENsITY FLUCTUATIONS IMPERFECT GAS JUST ABOVE THE CRITICAL 
imperfect gas departs from the ideal gas (whose equation state 
NkT) virtue intermolecular forces and the finite size the 
molecules. Such gas has critical point Ve, where 


little above the critical temperature can take 


The number molecules volume will have gaussian distribution 
with variance. 

var 

and this will always the form: const. where the constant 
order unity and determined the precise character the 
the gas. 
(i) For the van der Waals equation state 


NRT, 


Pe ’ c 3d, N kT, 27b ’ 


(ii) modification the van der Waals equation state consists 
replacing a’/T with the following result: 


both the van der Waals and the Berthelet equations the critical point 
parameter 8/3, which considerably below the values observed 
for real gases. Better agreement given by: 
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(iii) equation 


for which 


26, Pe ’ N ki 4b 
giving 
3.695 
and 


The variation the numerical constant with the equation state demon- 
strates the sensitivity the fluctuations the form the intermolecular 
forces, expressed through the compressibility. This analogous the suscepti- 
bility which, the ferroelectric crystal, represents the effect any one dipole 
the polarization due the other dipoles, and also directly enters into the 
magnitude the fluctuations the net dipole moment. 


NOTES 


MILLIMICROSECOND PULSES FOR CHECKING 
FAST TRANSIENT ANALYZERS! 


DENHOLM 


INTRODUCTION 


During study fast pulses generated corona discharge, generator 
was needed which could produce pulse with rise time shorter than those 
the discharge and amplitude sufficient give satisfactory deflection 
high writing speed oscilloscope. attenuating this pulse the level 
the comparatively small corona pulses being measured, and then amplifying 
through the pulse measuring system, indication the minimum rise 
time the system could obtained. 

The corona pulses which were being studied are the prime cause radio 
interference from high voltage transmission lines, but not television inter- 
ference, which generally traced bad contacts, i.e. the sparking small 
gaps. Apparently, the time for such gaps change from nonconducting 
conducting very short, which confirmed experiments this laboratory 
the generation step functions for checking the transient response high 
voltage dividers (Creed 1958). The spark gap has been used produce very 
fast pulses the higher voltages (e.g. kv, Fletcher 1949) and there 
evidence that under certain modes operation mercury relays produce fast 
pulses due spark discharge the contacts close (Lewis and Wells 1954). 
The above suggested that very small gap could used generate fast 
pulses using comparatively low voltage which would need only conventional 

components. 


CIRCUITRY 


The circuit shown Fig. When the hydrogen thyratron tripped, 
capacitor discharges and fast rising pulse developed across the 100-ohm 
resistor. Coaxial cable (D) charged through the 2-kilohm resistor the 
voltage which the small gap breaks down, whereupon the cable discharges 
through the 50-ohm resistor giving rectangular pulse with very fast rise 
and fall. The pulse duration equal twice the electrical length the 
cable. After the cable has discharged, the gap extinguishes, the cable re- 
charges, and again the gap breaks down produce another pulse. Pulses 
recur, shown Fig. until the cable fails charge voltage sufficient 
break down the gap. 

The relays the circuit allow the generation positive negative pulses 
without resorting components such transmission line transformers. 
Negative impulses are the more usual output the impulse generator circuit 


N.R.C. No. 4909. 
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Fig. and are given point when grounded. Grounding and 
taking the output from obviously gives positive pulse, but the presence 
the heater winding usually unfortunate, only because introduces 
capacitance ground which limits the rate rise voltage across AB. The 
circuit Fig. removes such effects using the high thermal time constant 
the hydrogen thyratron heaters. When positive pulse required, relay 
opens the heater circuit just before the thyratron fired and then closes 
afterwards. Relay operates reversing switch give the desired output 
polarity. 

Key controls the operation relays and the key 
pushed the left (positive pulse), contact broken, relay de- 
energized, and the heater tube isolated from its transformer. the 
end the key travel, contact made which initiates the tripping pulse 
for the hydrogen thyratron causing contact close and remove 
the bias tube (Contact opens ensure that only one tripping 
pulse occurs.) When the key pushed the right (negative pulse) 
de-energized reverse the output connections before contact made 
fire 


PERFORMANCE 


The rectangular pulse generated the circuit shown Fig. The 
rise time, measured the oscillogram, seconds, but this 
close the estimated time constant the CRO plate system, the actual 
rise time the pulse may less than this. The fall time the pulse about 
seconds. The steps the top the pulse seem due the 
oscilloscope rather than the pulse generator. 

The average time between the time base tripping and the first output pulse 
seconds, but the actual values are very inconsistent, apparently 
because the performance the gap, since the oscilloscope time base normally 
trips stably. the gap with ultraviolet light gave improvement. 
However, the time between the first pulse and the second was reasonably 
consistent (standard deviation, 5%) was the amplitude the second 
(standard deviation, 5.6% cf. first pulse, standard deviation, 9%), that 
relatively stable display the second pulse may obtained the first 
pulse used for time base triggering. seems that the initial pulse conditions 
the gap some way, perhaps residual ionization, that the state the 
gap prior each second pulse essentially the same. 

hope continue work this field. 


CONCLUSIONS 


simple circuit for generating millimicrosecond pulses has been given 
which uses only conventional components and two small spheres. expected 
that almost any electrode geometry would satisfactory, long the 
gap was small. 

Apart from its original purpose, which was for checking amplifiers and 


A 


rectangular pulses generated spark gap. Pulse duration, 
seconds. Pulse amplitudes, No. 1—145 volts, No. 2—80 volts, No. 3—80 volts, 
volts. 

Fic. Output pulse. Duration, seconds. Amplitude, 125 volts. 
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circuitry, this generator circuit could useful for the fast operation 
electro-optical shutters such the image converter. 


The careful assistance Lackey, who built the equipment and made 
the measurements the generator performance, gratefully acknowledged. 


FLETCHER, 1949. Mass. Inst. Technol. Research Laboratory for Insulation Research, 


Tech. Rept. No. 20. 
Lewis, and 1954. Millimicrosecond pulse techniques, 102. The 
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NEW TRANSITION MOLECULAR 


CARROLL? 


the course studying the system new headless 
band with origin 1589.745 was discovered. The source employed the 
experiments was hollow cathode discharge through pure nitrogen low 
pressure. The excitation conditions were vigorous and the spectrum the 
visible and ultraviolet regions was dominated the first negative system 
and the bands The new band was photographed the 
fourth order the three-meter vacuum spectrograph the National Research 
Council Canada. Details the instrument and the procedure employed 
are given recent publications Brix and Herzberg (1954) and Carroll 
(1956). 

reproduction the new band shown Fig. The zero gap can 
located with certainty and the rotational analysis follows straightforward 
manner. The data for the new band are presented Table and the usual 
graphical methods the analysis leads the following constants: 


62,903.18 


B’ 


first sight the band appears 0-0 transition, but the values quoted 
above show that the level none the many states and for 
which values are known involved. This raises the problem identifying 


the transition. 


N.R.C. No. 4930. 
Research Council Postdoctorate Fellow, 1953-55. Present address: Physics 


Department, University College, Dublin, Ireland. 


Can. J. Phys. Vol. 36 (1958) 


= 


CANADIAN JOURNAL PHYSICS. VOL. 36, 1958 


TABLE 

WAVE NUMBERS THE LINES THE NEW BAND 
P(N) R(N) P(N) 
14.29 13.45 80.54 
87.81 21.43 04.54 
84.04 25.10 00.18 89.40 
76.20 32.14 91.15 98.21 
72.24 35.59 86.55 .65 
60.07 45.74 72.57 05.79 
51.95 52.43 63.20 10.81 
35.07 65.34 44.51 
30.83 68.44 
26.48 71.51 


The clear alternation intensities leaves doubt that the emitter 
nitrogen. The absence any structure the zero gap requires that 
must zero for both states. Furthermore, triplet states are not involved 
sign spin splitting detected, even low values. therefore seems 
that the band must arise from transition. The fact that, 
the intensity alternation, the lines odd values are the weaker ones 
shows that the lower state either state. 

very unlikely that the band due transition between excited 
states this were the lower state could scarcely lie above 80,000 
otherwise the upper state would lie improbably high energy. The 
lower level would also have metastable has not been observed 
transitions the ground state. consideration molecular orbitals there 
seems only one state likely satisfy these requirements, namely the 
state arising from the configuration 
(cf. Herzberg 1946). But this state almost certainly the lower state the 
fifth positive bands (cf. Gaydon 1944 and Lofthus 1956), which definitely 
not the same the lower state the new transition. therefore seems that 
the band cannot assigned This consistent with the experimental 
conditions which, mentioned above, were such that ionized spectra were 
strongly favored. 

attempting assign the band will consider the following possibilities: 

transition between excited states identify the lower state 
consider the electron configuration and again confine our attention 
states which are not likely have excessively high energy. The configura- 
and two states the only one likely interest. One the 
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states undoubtedly the upper state the system; one the 
states could the lower level the new transition. this the lower 
state the new band should lie about 65,000 Furthermore 
state would metastable with respect all other known states 
which would explain the fact that hitherto has escaped observation. 

pointed out the author that the value from the above analysis (1.8801 
agrees very well with the value (1.879 for the ground state 
might therefore tentatively identify the band the 0-2 transition 
from state 67,220 the ground state, this interpretation 
one would expect such bands 0-0, 0-1, and 0-3 present 
search for these bands, however, proved negative. further argument 
against this interpretation arises from the fact that that one would 
expect the strongest bands the system lie the sequence. just 
possible, however, that the band the transition the ground state 
This would satisfy the requirement that the strongest bands lie the 
0-0 sequence and conceivable that the strongest these might 2-2. 
Failure observe other bands the sequence might due relative 
weakness and overlapping other structure. 

(ground state) transition This transition has been 
observed for the isoelectronic molecules and The Mulliken bands 
form headless system with 43,227.23 while the bands 
Douglas and Routly (1954) form similar system with origin 45,844.65 
tempting identify the present transition the analogue these 
two systems. This interpretation, however, requires that have not only 
remarkably stable ground state, but also almost equally stable excited 
state. Because Coulomb repulsion doubly ionized molecule one would 
expect the electronic states considerably less stable general than those 
the singly ionized neutral molecule. 

can seen from the present discussion that the assignment the band 
cannot made with any certainty. Other transitions the states involved 
must observed before final decision can taken. 


conclusion the author would like express his indebtedness Dr. 
Herzberg and Dr. Douglas for valuable discussions the above problem 
and Mr. Shoosmith for his expert assistance operating the vacuum 
spectrograph. 


Brix, and HERZBERG, 1954. Can. Phys. 32, 110. 
1956. Can. Phys. 34, 83. 

1944. Proc. Roy. Soc. 182, 286. 

HERZBERG, 1946. Phys. Rev. 69, 362. 

1956. Chem. Phys. 25, 494. 
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ERRATA AND SOME ADDITIONAL COMMENTS CRITICAL STUDY 


VIBRONIC INTERACTION CALCULATIONS* 


ANDREW LIEHR 


TABLE 


Intensities calc. basis of: 


Eq. (2) Eq. (6) Eq. (7) 
0.0636 0.0121 0.0167 
0.00591 

199 2.06 2.83 
TABLE 


14 lg ~ 1Boy 
1Bsy 


Intensities calculated: 


Using Using 
screening assumption 
assumption zero 
With eqs. (10) and overlap, etc. 
eqs. (10) (7), and eqs. (10) and 
0.0065 0.040 
0.0040 0.0015 
1.6 2.9 
TABLE III 


fo 


f® by: 
Eq. (2) Eq. (6) 


The use incorrect co-ordinate transformation relate the cartesian 
symmetry co-ordinates benzene the appropriate internal symmetry 
co-ordinates invalidates the numerical results reported Liehr (1955, 1957, 
The corrected intensities are given Tables through readily 
apparent, although the quantitative results are somewhat improved, the 
previous qualitative conclusions remain unchanged. 


Experimental 
intensities 


0.094 
0.0014 


Liehr (1955) and 
Liehr and Moffitt 
(1958) using either 
the normal the 
Wilsonian sym- 
metry co-ordinates 


0.10 
~100 


Eq. (7) 


0.184 
4.43 


*Liehr (1957). 
the algebraic details the correction see Liehr Table this paper please 
read for Dys the number —0.8836. 
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TABLE 
by: 
(10) Liehr (1955) and 
and (7) (10) Liehr and Moffitt 
with and (7) with (1958) with internal 
Egs. screening zero overlap, symmetry 
(10) and (7) assumption etc. co-ordinates 
0.32 0.90 0.85 2.8 


Recently Pople and Sidman (1957) have utilized the Murrell and Pople 
(1956) semiempirical model transcribe Herzberg and Teller’s (1933) 
qualitative discussion the forbidden bands formaldehyde into quanti- 
tative form. The success their treatment quite spectacular. therefore 
appears that unduly pessimistic attitude was expressed the present 
author upon the basis his benzene intensity studies (Liehr 1957). 

illustration the sensitivity the calculations the numerical 
values the benzene electronic energy differences and the in- 
tensity utilized the computation, have presented Table the result 
using the experimental energy differences and intensity 
equations (2), (6), and (7) Liehr (1957). 


TABLE 
COMPARISON THE VARIOUS NON-EMPIRICAL INTENSITIES USING THE EXPERIMENTAL ENERGY 
DIFFERENCES (PARR, CRAIG, AND Ross 1950) AND THE EXPERIMENTAL 
INTENSITY (HAMMOND AND PRICE 1955) 


Intensities calc. basis of: 


Eq. (2) Eq. (6) Eq. (7) 
0.171 0.0480 0.0527 
1 1 
452 6.91 7.54 


14 


HERZBERG, and TELLER, 1933. Physik. Chem. 21, 
‘conjugated hydrocarbons. Thesis, Harvard University, Cambridge, Mass. (Corrected 
copies this thesis are available from the author upon request.) 
1957. Can. Phys. 35, 1123. 
Naturforsch 13a, 311. 
Naturforsch 13a, 596. 
Parr, G., P., and Ross, 1950. Chem. Phys. 18, 1561. 
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RADICAL SPECTRA LIQUID HELIUM TEMPERATURES* 


WILSE AND McCarty, Jr. 


Recently reported the electronic absorption spectra certain radicals 
trapped argon matrix (Robinson and McCarty 1958 
Figures and show the results some this preliminary work regard 
the HNO, OH, and radicals. The photographs were made the 
first order 2-meter, 15000-line per inch grating instrument. Many the 
observed absorption are surprisingly sharp, especially those 
which are only about wide. Because environmental perturbations 
the ground and excited electronic states, the spectra all the molecules 
studied exhibit ‘‘matrix upward 300 These shifts may 
either the red the blue with respect the gas phase bands. The 
spectrum trapped HNO clearly shows number fine structure features 
which are result interactions the molecular motion with the lattice. 
Besides the three strong bands which reported previously, three additional 
bands belonging HNO have been found the solid phase spectrum 
6445.9 6255.6 and 5899.6 These transitions involve the excited state 
vibrational levels 2v2, and 

Thirteen the sharp absorption lines have now been correlated with the 
flash photolysis spectrum NH2 (Ramsay 1957), and these account for all the 
lines strong and medium intensity expected the 4000-7000 region, 
well for few weak ones. Rotational structure the electronic spectrum 
trapped indicates that this molecule experiencing nearly free 
end-over-end rotation the matrix phase. Fine structure the 
somewhat more hindered motions these molecules, since there poorer 
correlation with the gas phase rotational levels. the relative strength 
the absorption from the first excited rotational level shows 
actual population which about 100-fold that expected for thermal 
equilibrium 4.2° The level metastable because the selection 
rule for exchange identical nuclei, but the rotational part this selection 
rule expected somewhat relaxed because electron and nuclear spin 
interactions with the molecular rotation. weak absorption line occurs 
30466 which attributed tentatively trapped radicals the 
metastable electronic state. The photographs show the presence 
emission lines associated with the *II absorption band trapped 
NH. The emission results from resonance fluorescence from the low-lying 
rotational levels the state. evidence has been found yet for 
HCO, but current efforts are expected yield fruitful 
results for least some these species. 


*Work supported part the OOR, U.S. Army under Project No. TB2-0001 Contract 
No. DA-36-034-OR D-2169. 
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Because the relatively unperturbed condition which small molecules 
can trapped liquid helium temperatures, the technique should prove 
very useful one for the study diatomics and simple polyatomics through- 
out the entire optical region the spectrum. Studies this type should 
especially valuable the vacuum ultraviolet and the infrared regions where 
little known about free radical spectra. 

Details concerning some our low temperature free radical work will 
published shortly another journal. 


Ramsay, soc. roy. sci. 18, 471. 
and McCarty, M., Chem. Phys. 28, 349. 
Chem. Phys. 28, 350. 
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